Mercurial > repos > public > sbplib

diff +sbp/+implementations/d1_upwind_3.m @ 820:501750fbbfdb

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with feature/grids
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 07 Sep 2018 14:40:58 +0200
parents 52b4cdf27633
children
line wrap: on
line diff
--- a/+sbp/+implementations/d1_upwind_3.m	Fri Sep 07 14:39:38 2018 +0200
+++ b/+sbp/+implementations/d1_upwind_3.m	Fri Sep 07 14:40:58 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);