Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d1_upwind_4.m @ 425:e56dbd9e4196 feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge feature/beams
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 07 Feb 2017 16:09:02 +0100
parents 88584b0cfba1
children
comparison
equal deleted inserted replaced
423:a2cb0d4f4a02 425:e56dbd9e4196
21 0.11e2/0.144e3 -0.13e2/0.32e2 -0.29e2/0.48e2 0.389e3/0.288e3; 21 0.11e2/0.144e3 -0.13e2/0.32e2 -0.29e2/0.48e2 0.389e3/0.288e3;
22 0.1e1/0.32e2 -0.11e2/0.144e3 -0.65e2/0.288e3 -0.13e2/0.16e2; 22 0.1e1/0.32e2 -0.11e2/0.144e3 -0.65e2/0.288e3 -0.13e2/0.16e2;
23 ]; 23 ];
24 24
25 Qp(1:4,1:4)=Q_U; 25 Qp(1:4,1:4)=Q_U;
26 Qp(m-3:m,m-3:m)=rot90(Q_U,2)'; %%% This is different from standard SBP 26 Qp(m-3:m,m-3:m)=rot90(Q_U, 2)'; %%% This is different from standard SBP
27 27
28 Qm=-Qp'; 28 Qm=-Qp';
29 29
30 e_1=sparse(m,1);e_1(1)=1; 30 e_1=sparse(m,1);e_1(1)=1;
31 e_m=sparse(m,1);e_m(m)=1; 31 e_m=sparse(m,1);e_m(m)=1;