Mercurial > repos > public > sbplib
diff +sbp/D1Upwind.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 | f39f98b59f61 |
line wrap: on
line diff
--- a/+sbp/D1Upwind.m Fri Sep 09 14:53:41 2016 +0200 +++ b/+sbp/D1Upwind.m Fri Sep 09 15:06:03 2016 +0200 @@ -4,8 +4,8 @@ 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 m % Number of grid points. h % Step size x % grid @@ -23,28 +23,28 @@ switch order case 2 - [obj.H, obj.HI, obj.Dp, obj.Dm, obj.e_1, obj.e_m] = ... + [obj.H, obj.HI, obj.Dp, obj.Dm, obj.e_l, obj.e_r] = ... sbp.implementations.d1_upwind_2(m,obj.h); case 3 - [obj.H, obj.HI, obj.Dp, obj.Dm, obj.e_1, obj.e_m] = ... + [obj.H, obj.HI, obj.Dp, obj.Dm, obj.e_l, obj.e_r] = ... sbp.implementations.d1_upwind_3(m,obj.h); case 4 - [obj.H, obj.HI, obj.Dp, obj.Dm, obj.e_1, obj.e_m] = ... + [obj.H, obj.HI, obj.Dp, obj.Dm, obj.e_l, obj.e_r] = ... sbp.implementations.d1_upwind_4(m,obj.h); case 5 - [obj.H, obj.HI, obj.Dp, obj.Dm, obj.e_1, obj.e_m] = ... + [obj.H, obj.HI, obj.Dp, obj.Dm, obj.e_l, obj.e_r] = ... sbp.implementations.d1_upwind_5(m,obj.h); case 6 - [obj.H, obj.HI, obj.Dp, obj.Dm, obj.e_1, obj.e_m] = ... + [obj.H, obj.HI, obj.Dp, obj.Dm, obj.e_l, obj.e_r] = ... sbp.implementations.d1_upwind_6(m,obj.h); case 7 - [obj.H, obj.HI, obj.Dp, obj.Dm, obj.e_1, obj.e_m] = ... + [obj.H, obj.HI, obj.Dp, obj.Dm, obj.e_l, obj.e_r] = ... sbp.implementations.d1_upwind_7(m,obj.h); case 8 - [obj.H, obj.HI, obj.Dp, obj.Dm, obj.e_1, obj.e_m] = ... + [obj.H, obj.HI, obj.Dp, obj.Dm, obj.e_l, obj.e_r] = ... sbp.implementations.d1_upwind_8(m,obj.h); case 9 - [obj.H, obj.HI, obj.Dp, obj.Dm, obj.e_1, obj.e_m] = ... + [obj.H, obj.HI, obj.Dp, obj.Dm, obj.e_l, obj.e_r] = ... sbp.implementations.d1_upwind_9(m,obj.h); otherwise error('Invalid operator order %d.',order);