comparison +sbp/+implementations/d4_variable_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 88584b0cfba1
children bf801c3709be
comparison
equal deleted inserted replaced
325:72468bc9b63f 326:b19e142fcae1
27 e_l(1) = 1; 27 e_l(1) = 1;
28 e_r = rot90(e_l, 2); 28 e_r = rot90(e_l, 2);
29 29
30 d1_l = sparse(m,1); 30 d1_l = sparse(m,1);
31 d1_l(1:3) = 1/h*[-3/2 2 -1/2]; 31 d1_l(1:3) = 1/h*[-3/2 2 -1/2];
32 d1_r = -rot90(d1_l); 32 d1_r = -rot90(d1_l, 2);
33 33
34 d2_l = sparse(m,1); 34 d2_l = sparse(m,1);
35 d2_l(1:3) = 1/h^2*[1 -2 1]; 35 d2_l(1:3) = 1/h^2*[1 -2 1];
36 d2_r = rot90(d2_l, 2); 36 d2_r = rot90(d2_l, 2);
37 37