Mercurial > repos > public > sbplib

diff +sbp/+implementations/d1_upwind_3.m @ 427:a613960a157b feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
merged with feature/beams
author Ylva Rydin <ylva.rydin@telia.com>
date Thu, 26 Jan 2017 15:59:25 +0100
parents 52b4cdf27633
children
line wrap: on
line diff
--- a/+sbp/+implementations/d1_upwind_3.m	Wed Jan 25 17:14:17 2017 +0100
+++ b/+sbp/+implementations/d1_upwind_3.m	Thu Jan 26 15:59:25 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);