Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d4_variable_6.m @ 326:b19e142fcae1 feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed bug in setting of boundary derivative.
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 27 Sep 2016 08:44:17 +0200
parents c0cbffcf6513
children bf801c3709be
comparison
equal deleted inserted replaced
325:72468bc9b63f 326:b19e142fcae1
35 e_l(1) = 1; 35 e_l(1) = 1;
36 e_r = rot90(e_l, 2); 36 e_r = rot90(e_l, 2);
37 37
38 d1_l = sparse(m,1); 38 d1_l = sparse(m,1);
39 d1_l(1:5) = [-25/12, 4, -3, 4/3, -1/4]/h; 39 d1_l(1:5) = [-25/12, 4, -3, 4/3, -1/4]/h;
40 d1_r = -rot90(d1_l); 40 d1_r = -rot90(d1_l, 2);
41 41
42 d2_l = sparse(m,1); 42 d2_l = sparse(m,1);
43 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; 43 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;
44 d2_r = rot90(d2_l, 2); 44 d2_r = rot90(d2_l, 2);
45 45