Mercurial > repos > public > sbplib
diff +sbp/D4Standard.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 | 8a625c5a3633 |
| children | e1d11b6a68d8 |
line wrap: on
line diff
--- a/+sbp/D4Standard.m Fri Sep 09 14:53:41 2016 +0200 +++ b/+sbp/D4Standard.m Fri Sep 09 15:06:03 2016 +0200 @@ -4,20 +4,20 @@ 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 D3 % SBP operator for third derivative Q3 % Skew-symmetric matrix in third derivative - S2_1 % Left boundary second derivative - S2_m % Right boundary second derivative + d2_l % Left boundary second derivative + d2_r % Right boundary second derivative D4 % SBP operator for fourth derivative M4 % Norm matrix, fourth derivative - S3_1 % Left boundary third derivative - S3_m % Right boundary third derivative + d3_l % Left boundary third derivative + d3_r % Right boundary third derivative m % Number of grid points. h % Step size x % grid @@ -37,15 +37,15 @@ if order == 4 [obj.H, obj.HI, obj.D1, obj.D2, obj.D3, obj.D4,... - obj.e_1, obj.e_m, obj.M, obj.M4, obj.Q, obj.Q3, obj.S2_1,... - obj.S2_m, obj.S3_1, obj.S3_m, obj.S_1, obj.S_m] = ... + obj.e_l, obj.e_r, obj.M, obj.M4, obj.Q, obj.Q3, obj.d2_l,... + obj.d2_r, obj.d3_l, obj.d3_r, obj.d1_l, obj.d1_r] = ... sbp.implementations.d4_4(m,obj.h); obj.borrowing.N.S2 = 0.5485; obj.borrowing.N.S3 = 1.0882; elseif order == 6 [obj.H, obj.HI, obj.D1, obj.D2, obj.D3, obj.D4,... - obj.e_1, obj.e_m, obj.M, obj.M4, obj.Q, obj.Q3, obj.S2_1,... - obj.S2_m, obj.S3_1, obj.S3_m, obj.S_1, obj.S_m] = ... + obj.e_l, obj.e_r, obj.M, obj.M4, obj.Q, obj.Q3, obj.d2_l,... + obj.d2_r, obj.d3_l, obj.d3_r, obj.d1_l, obj.d1_r] = ... sbp.implementations.d4_6(m,obj.h); obj.borrowing.N.S2 = 0.3227; obj.borrowing.N.S3 = 0.1568;