Mercurial > repos > public > sbplib
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 |