Mercurial > repos > public > sbplib

comparison +sbp/D2.m @ 255:df3cc9c5dffc operator_remake

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added ordinary 12th order accurate, with D1*D1 as 2nd derivative.
author Martin Almquist <martin.almquist@it.uu.se>
date Wed, 07 Sep 2016 15:54:41 +0200
parents 07fa0d6a05bb
children 87504db14c6d
comparison
equal deleted inserted replaced
254:9c97eb0d61bf 255:df3cc9c5dffc
10 10
11 methods 11 methods
12 function obj = D2(m,h,order) 12 function obj = D2(m,h,order)
13 13
14 if order == 2 14 if order == 2
15 [H, HI, D1, D2, e_1, e_m, M,Q S_1, S_m] = sbp.ordinary2(m,h); 15 [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = sbp.ordinary2(m,h);
16 obj.borrowing.M.S = 0.4000; 16 obj.borrowing.M.S = 0.4000;
17 elseif order == 4 17 elseif order == 4
18 [H, HI, D1, D2, e_1, e_m, M,Q S_1, S_m] = sbp.ordinary4(m,h); 18 [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = sbp.ordinary4(m,h);
19 obj.borrowing.M.S = 0.2508; 19 obj.borrowing.M.S = 0.2508;
20 elseif order == 6 20 elseif order == 6
21 [H, HI, D1, D2, e_1, e_m, M,Q S_1, S_m] = sbp.ordinary6(m,h); 21 [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = sbp.ordinary6(m,h);
22 obj.borrowing.M.S = 0.1878; 22 obj.borrowing.M.S = 0.1878;
23 elseif order == 8 23 elseif order == 8
24 [H, HI, D1, D2, e_1, e_m, M,Q S_1, S_m] = sbp.ordinary8(m,h); 24 [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = sbp.ordinary8(m,h);
25 obj.borrowing.M.S = 0.0015; 25 obj.borrowing.M.S = 0.0015;
26 elseif order == 10 26 elseif order == 10
27 [H, HI, D1, D2, e_1, e_m, M,Q S_1, S_m] = sbp.ordinary10(m,h); 27 [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = sbp.ordinary10(m,h);
28 obj.borrowing.M.S = 0.0351; 28 obj.borrowing.M.S = 0.0351;
29 else 29 else
30 error('Invalid operator order %d.',order); 30 error('Invalid operator order %d.',order);
31 end 31 end
32 32