Mercurial > repos > public > sbplib

comparison +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
comparison
equal deleted inserted replaced
426:29944ea7674b 427:a613960a157b
1 function [H, HI, Dp, Dm, e_1, e_m] = d1_upwind_3(m,h) 1 function [H, HI, Dp, Dm, e_1, e_m] = d1_upwind_3(m,h)
2 2
3 if(m<6) 3 if(m<6)
4 error('Operator requires at least 6 grid points'); 4 error('Operator requires at least 6 grid points');
5 end 5 end
6 6
7 Hv = ones(m,1); 7 Hv = ones(m,1);
8 Hv(1:3) = [3/8; 7/6; 23/24]; 8 Hv(1:3) = [3/8; 7/6; 23/24];
9 Hv(m-2:m) = rot90(Hv(1:3),2); 9 Hv(m-2:m) = rot90(Hv(1:3), 2);
10 Hv = Hv*h; 10 Hv = Hv*h;
11 H = spdiag(Hv,0); 11 H = spdiag(Hv,0);
12 HI = spdiag(1./Hv,0); 12 HI = spdiag(1./Hv,0);
13 13
14 q_diags = [-1, 0, 1, 2]; 14 q_diags = [-1, 0, 1, 2];