Mercurial > repos > public > sbplib

comparison +sbp/D2.m @ 261:6009f2712d13 operator_remake

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Moved and renamned all implementations.
author Martin Almquist <martin.almquist@it.uu.se>
date Thu, 08 Sep 2016 15:35:45 +0200
parents 87504db14c6d
children 23051a86faa4
comparison
equal deleted inserted replaced
260:b4116ce49ac4 261:6009f2712d13
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.implementations.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.implementations.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.implementations.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.implementations.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.implementations.ordinary10(m,h);
28 obj.borrowing.M.S = 0.0351; 28 obj.borrowing.M.S = 0.0351;
29 elseif order == 12 29 elseif order == 12
30 [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = sbp.ordinary12(m,h); 30 [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = sbp.implementations.ordinary12(m,h);
31 obj.borrowing.M.S = []; % Need to compute! 31 obj.borrowing.M.S = []; % Need to compute!
32 else 32 else
33 error('Invalid operator order %d.',order); 33 error('Invalid operator order %d.',order);
34 end 34 end
35 35