Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d4_lonely_8_higher_boundary_order.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 72468bc9b63f
children bf801c3709be
comparison
equal deleted inserted replaced
325:72468bc9b63f 326:b19e142fcae1
36 e_l(1) = 1; 36 e_l(1) = 1;
37 e_r = rot90(e_l, 2); 37 e_r = rot90(e_l, 2);
38 38
39 d1_l = sparse(m,1); 39 d1_l = sparse(m,1);
40 d1_l(1:7) = [-0.49e2/0.20e2 6 -0.15e2/0.2e1 0.20e2/0.3e1 -0.15e2/0.4e1 0.6e1/0.5e1 -0.1e1/0.6e1]/h; 40 d1_l(1:7) = [-0.49e2/0.20e2 6 -0.15e2/0.2e1 0.20e2/0.3e1 -0.15e2/0.4e1 0.6e1/0.5e1 -0.1e1/0.6e1]/h;
41 d1_r = -rot90(d1_l); 41 d1_r = -rot90(d1_l, 2);
42 42
43 d2_l = sparse(m,1); 43 d2_l = sparse(m,1);
44 d2_l(1:7) = [0.203e3/0.45e2 -0.87e2/0.5e1 0.117e3/0.4e1 -0.254e3/0.9e1 0.33e2/0.2e1 -0.27e2/0.5e1 0.137e3/0.180e3]/h^2; 44 d2_l(1:7) = [0.203e3/0.45e2 -0.87e2/0.5e1 0.117e3/0.4e1 -0.254e3/0.9e1 0.33e2/0.2e1 -0.27e2/0.5e1 0.137e3/0.180e3]/h^2;
45 d2_r = rot90(d2_l, 2); 45 d2_r = rot90(d2_l, 2);
46 46