Mercurial > repos > public > sbplib

comparison +util/calc_borrowing.m @ 0:48b6fb693025

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Initial commit.
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 17 Sep 2015 10:12:50 +0200
parents
children 03151dd46534
comparison
equal deleted inserted replaced
-1:000000000000 0:48b6fb693025
1
2 m = 30;
3 h = 1;
4
5
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);
8 S1 = S_1*S_1' + S_m*S_m';
9 S2 = S2_1*S2_1' + S2_m*S2_m';
10 S3 = S3_1*S3_1' + S3_m*S3_m';
11
12 alpha_II = util.matrixborrow(M4, h*S2 );
13 alpha_III = util.matrixborrow(M4, h^3*S3);
14 fprintf('4th order non-compatible\n')
15 fprintf('alpha_II: %.10f\n',alpha_II)
16 fprintf('alpha_III: %.10f\n',alpha_III)
17 fprintf('\n')
18
19
20 %% 6th order non-compatible
21 [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.higher6(m,h);
22 S1 = S_1*S_1' + S_m*S_m';
23 S2 = S2_1*S2_1' + S2_m*S2_m';
24 S3 = S3_1*S3_1' + S3_m*S3_m';
25
26 alpha_II = util.matrixborrow(M4, h*S2 );
27 alpha_III = util.matrixborrow(M4, h^3*S3);
28 fprintf('6th order non-compatible\n')
29 fprintf('alpha_II: %.10f\n',alpha_II)
30 fprintf('alpha_III: %.10f\n',alpha_III)
31 fprintf('\n')
32
33
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);
36 S1 = S_1*S_1' + S_m*S_m';
37 S2 = S2_1*S2_1' + S2_m*S2_m';
38 S3 = S3_1*S3_1' + S3_m*S3_m';
39
40 alpha_II = util.matrixborrow(M4, h*S2 );
41 alpha_III = util.matrixborrow(M4, h^3*S3);
42 fprintf('2nd order compatible\n')
43 fprintf('alpha_II: %.10f\n',alpha_II)
44 fprintf('alpha_III: %.10f\n',alpha_III)
45 fprintf('\n')
46
47
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);
50 S1 = S_1*S_1' + S_m*S_m';
51 S2 = S2_1*S2_1' + S2_m*S2_m';
52 S3 = S3_1*S3_1' + S3_m*S3_m';
53
54 alpha_II = util.matrixborrow(M4, h*S2 );
55 alpha_III = util.matrixborrow(M4, h^3*S3);
56 fprintf('4th order compatible\n')
57 fprintf('alpha_II: %.10f\n',alpha_II)
58 fprintf('alpha_III: %.10f\n',alpha_III)
59 fprintf('\n')
60
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);
63 S1 = S_1*S_1' + S_m*S_m';
64 S2 = S2_1*S2_1' + S2_m*S2_m';
65 S3 = S3_1*S3_1' + S3_m*S3_m';
66
67 alpha_II = util.matrixborrow(M4, h*S2 );
68 alpha_III = util.matrixborrow(M4, h^3*S3);
69 fprintf('6th order compatible\n')
70 fprintf('alpha_II: %.10f\n',alpha_II)
71 fprintf('alpha_III: %.10f\n',alpha_III)