Mercurial > repos > public > sbplib
diff +sbp/+implementations/d1_upwind_3.m @ 832:5573913a0949 feature/burgers1d
Merged with default, and updated +scheme/Burgers1D accordingly
author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
---|---|
date | Tue, 11 Sep 2018 15:58:35 +0200 |
parents | 52b4cdf27633 |
children |
line wrap: on
line diff
--- a/+sbp/+implementations/d1_upwind_3.m Tue Sep 11 13:24:08 2018 +0200 +++ b/+sbp/+implementations/d1_upwind_3.m Tue Sep 11 15:58:35 2018 +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);