Mercurial > repos > public > sbplib
diff +sbp/+implementations/d1_upwind_3.m @ 313:52b4cdf27633 feature/beams
Cleaning order 2.
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Fri, 23 Sep 2016 21:58:52 +0200 |
parents | f7ac3cd6eeaa |
children |
line wrap: on
line diff
--- a/+sbp/+implementations/d1_upwind_3.m Fri Sep 23 19:14:04 2016 +0200 +++ b/+sbp/+implementations/d1_upwind_3.m Fri Sep 23 21:58:52 2016 +0200 @@ -1,12 +1,12 @@ function [H, HI, Dp, Dm, e_1, e_m] = d1_upwind_3(m,h) - + if(m<6) error('Operator requires at least 6 grid points'); end Hv = ones(m,1); Hv(1:3) = [3/8; 7/6; 23/24]; - Hv(m-2:m) = rot90(Hv(1:3),2); + Hv(m-2:m) = rot90(Hv(1:3), 2); Hv = Hv*h; H = spdiag(Hv,0); HI = spdiag(1./Hv,0);