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