Mercurial > repos > public > sbplib
diff +sbp/D2Variable.m @ 1071:92cb03e64ca4 feature/grids/LaplaceSquared
Merge with default
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Thu, 20 Sep 2018 12:05:20 +0200 |
| parents | 52a9dedb9171 |
| children | 738de3a4058b |
line wrap: on
line diff
--- a/+sbp/D2Variable.m Tue Jul 24 21:08:08 2018 -0700 +++ b/+sbp/D2Variable.m Thu Sep 20 12:05:20 2018 +0200 @@ -26,22 +26,39 @@ obj.x = linspace(x_l,x_r,m)'; switch order + + case 6 + + [obj.H, obj.HI, obj.D1, obj.D2, ... + ~, 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; + % 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); obj.borrowing.M.d1 = 0.2505765857; + + obj.borrowing.R.delta_D = 0.577587500088313; + % 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; % Borrowing const taken from Virta 2014 + + obj.borrowing.R.delta_D = 1.000000538455350; + % Borrowing e^T*D1 - d1 from R otherwise error('Invalid operator order %d.',order); end - + obj.borrowing.H11 = obj.H(1,1)/obj.h; % First element in H/h, obj.m = m; obj.M = []; end