Mercurial > repos > public > sbplib

diff +sbp/D4Variable.m @ 389:42c89b5eedc0 feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add borrowing constants for D2 operators in D4Variable
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 10 Jan 2017 17:31:28 +0100
parents 25b01643e438
children
line wrap: on
line diff
--- a/+sbp/D4Variable.m	Mon Jan 09 11:20:28 2017 +0100
+++ b/+sbp/D4Variable.m	Tue Jan 10 17:31:28 2017 +0100
@@ -30,16 +30,19 @@
             if order == 2
                 [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
                     sbp.implementations.d4_variable_2(m, obj.h);
+                obj.borrowing.M.d1 = 0.4000;
                 obj.borrowing.N.S2 = 1.2500;
                 obj.borrowing.N.S3 = 0.4000;
             elseif order == 4
                 [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
                     sbp.implementations.d4_variable_4(m, obj.h);
+                obj.borrowing.M.d1 = 0.2508;
                 obj.borrowing.N.S2 = 0.5055;
                 obj.borrowing.N.S3 = 0.9290;
             elseif order == 6
                 [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
                     sbp.implementations.d4_variable_6(m, obj.h);
+                obj.borrowing.M.d1 = 0.1878;
                 obj.borrowing.N.S2 = 0.3259;
                 obj.borrowing.N.S3 = 0.1580;
             else