Mercurial > repos > public > sbplib
diff +sbp/+implementations/d4_compatible_6.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/d4_compatible_6.m Fri Sep 09 10:41:31 2016 +0200 +++ b/+sbp/+implementations/d4_compatible_6.m Fri Sep 09 11:03:13 2016 +0200 @@ -25,6 +25,11 @@ % Vi b?rjar med normen. Notera att alla SBP operatorer delar samma norm, % vilket ?r n?dv?ndigt f?r stabilitet + + BP = 8; + if(m<2*BP) + error(['Operator requires at least ' num2str(2*BP) ' grid points']); + end H=diag(ones(m,1),0); H_U=[0.7493827e7 / 0.25401600e8 0 0 0 0 0 0 0; 0 0.5534051e7 / 0.3628800e7 0 0 0 0 0 0; 0 0 0.104561e6 / 0.403200e6 0 0 0 0 0; 0 0 0 0.260503e6 / 0.145152e6 0 0 0 0; 0 0 0 0 0.43237e5 / 0.103680e6 0 0 0; 0 0 0 0 0 0.514081e6 / 0.403200e6 0 0; 0 0 0 0 0 0 0.3356179e7 / 0.3628800e7 0; 0 0 0 0 0 0 0 0.25631027e8 / 0.25401600e8;];