Mercurial > repos > public > sbplib

comparison +sbp/D4Variable.m @ 389:42c89b5eedc0 feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add borrowing constants for D2 operators in D4Variable
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 10 Jan 2017 17:31:28 +0100
parents 25b01643e438
children
comparison
equal deleted inserted replaced
388:288138d9586d 389:42c89b5eedc0
28 obj.x = linspace(x_l, x_r,m)'; 28 obj.x = linspace(x_l, x_r,m)';
29 29
30 if order == 2 30 if order == 2
31 [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... 31 [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
32 sbp.implementations.d4_variable_2(m, obj.h); 32 sbp.implementations.d4_variable_2(m, obj.h);
33 obj.borrowing.M.d1 = 0.4000;
33 obj.borrowing.N.S2 = 1.2500; 34 obj.borrowing.N.S2 = 1.2500;
34 obj.borrowing.N.S3 = 0.4000; 35 obj.borrowing.N.S3 = 0.4000;
35 elseif order == 4 36 elseif order == 4
36 [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... 37 [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
37 sbp.implementations.d4_variable_4(m, obj.h); 38 sbp.implementations.d4_variable_4(m, obj.h);
39 obj.borrowing.M.d1 = 0.2508;
38 obj.borrowing.N.S2 = 0.5055; 40 obj.borrowing.N.S2 = 0.5055;
39 obj.borrowing.N.S3 = 0.9290; 41 obj.borrowing.N.S3 = 0.9290;
40 elseif order == 6 42 elseif order == 6
41 [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... 43 [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
42 sbp.implementations.d4_variable_6(m, obj.h); 44 sbp.implementations.d4_variable_6(m, obj.h);
45 obj.borrowing.M.d1 = 0.1878;
43 obj.borrowing.N.S2 = 0.3259; 46 obj.borrowing.N.S2 = 0.3259;
44 obj.borrowing.N.S3 = 0.1580; 47 obj.borrowing.N.S3 = 0.1580;
45 else 48 else
46 error('Invalid operator order.'); 49 error('Invalid operator order.');
47 end 50 end