Mercurial > repos > public > sbplib
diff +sbp/D2VariablePeriodic.m @ 1337:bf2554f1825d feature/D2_boundary_opt
Add periodic D1 and D2 operators for orders 8,10,12
| author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
|---|---|
| date | Fri, 13 May 2022 13:28:10 +0200 |
| parents | 5ccf6aaf6d6b |
| children |
line wrap: on
line diff
--- a/+sbp/D2VariablePeriodic.m Fri May 13 13:26:16 2022 +0200 +++ b/+sbp/D2VariablePeriodic.m Fri May 13 13:28:10 2022 +0200 @@ -28,6 +28,27 @@ switch order + case 12 + [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_12(m,obj.h); + obj.borrowing.M.d1 = obj.H(1, 1) / obj.h; % First element in H/h is borrowing also for M + obj.borrowing.R.delta_D = inf; + + case 10 + [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_10(m,obj.h); + obj.borrowing.M.d1 = obj.H(1, 1) / obj.h; % First element in H/h is borrowing also for M + obj.borrowing.R.delta_D = inf; + + case 8 + [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_8(m,obj.h); + obj.borrowing.M.d1 = obj.H(1, 1) / obj.h; % First element in H/h is borrowing also for M + obj.borrowing.R.delta_D = inf; + case 6 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... obj.e_r, obj.d1_l, obj.d1_r] = ...