Mercurial > repos > public > sbplib
comparison +sbp/+implementations/d1_upwind_4.m @ 832:5573913a0949 feature/burgers1d
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 | 88584b0cfba1 |
children |
comparison
equal
deleted
inserted
replaced
831:d0934d1143b7 | 832:5573913a0949 |
---|---|
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; |