Mercurial > repos > public > sbplib
diff +sbp/+implementations/d2_blocknorm_4.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/d2_blocknorm_4.m Fri Sep 09 10:41:31 2016 +0200 +++ b/+sbp/+implementations/d2_blocknorm_4.m Fri Sep 09 11:03:13 2016 +0200 @@ -1,4 +1,10 @@ function [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = d2_blocknorm_4(m,h) + + BP = 4; + if(m<2*BP) + error(['Operator requires at least ' num2str(2*BP) ' grid points']); + end + H_U=[0.751e3 / 0.3456e4 0.661e3 / 0.3456e4 -0.515e3 / 0.3456e4 0.5e1 / 0.128e3; 0.661e3 / 0.3456e4 0.1405e4 / 0.1152e4 -0.3e1 / 0.128e3 0.29e2 / 0.3456e4; -0.515e3 / 0.3456e4 -0.3e1 / 0.128e3 0.989e3 / 0.1152e4 0.149e3 / 0.3456e4; 0.5e1 / 0.128e3 0.29e2 / 0.3456e4 0.149e3 / 0.3456e4 0.3407e4 / 0.3456e4;];