Mercurial > repos > public > sbplib

annotate +util/calc_borrowing.m @ 148:1f5f988e9682

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added functions for printing and timing tasks in a script.
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 05 Apr 2016 11:26:04 +0200
parents 03151dd46534
children 8b4993d53663
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
1
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
2 m = 30;
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
3 h = 1;
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
4
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
5
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
6 %% 4th order non-compatible
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
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);
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
8 S1 = S_1*S_1' + S_m*S_m';
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
9 S2 = S2_1*S2_1' + S2_m*S2_m';
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
10 S3 = S3_1*S3_1' + S3_m*S3_m';
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
11
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
12 alpha_II = util.matrixborrow(M4, h*S2 );
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
13 alpha_III = util.matrixborrow(M4, h^3*S3);
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
14 fprintf('4th order non-compatible\n')
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
15 fprintf('alpha_II: %.10f\n',alpha_II)
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
16 fprintf('alpha_III: %.10f\n',alpha_III)
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
17 fprintf('\n')
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
18
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
19
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
20 %% 6th order non-compatible
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
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);
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
22 S1 = S_1*S_1' + S_m*S_m';
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
23 S2 = S2_1*S2_1' + S2_m*S2_m';
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
24 S3 = S3_1*S3_1' + S3_m*S3_m';
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
25
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
26 alpha_II = util.matrixborrow(M4, h*S2 );
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
27 alpha_III = util.matrixborrow(M4, h^3*S3);
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
28 fprintf('6th order non-compatible\n')
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
29 fprintf('alpha_II: %.10f\n',alpha_II)
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
30 fprintf('alpha_III: %.10f\n',alpha_III)
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
31 fprintf('\n')
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
32
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
33
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
34 %% 2nd order compatible
147
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
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);
0
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
36 S1 = S_1*S_1' + S_m*S_m';
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
37 S2 = S2_1*S2_1' + S2_m*S2_m';
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
38 S3 = S3_1*S3_1' + S3_m*S3_m';
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
39
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
40 alpha_II = util.matrixborrow(M4, h*S2 );
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
41 alpha_III = util.matrixborrow(M4, h^3*S3);
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
42 fprintf('2nd order compatible\n')
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
43 fprintf('alpha_II: %.10f\n',alpha_II)
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
44 fprintf('alpha_III: %.10f\n',alpha_III)
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
45 fprintf('\n')
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
46
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
47
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
48 %% 4th order compatible
147
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
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);
0
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
50 S1 = S_1*S_1' + S_m*S_m';
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
51 S2 = S2_1*S2_1' + S2_m*S2_m';
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
52 S3 = S3_1*S3_1' + S3_m*S3_m';
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
53
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
54 alpha_II = util.matrixborrow(M4, h*S2 );
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
55 alpha_III = util.matrixborrow(M4, h^3*S3);
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
56 fprintf('4th order compatible\n')
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
57 fprintf('alpha_II: %.10f\n',alpha_II)
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
58 fprintf('alpha_III: %.10f\n',alpha_III)
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
59 fprintf('\n')
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
60
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
61 %% 6th order compatible
147
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
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);
0
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
63 S1 = S_1*S_1' + S_m*S_m';
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
64 S2 = S2_1*S2_1' + S2_m*S2_m';
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
65 S3 = S3_1*S3_1' + S3_m*S3_m';
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
66
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
67 alpha_II = util.matrixborrow(M4, h*S2 );
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
68 alpha_III = util.matrixborrow(M4, h^3*S3);
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
69 fprintf('6th order compatible\n')
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
70 fprintf('alpha_II: %.10f\n',alpha_II)
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
71 fprintf('alpha_III: %.10f\n',alpha_III)
147
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
72 fprintf('\n')3
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
73
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
74
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
75
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
76
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
77
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
78 % Ordinary
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
79
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
80 for order = [2 4 6 8 10]
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
81 op = sbp.Ordinary(m,h, order);
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
82
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
83 S_1 = op.boundary.S_1;
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
84 S_m = op.boundary.S_m;
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
85
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
86 M = op.norms.M;
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
87
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
88 S1 = S_1*S_1' + S_m*S_m';
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
89 alpha = util.matrixborrow(M, h*S1);
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
90 fprintf('%dth order Ordinary\n', order)
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
91 fprintf('alpha: %.10f\n', alpha)
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
92 fprintf('\n')
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
93 end
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
94
03151dd46534 Fixed function name bug in calc_borrowing and added ordinary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 0
diff changeset
95