Mercurial > repos > public > sbplib
diff +sbp/D2Variable.m @ 1140:738de3a4058b feature/laplace_curvilinear_test
Update borrowing constants so that M constants are from Virta2014 and R constants are for bp=2,4,7.
| author | Martin Almquist <malmquist@stanford.edu> |
|---|---|
| date | Fri, 21 Jun 2019 17:28:02 +0200 |
| parents | 52a9dedb9171 |
| children |
line wrap: on
line diff
--- a/+sbp/D2Variable.m Fri Jun 21 16:27:49 2019 +0200 +++ b/+sbp/D2Variable.m Fri Jun 21 17:28:02 2019 +0200 @@ -33,28 +33,34 @@ ~, obj.e_l, obj.e_r, ~, ~, ~, ~, ~,... obj.d1_l, obj.d1_r] = ... sbp.implementations.d4_variable_6(m, obj.h); - obj.borrowing.M.d1 = 0.1878; - obj.borrowing.R.delta_D = 0.3696; + % obj.borrowing.M.d1 = 0.1878; + + % From Virta 2014: + obj.borrowing.M.d1 = 0.1878687080; + + obj.borrowing.R.delta_D = 0.369684816368150; % Borrowing e^T*D1 - d1 from R case 4 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... obj.e_r, obj.d1_l, obj.d1_r] = ... sbp.implementations.d2_variable_4(m,obj.h); + + % From Virta 2014: obj.borrowing.M.d1 = 0.2505765857; - obj.borrowing.R.delta_D = 0.577587500088313; + obj.borrowing.R.delta_D = 0.577587500088252; % Borrowing e^T*D1 - d1 from R case 2 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... obj.e_r, obj.d1_l, obj.d1_r] = ... sbp.implementations.d2_variable_2(m,obj.h); - obj.borrowing.M.d1 = 0.3636363636; + obj.borrowing.M.d1 = 0.3636363636; % Borrowing const taken from Virta 2014 - obj.borrowing.R.delta_D = 1.000000538455350; + obj.borrowing.R.delta_D = 1.000000000000000; % Borrowing e^T*D1 - d1 from R - + otherwise error('Invalid operator order %d.',order); end