Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d1_upwind_3.m @ 313:52b4cdf27633 feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Cleaning order 2.
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 23 Sep 2016 21:58:52 +0200
parents f7ac3cd6eeaa
children
comparison
equal deleted inserted replaced
312:9230c056a574 313:52b4cdf27633
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];