Mercurial > repos > public > sbplib

comparison +sbp/D4Lonely.m @ 325:72468bc9b63f feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Renamed some operator implementations.
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 26 Sep 2016 09:55:16 +0200
parents c0cbffcf6513
children d24869abc7cd
comparison
equal deleted inserted replaced
324:c0cbffcf6513 325:72468bc9b63f
32 32
33 elseif order == 4 33 elseif order == 4
34 switch opt 34 switch opt
35 case 'min_boundary_points' 35 case 'min_boundary_points'
36 [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... 36 [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
37 sbp.implementations.d4_variable_4_min_boundary_points(m, obj.h); 37 sbp.implementations.d4_lonely_4_min_boundary_points(m, obj.h);
38 % obj.borrowing.N.S2 = 0.5055; 38 % obj.borrowing.N.S2 = 0.5055;
39 % obj.borrowing.N.S3 = 0.9290; 39 % obj.borrowing.N.S3 = 0.9290;
40 otherwise 40 otherwise
41 [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... 41 [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_4(m, obj.h); 42 sbp.implementations.d4_variable_4(m, obj.h);
46 46
47 elseif order == 6 47 elseif order == 6
48 switch opt 48 switch opt
49 case '2' 49 case '2'
50 [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... 50 [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
51 sbp.implementations.d4_variable_6_2(m, obj.h); 51 sbp.implementations.d4_lonely_6_2(m, obj.h);
52 % obj.borrowing.N.S2 = 0.3259; 52 % obj.borrowing.N.S2 = 0.3259;
53 % obj.borrowing.N.S3 = 0.1580; 53 % obj.borrowing.N.S3 = 0.1580;
54 case '3' 54 case '3'
55 [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... 55 [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
56 sbp.implementations.d4_variable_6_3(m, obj.h); 56 sbp.implementations.d4_lonely_6_3(m, obj.h);
57 % obj.borrowing.N.S2 = 0.3259; 57 % obj.borrowing.N.S2 = 0.3259;
58 % obj.borrowing.N.S3 = 0.1580; 58 % obj.borrowing.N.S3 = 0.1580;
59 case 'min_boundary_points' 59 case 'min_boundary_points'
60 [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... 60 [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
61 sbp.implementations.d4_variable_6_min_boundary_points(m, obj.h); 61 sbp.implementations.d4_lonely_6_min_boundary_points(m, obj.h);
62 % obj.borrowing.N.S2 = 0.3259; 62 % obj.borrowing.N.S2 = 0.3259;
63 % obj.borrowing.N.S3 = 0.1580; 63 % obj.borrowing.N.S3 = 0.1580;
64 otherwise 64 otherwise
65 [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... 65 [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
66 sbp.implementations.d4_variable_6(m, obj.h); 66 sbp.implementations.d4_variable_6(m, obj.h);
70 70
71 elseif order == 8 71 elseif order == 8
72 switch opt 72 switch opt
73 case 'min_boundary_points' 73 case 'min_boundary_points'
74 [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... 74 [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
75 sbp.implementations.d4_variable_8_min_boundary_points(m, obj.h); 75 sbp.implementations.d4_lonely_8_min_boundary_points(m, obj.h);
76 % obj.borrowing.N.S2 = 0.3259; 76 % obj.borrowing.N.S2 = 0.3259;
77 % obj.borrowing.N.S3 = 0.1580; 77 % obj.borrowing.N.S3 = 0.1580;
78 otherwise 78 otherwise
79 [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... 79 [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
80 sbp.implementations.d4_variable_8_higher_boundary_order(m, obj.h); 80 sbp.implementations.d4_lonely_8_higher_boundary_order(m, obj.h);
81 % obj.borrowing.N.S2 = 0.3259; 81 % obj.borrowing.N.S2 = 0.3259;
82 % obj.borrowing.N.S3 = 0.1580; 82 % obj.borrowing.N.S3 = 0.1580;
83 end 83 end
84 else 84 else
85 error('Invalid operator order.'); 85 error('Invalid operator order.');