Mercurial > repos > public > sbplib

comparison +util/calc_borrowing.m @ 219:f66513508c75 feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merged with feature/grids.
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 27 Jun 2016 13:26:02 +0200
parents 8b4993d53663
children d24869abc7cd
comparison
equal deleted inserted replaced
217:fc07ebc49412 219:f66513508c75
1 1
2 m = 30; 2 m = 100;
3 h = 1; 3 h = 1;
4 4
5 5
6 %% 4th order non-compatible 6 %% 4th order non-compatible
7 [H, HI, D1, D2, D3, D4, e_1, e_m, M, M4,Q, Q3, S2_1, S2_m, S3_1, S3_m, S_1, S_m] = sbp.higher4(m,h); 7 [H, HI, D1, D2, D3, D4, e_1, e_m, M, M4,Q, Q3, S2_1, S2_m, S3_1, S3_m, S_1, S_m] = sbp.higher4(m,h);
8 S1 = S_1*S_1' + S_m*S_m'; 8 S1 = S_1*S_1' + S_m*S_m';
9 S2 = S2_1*S2_1' + S2_m*S2_m'; 9 S2 = S2_1*S2_1' + S2_m*S2_m';
10 S3 = S3_1*S3_1' + S3_m*S3_m'; 10 S3 = S3_1*S3_1' + S3_m*S3_m';
11 11
12 alpha_I = util.matrixborrow(M4, h^-1*S1 );
12 alpha_II = util.matrixborrow(M4, h*S2 ); 13 alpha_II = util.matrixborrow(M4, h*S2 );
13 alpha_III = util.matrixborrow(M4, h^3*S3); 14 alpha_III = util.matrixborrow(M4, h^3*S3);
14 fprintf('4th order non-compatible\n') 15 fprintf('4th order non-compatible\n')
16 fprintf('alpha_I1: %.10f\n',alpha_I)
15 fprintf('alpha_II: %.10f\n',alpha_II) 17 fprintf('alpha_II: %.10f\n',alpha_II)
16 fprintf('alpha_III: %.10f\n',alpha_III) 18 fprintf('alpha_III: %.10f\n',alpha_III)
17 fprintf('\n') 19 fprintf('\n')
18 20
19 21
67 alpha_II = util.matrixborrow(M4, h*S2 ); 69 alpha_II = util.matrixborrow(M4, h*S2 );
68 alpha_III = util.matrixborrow(M4, h^3*S3); 70 alpha_III = util.matrixborrow(M4, h^3*S3);
69 fprintf('6th order compatible\n') 71 fprintf('6th order compatible\n')
70 fprintf('alpha_II: %.10f\n',alpha_II) 72 fprintf('alpha_II: %.10f\n',alpha_II)
71 fprintf('alpha_III: %.10f\n',alpha_III) 73 fprintf('alpha_III: %.10f\n',alpha_III)
72 fprintf('\n')3 74 fprintf('\n')
73 75
74 76
75 77
76 78
77 79