Mercurial > repos > public > sbplib
diff +util/calc_borrowing.m @ 0:48b6fb693025
Initial commit.
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Thu, 17 Sep 2015 10:12:50 +0200 |
| parents | |
| children | 03151dd46534 |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/+util/calc_borrowing.m Thu Sep 17 10:12:50 2015 +0200 @@ -0,0 +1,71 @@ + +m = 30; +h = 1; + + +%% 4th order non-compatible +[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); +S1 = S_1*S_1' + S_m*S_m'; +S2 = S2_1*S2_1' + S2_m*S2_m'; +S3 = S3_1*S3_1' + S3_m*S3_m'; + +alpha_II = util.matrixborrow(M4, h*S2 ); +alpha_III = util.matrixborrow(M4, h^3*S3); +fprintf('4th order non-compatible\n') +fprintf('alpha_II: %.10f\n',alpha_II) +fprintf('alpha_III: %.10f\n',alpha_III) +fprintf('\n') + + +%% 6th order non-compatible +[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); +S1 = S_1*S_1' + S_m*S_m'; +S2 = S2_1*S2_1' + S2_m*S2_m'; +S3 = S3_1*S3_1' + S3_m*S3_m'; + +alpha_II = util.matrixborrow(M4, h*S2 ); +alpha_III = util.matrixborrow(M4, h^3*S3); +fprintf('6th order non-compatible\n') +fprintf('alpha_II: %.10f\n',alpha_II) +fprintf('alpha_III: %.10f\n',alpha_III) +fprintf('\n') + + +%% 2nd order compatible +[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); +S1 = S_1*S_1' + S_m*S_m'; +S2 = S2_1*S2_1' + S2_m*S2_m'; +S3 = S3_1*S3_1' + S3_m*S3_m'; + +alpha_II = util.matrixborrow(M4, h*S2 ); +alpha_III = util.matrixborrow(M4, h^3*S3); +fprintf('2nd order compatible\n') +fprintf('alpha_II: %.10f\n',alpha_II) +fprintf('alpha_III: %.10f\n',alpha_III) +fprintf('\n') + + +%% 4th order compatible +[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); +S1 = S_1*S_1' + S_m*S_m'; +S2 = S2_1*S2_1' + S2_m*S2_m'; +S3 = S3_1*S3_1' + S3_m*S3_m'; + +alpha_II = util.matrixborrow(M4, h*S2 ); +alpha_III = util.matrixborrow(M4, h^3*S3); +fprintf('4th order compatible\n') +fprintf('alpha_II: %.10f\n',alpha_II) +fprintf('alpha_III: %.10f\n',alpha_III) +fprintf('\n') + +%% 6th order compatible +[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); +S1 = S_1*S_1' + S_m*S_m'; +S2 = S2_1*S2_1' + S2_m*S2_m'; +S3 = S3_1*S3_1' + S3_m*S3_m'; + +alpha_II = util.matrixborrow(M4, h*S2 ); +alpha_III = util.matrixborrow(M4, h^3*S3); +fprintf('6th order compatible\n') +fprintf('alpha_II: %.10f\n',alpha_II) +fprintf('alpha_III: %.10f\n',alpha_III)