Mercurial > repos > public > sbplib
comparison +sbp/+implementations/d4_lonely_6_3.m @ 326:b19e142fcae1 feature/beams
Fixed bug in setting of boundary derivative.
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Tue, 27 Sep 2016 08:44:17 +0200 |
| parents | 72468bc9b63f |
| children |
comparison
equal
deleted
inserted
replaced
| 325:72468bc9b63f | 326:b19e142fcae1 |
|---|---|
| 34 e_l(1) = 1; | 34 e_l(1) = 1; |
| 35 e_r = rot90(e_l, 2); | 35 e_r = rot90(e_l, 2); |
| 36 | 36 |
| 37 d1_l = sparse(m,1); | 37 d1_l = sparse(m,1); |
| 38 d1_l(1:6) = [-0.137e3/0.60e2 5 -5 0.10e2/0.3e1 -0.5e1/0.4e1 0.1e1/0.5e1;]/h; | 38 d1_l(1:6) = [-0.137e3/0.60e2 5 -5 0.10e2/0.3e1 -0.5e1/0.4e1 0.1e1/0.5e1;]/h; |
| 39 d1_r = -rot90(d1_l); | 39 d1_r = -rot90(d1_l, 2); |
| 40 | 40 |
| 41 d2_l = sparse(m,1); | 41 d2_l = sparse(m,1); |
| 42 d2_l(1:6) = [0.15e2/0.4e1 -0.77e2/0.6e1 0.107e3/0.6e1 -13 0.61e2/0.12e2 -0.5e1/0.6e1;]/h^2; | 42 d2_l(1:6) = [0.15e2/0.4e1 -0.77e2/0.6e1 0.107e3/0.6e1 -13 0.61e2/0.12e2 -0.5e1/0.6e1;]/h^2; |
| 43 d2_r = rot90(d2_l, 2); | 43 d2_r = rot90(d2_l, 2); |
| 44 | 44 |