Mercurial > repos > public > sbplib

comparison +sbp/D2Variable.m @ 362:ded4156e53e2

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added 2nd order accurate 2nd derivative with variable coefficents in a separate implementation file, used by the class D2Variable.
author Martin Almquist <martin.almquist@it.uu.se>
date Tue, 20 Dec 2016 15:04:02 +0100
parents 4b9310edcdf8
children 359861563866
comparison
equal deleted inserted replaced
303:f18142c1530b 362:ded4156e53e2
29 case 4 29 case 4
30 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... 30 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,...
31 obj.e_r, obj.d1_l, obj.d1_r] = ... 31 obj.e_r, obj.d1_l, obj.d1_r] = ...
32 sbp.implementations.d2_variable_4(m,obj.h); 32 sbp.implementations.d2_variable_4(m,obj.h);
33 obj.borrowing.M.S = 0.2505765857; 33 obj.borrowing.M.S = 0.2505765857;
34 case 2
35 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,...
36 obj.e_r, obj.d1_l, obj.d1_r] = ...
37 sbp.implementations.d2_variable_2(m,obj.h);
38 obj.borrowing.M.S = 0.3636363636;
39 % Borrowing const taken from Virta 2014
40
34 otherwise 41 otherwise
35 error('Invalid operator order %d.',order); 42 error('Invalid operator order %d.',order);
36 end 43 end
37 44
38 obj.m = m; 45 obj.m = m;