Mercurial > repos > public > sbplib
diff +sbp/D2Standard.m @ 268:4b9310edcdf8 operator_remake
Renamned boundary operators!
| author | Martin Almquist <martin.almquist@it.uu.se> |
|---|---|
| date | Fri, 09 Sep 2016 15:06:03 +0200 |
| parents | 3e0432b3cb68 |
| children | 479194747f9c |
line wrap: on
line diff
--- a/+sbp/D2Standard.m Fri Sep 09 14:53:41 2016 +0200 +++ b/+sbp/D2Standard.m Fri Sep 09 15:06:03 2016 +0200 @@ -4,12 +4,12 @@ H % Norm matrix HI % H^-1 Q % Skew-symmetric matrix - e_1 % Left boundary operator - e_m % Right boundary operator + e_l % Left boundary operator + e_r % Right boundary operator D2 % SBP operator for second derivative M % Norm matrix, second derivative - S_1 % Left boundary first derivative - S_m % Right boundary first derivative + d1_l % Left boundary first derivative + d1_r % Right boundary first derivative m % Number of grid points. h % Step size x % grid @@ -27,33 +27,33 @@ obj.x = linspace(x_l,x_r,m)'; if order == 2 - [obj.H, obj.HI, obj.D1, obj.D2, obj.e_1,... - obj.e_m, obj.M, obj.Q, obj.S_1, obj.S_m] =... + [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... + obj.e_r, obj.M, obj.Q, obj.d1_l, obj.d1_r] =... sbp.implementations.ordinary2(m,obj.h); obj.borrowing.M.S = 0.4000; elseif order == 4 - [obj.H, obj.HI, obj.D1, obj.D2, obj.e_1,... - obj.e_m, obj.M, obj.Q, obj.S_1, obj.S_m] =... + [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... + obj.e_r, obj.M, obj.Q, obj.d1_l, obj.d1_r] =... sbp.implementations.ordinary4(m,obj.h); obj.borrowing.M.S = 0.2508; elseif order == 6 - [obj.H, obj.HI, obj.D1, obj.D2, obj.e_1,... - obj.e_m, obj.M, obj.Q, obj.S_1, obj.S_m] =... + [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... + obj.e_r, obj.M, obj.Q, obj.d1_l, obj.d1_r] =... sbp.implementations.ordinary6(m,obj.h); obj.borrowing.M.S = 0.1878; elseif order == 8 - [obj.H, obj.HI, obj.D1, obj.D2, obj.e_1,... - obj.e_m, obj.M, obj.Q, obj.S_1, obj.S_m] =... + [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... + obj.e_r, obj.M, obj.Q, obj.d1_l, obj.d1_r] =... sbp.implementations.ordinary8(m,obj.h); obj.borrowing.M.S = 0.0015; elseif order == 10 - [obj.H, obj.HI, obj.D1, obj.D2, obj.e_1,... - obj.e_m, obj.M, obj.Q, obj.S_1, obj.S_m] =... + [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... + obj.e_r, obj.M, obj.Q, obj.d1_l, obj.d1_r] =... sbp.implementations.ordinary10(m,obj.h); obj.borrowing.M.S = 0.0351; elseif order == 12 - [obj.H, obj.HI, obj.D1, obj.D2, obj.e_1,... - obj.e_m, obj.M, obj.Q, obj.S_1, obj.S_m] =... + [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... + obj.e_r, obj.M, obj.Q, obj.d1_l, obj.d1_r] =... sbp.implementations.ordinary12(m,obj.h); obj.borrowing.M.S = []; % Need to compute! warning('D2=D1*D1 (wide stencil).');