Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d1_upwind_3.m @ 832:5573913a0949 feature/burgers1d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merged with default, and updated +scheme/Burgers1D accordingly
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Tue, 11 Sep 2018 15:58:35 +0200
parents 52b4cdf27633
children
comparison
equal deleted inserted replaced
831:d0934d1143b7 832:5573913a0949
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];