Mercurial > repos > public > sbplib

diff +sbp/+implementations/d1_upwind_3.m @ 425:e56dbd9e4196 feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge feature/beams
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 07 Feb 2017 16:09:02 +0100
parents 52b4cdf27633
children
line wrap: on
line diff
--- a/+sbp/+implementations/d1_upwind_3.m	Tue Feb 07 15:47:51 2017 +0100
+++ b/+sbp/+implementations/d1_upwind_3.m	Tue Feb 07 16:09:02 2017 +0100
@@ -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);