Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d1_upwind_2.m @ 261:6009f2712d13 operator_remake

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Moved and renamned all implementations.
author Martin Almquist <martin.almquist@it.uu.se>
date Thu, 08 Sep 2016 15:35:45 +0200
parents
children bfa130b7abf6
comparison
equal deleted inserted replaced
260:b4116ce49ac4 261:6009f2712d13
1 function [H, HI, Dp, Dm, e_1, e_m] = d1_upwind_2(m,h)
2 Hv=ones(m,1);
3 Hv(1:2)=[0.1e1/0.4e1; 0.5e1/0.4e1;];
4 Hv(m-1:m)=rot90(Hv(1:2),2);
5 Hv = Hv*h;
6 H = spdiag(Hv,0);
7 HI = spdiag(1./Hv,0);
8
9 q_diags = [0 1 2];
10 q_stencil = [-3/2 +2 -1/2];
11 Qp = stripeMatrix(q_stencil, q_diags,m);
12
13 Q_U = [
14 -0.1e1/0.4e1 0.5e1/0.4e1;
15 -0.1e1/0.4e1 -0.5e1/0.4e1;
16 ];
17
18 Qp(1:2,1:2)=Q_U;
19 Qp(m-1:m,m-1:m)=rot90(Q_U,2)'; %%% This is different from standard SBP
20
21 Qm=-Qp';
22
23 e_1=zeros(m,1);e_1(1)=1;
24 e_m=zeros(m,1);e_m(m)=1;
25
26 Dp=HI*(Qp-1/2*e_1*e_1'+1/2*e_m*e_m') ;
27
28 Dm=HI*(Qm-1/2*e_1*e_1'+1/2*e_m*e_m') ;
29 end
30