Mercurial > repos > public > sbplib

comparison +sbp/D4Lonely.m @ 327:d24869abc7cd feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Calculated borrowin for D4Lonely.
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 27 Sep 2016 09:46:58 +0200
parents 72468bc9b63f
children 31d6698c1edf
comparison
equal deleted inserted replaced
326:b19e142fcae1 327:d24869abc7cd
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_lonely_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.6244;
39 % obj.borrowing.N.S3 = 0.9290; 39 obj.borrowing.N.S3 = 1.3961;
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);
43 obj.borrowing.N.S2 = 0.5055; 43 obj.borrowing.N.S2 = 0.5055;
44 obj.borrowing.N.S3 = 0.9290; 44 obj.borrowing.N.S3 = 0.9290;
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_lonely_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.2931;
53 % obj.borrowing.N.S3 = 0.1580; 53 obj.borrowing.N.S3 = 0.0807;
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_lonely_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.2842;
58 % obj.borrowing.N.S3 = 0.1580; 58 obj.borrowing.N.S3 = 0.0709;
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_lonely_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.3569;
63 % obj.borrowing.N.S3 = 0.1580; 63 obj.borrowing.N.S3 = 0.1908;
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);
67 obj.borrowing.N.S2 = 0.3259; 67 obj.borrowing.N.S2 = 0.3259;
68 obj.borrowing.N.S3 = 0.1580; 68 obj.borrowing.N.S3 = 0.1580;
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_lonely_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.2804;
77 % obj.borrowing.N.S3 = 0.1580; 77 obj.borrowing.N.S3 = 0.0740;
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_lonely_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.2475;
82 % obj.borrowing.N.S3 = 0.1580; 82 obj.borrowing.N.S3 = 0.0401;
83 end 83 end
84 else 84 else
85 error('Invalid operator order.'); 85 error('Invalid operator order.');
86 end 86 end
87 87