Mercurial > repos > public > sbplib
comparison +sbp/+implementations/d4_lonely_6_2.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 |
|---|---|
| 39 e_l(1) = 1; | 39 e_l(1) = 1; |
| 40 e_r = rot90(e_l, 2); | 40 e_r = rot90(e_l, 2); |
| 41 | 41 |
| 42 d1_l = sparse(m,1); | 42 d1_l = sparse(m,1); |
| 43 d1_l(1:6) = [-0.137e3/0.60e2 5 -5 0.10e2/0.3e1 -0.5e1/0.4e1 0.1e1/0.5e1;]/h; | 43 d1_l(1:6) = [-0.137e3/0.60e2 5 -5 0.10e2/0.3e1 -0.5e1/0.4e1 0.1e1/0.5e1;]/h; |
| 44 d1_r = -rot90(d1_l); | 44 d1_r = -rot90(d1_l, 2); |
| 45 | 45 |
| 46 d2_l = sparse(m,1); | 46 d2_l = sparse(m,1); |
| 47 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; | 47 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; |
| 48 d2_r = rot90(d2_l, 2); | 48 d2_r = rot90(d2_l, 2); |
| 49 | 49 |