Mercurial > repos > public > sbplib

comparison +util/calc_borrowing.m @ 147:03151dd46534

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed function name bug in calc_borrowing and added ordinary operators.
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 04 Apr 2016 18:34:12 +0200
parents 48b6fb693025
children 8b4993d53663
comparison
equal deleted inserted replaced
146:2ccfe80e9b58 147:03151dd46534
30 fprintf('alpha_III: %.10f\n',alpha_III) 30 fprintf('alpha_III: %.10f\n',alpha_III)
31 fprintf('\n') 31 fprintf('\n')
32 32
33 33
34 %% 2nd order compatible 34 %% 2nd order compatible
35 [H, HI, D1, D4, e_1, e_m, M4, Q, S2_1, S2_m, S3_1, S3_m, S_1, S_m] = sbp.higher2_compatible(m,h); 35 [H, HI, D1, D4, e_1, e_m, M4, Q, S2_1, S2_m, S3_1, S3_m, S_1, S_m] = sbp.higher_compatible2(m,h);
36 S1 = S_1*S_1' + S_m*S_m'; 36 S1 = S_1*S_1' + S_m*S_m';
37 S2 = S2_1*S2_1' + S2_m*S2_m'; 37 S2 = S2_1*S2_1' + S2_m*S2_m';
38 S3 = S3_1*S3_1' + S3_m*S3_m'; 38 S3 = S3_1*S3_1' + S3_m*S3_m';
39 39
40 alpha_II = util.matrixborrow(M4, h*S2 ); 40 alpha_II = util.matrixborrow(M4, h*S2 );
44 fprintf('alpha_III: %.10f\n',alpha_III) 44 fprintf('alpha_III: %.10f\n',alpha_III)
45 fprintf('\n') 45 fprintf('\n')
46 46
47 47
48 %% 4th order compatible 48 %% 4th order compatible
49 [H, HI, D1, D4, e_1, e_m, M4, Q, S2_1, S2_m, S3_1, S3_m, S_1, S_m] = sbp.higher4_compatible(m,h); 49 [H, HI, D1, D4, e_1, e_m, M4, Q, S2_1, S2_m, S3_1, S3_m, S_1, S_m] = sbp.higher_compatible4(m,h);
50 S1 = S_1*S_1' + S_m*S_m'; 50 S1 = S_1*S_1' + S_m*S_m';
51 S2 = S2_1*S2_1' + S2_m*S2_m'; 51 S2 = S2_1*S2_1' + S2_m*S2_m';
52 S3 = S3_1*S3_1' + S3_m*S3_m'; 52 S3 = S3_1*S3_1' + S3_m*S3_m';
53 53
54 alpha_II = util.matrixborrow(M4, h*S2 ); 54 alpha_II = util.matrixborrow(M4, h*S2 );
57 fprintf('alpha_II: %.10f\n',alpha_II) 57 fprintf('alpha_II: %.10f\n',alpha_II)
58 fprintf('alpha_III: %.10f\n',alpha_III) 58 fprintf('alpha_III: %.10f\n',alpha_III)
59 fprintf('\n') 59 fprintf('\n')
60 60
61 %% 6th order compatible 61 %% 6th order compatible
62 [H, HI, D1, D4, e_1, e_m, M4, Q, S2_1, S2_m, S3_1, S3_m, S_1, S_m] = sbp.higher6_compatible(m,h); 62 [H, HI, D1, D4, e_1, e_m, M4, Q, S2_1, S2_m, S3_1, S3_m, S_1, S_m] = sbp.higher_compatible6(m,h);
63 S1 = S_1*S_1' + S_m*S_m'; 63 S1 = S_1*S_1' + S_m*S_m';
64 S2 = S2_1*S2_1' + S2_m*S2_m'; 64 S2 = S2_1*S2_1' + S2_m*S2_m';
65 S3 = S3_1*S3_1' + S3_m*S3_m'; 65 S3 = S3_1*S3_1' + S3_m*S3_m';
66 66
67 alpha_II = util.matrixborrow(M4, h*S2 ); 67 alpha_II = util.matrixborrow(M4, h*S2 );
68 alpha_III = util.matrixborrow(M4, h^3*S3); 68 alpha_III = util.matrixborrow(M4, h^3*S3);
69 fprintf('6th order compatible\n') 69 fprintf('6th order compatible\n')
70 fprintf('alpha_II: %.10f\n',alpha_II) 70 fprintf('alpha_II: %.10f\n',alpha_II)
71 fprintf('alpha_III: %.10f\n',alpha_III) 71 fprintf('alpha_III: %.10f\n',alpha_III)
72 fprintf('\n')3
73
74
75
76
77
78 % Ordinary
79
80 for order = [2 4 6 8 10]
81 op = sbp.Ordinary(m,h, order);
82
83 S_1 = op.boundary.S_1;
84 S_m = op.boundary.S_m;
85
86 M = op.norms.M;
87
88 S1 = S_1*S_1' + S_m*S_m';
89 alpha = util.matrixborrow(M, h*S1);
90 fprintf('%dth order Ordinary\n', order)
91 fprintf('alpha: %.10f\n', alpha)
92 fprintf('\n')
93 end
94
95