Mercurial > repos > public > sbplib
diff +sbp/D2VariablePeriodic.m @ 686:5ccf6aaf6d6b feature/poroelastic
Add D2VariablePeriodic orders 4 and 6.
| author | Martin Almquist <malmquist@stanford.edu> |
|---|---|
| date | Thu, 08 Feb 2018 16:44:46 -0800 |
| parents | 7368affc8f78 |
| children | bf2554f1825d |
line wrap: on
line diff
--- a/+sbp/D2VariablePeriodic.m Thu Feb 08 16:43:43 2018 -0800 +++ b/+sbp/D2VariablePeriodic.m Thu Feb 08 16:44:46 2018 -0800 @@ -29,19 +29,14 @@ switch order case 6 - error('Not impl') - - [obj.H, obj.HI, obj.D1, obj.D2, ... - ~, obj.e_l, obj.e_r, ~, ~, ~, ~, ~,... - obj.d1_l, obj.d1_r] = ... - sbp.implementations.d4_variable_periodic_6(m, obj.h); + [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... + obj.e_r, obj.d1_l, obj.d1_r] = ... + sbp.implementations.d2_variable_periodic_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 - error('Not impl') - [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... obj.e_r, obj.d1_l, obj.d1_r] = ... sbp.implementations.d2_variable_periodic_4(m,obj.h); @@ -63,6 +58,7 @@ 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