Mercurial > repos > public > sbplib
diff +sbp/+implementations/d1_upwind_9.m @ 266:bfa130b7abf6 operator_remake
Added error message for too few grid points to all implementation files.
| author | Martin Almquist <martin.almquist@it.uu.se> |
|---|---|
| date | Fri, 09 Sep 2016 11:03:13 +0200 |
| parents | 6009f2712d13 |
| children | f7ac3cd6eeaa |
line wrap: on
line diff
--- a/+sbp/+implementations/d1_upwind_9.m Fri Sep 09 10:41:31 2016 +0200 +++ b/+sbp/+implementations/d1_upwind_9.m Fri Sep 09 11:03:13 2016 +0200 @@ -1,5 +1,10 @@ function [H, HI, Dp, Dm, e_1, e_m] = d1_upwind_9(m,h) + BP = 8; + if(m<2*BP) + error(['Operator requires at least ' num2str(2*BP) ' grid points']); + end + Hv=ones(m,1); Hv(1:8) = [0.1070017e7/0.3628800e7; 0.5537111e7/0.3628800e7; 0.103613e6/0.403200e6; 0.261115e6/0.145152e6; 0.298951e6/0.725760e6; 0.515677e6/0.403200e6; 0.3349879e7/0.3628800e7; 0.3662753e7/0.3628800e7;]; Hv(m-7:m)=rot90(Hv(1:8),2);