Mercurial > repos > public > sbplib
comparison +sbp/+implementations/d4_lonely_4_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 |
---|---|
30 e_l(1) = 1; | 30 e_l(1) = 1; |
31 e_r = rot90(e_l, 2); | 31 e_r = rot90(e_l, 2); |
32 | 32 |
33 d1_l = sparse(m,1); | 33 d1_l = sparse(m,1); |
34 d1_l(1:4) = 1/h*[-11/6 3 -3/2 1/3]; | 34 d1_l(1:4) = 1/h*[-11/6 3 -3/2 1/3]; |
35 d1_r = -rot90(d1_l); | 35 d1_r = -rot90(d1_l, 2); |
36 | 36 |
37 d2_l = sparse(m,1); | 37 d2_l = sparse(m,1); |
38 d2_l(1:4) = 1/h^2*[2 -5 4 -1]; | 38 d2_l(1:4) = 1/h^2*[2 -5 4 -1]; |
39 d2_r = rot90(d2_l, 2); | 39 d2_r = rot90(d2_l, 2); |
40 | 40 |