Mercurial > repos > public > sbplib
diff +sbp/D2.m @ 255:df3cc9c5dffc operator_remake
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 |
line wrap: on
line diff
--- a/+sbp/D2.m Wed Sep 07 14:49:21 2016 +0200 +++ b/+sbp/D2.m Wed Sep 07 15:54:41 2016 +0200 @@ -12,19 +12,19 @@ function obj = D2(m,h,order) if order == 2 - [H, HI, D1, D2, e_1, e_m, M,Q S_1, S_m] = sbp.ordinary2(m,h); + [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = sbp.ordinary2(m,h); obj.borrowing.M.S = 0.4000; elseif order == 4 - [H, HI, D1, D2, e_1, e_m, M,Q S_1, S_m] = sbp.ordinary4(m,h); + [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = sbp.ordinary4(m,h); obj.borrowing.M.S = 0.2508; elseif order == 6 - [H, HI, D1, D2, e_1, e_m, M,Q S_1, S_m] = sbp.ordinary6(m,h); + [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = sbp.ordinary6(m,h); obj.borrowing.M.S = 0.1878; elseif order == 8 - [H, HI, D1, D2, e_1, e_m, M,Q S_1, S_m] = sbp.ordinary8(m,h); + [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = sbp.ordinary8(m,h); obj.borrowing.M.S = 0.0015; elseif order == 10 - [H, HI, D1, D2, e_1, e_m, M,Q S_1, S_m] = sbp.ordinary10(m,h); + [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = sbp.ordinary10(m,h); obj.borrowing.M.S = 0.0351; else error('Invalid operator order %d.',order);