Mercurial > repos > public > sbplib
comparison +sbp/+implementations/d4_lonely_6_min_boundary_points.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 |
|---|---|
| 37 e_l(1) = 1; | 37 e_l(1) = 1; |
| 38 e_r = rot90(e_l, 2); | 38 e_r = rot90(e_l, 2); |
| 39 | 39 |
| 40 d1_l = sparse(m,1); | 40 d1_l = sparse(m,1); |
| 41 d1_l(1:5) = [-25/12, 4, -3, 4/3, -1/4]/h; | 41 d1_l(1:5) = [-25/12, 4, -3, 4/3, -1/4]/h; |
| 42 d1_r = -rot90(d1_l); | 42 d1_r = -rot90(d1_l, 2); |
| 43 | 43 |
| 44 d2_l = sparse(m,1); | 44 d2_l = sparse(m,1); |
| 45 d2_l(1:5) = [0.35e2/0.12e2 -0.26e2/0.3e1 0.19e2/0.2e1 -0.14e2/0.3e1 0.11e2/0.12e2;]/h^2; | 45 d2_l(1:5) = [0.35e2/0.12e2 -0.26e2/0.3e1 0.19e2/0.2e1 -0.14e2/0.3e1 0.11e2/0.12e2;]/h^2; |
| 46 d2_r = rot90(d2_l, 2); | 46 d2_r = rot90(d2_l, 2); |
| 47 | 47 |