Mercurial > repos > public > sbplib

annotate +sbp/upwind2.m @ 163:51aaf67a7df5 feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed naming and added functions for converting vectors to matrices for structured grid.
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 22 Feb 2016 13:55:19 +0100
parents c56437d097de
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
119
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
1 function [H, HI, Dp, Dm, e_1, e_m] = upwind2(m,h)
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
2 Hv=ones(m,1);
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
3 Hv(1:2)=[0.1e1/0.4e1; 0.5e1/0.4e1;];
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
4 Hv(m-1:m)=rot90(Hv(1:2),2);
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
5 Hv = Hv*h;
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
6 H = spdiag(Hv,0);
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
7 HI = spdiag(1./Hv,0);
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
8
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
9 q_diags = [0 1 2];
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
10 q_stencil = [-3/2 +2 -1/2];
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
11 Qp = stripeMatrix(q_stencil, q_diags,m);
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
12
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
13 Q_U = [
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
14 -0.1e1/0.4e1 0.5e1/0.4e1;
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
15 -0.1e1/0.4e1 -0.5e1/0.4e1;
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
16 ];
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
17
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
18 Qp(1:2,1:2)=Q_U;
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
19 Qp(m-1:m,m-1:m)=rot90(Q_U,2)'; %%% This is different from standard SBP
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
20
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
21 Qm=-Qp';
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
22
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
23 e_1=zeros(m,1);e_1(1)=1;
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
24 e_m=zeros(m,1);e_m(m)=1;
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
25
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
26 Dp=HI*(Qp-1/2*e_1*e_1'+1/2*e_m*e_m') ;
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
27
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
28 Dm=HI*(Qm-1/2*e_1*e_1'+1/2*e_m*e_m') ;
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
29 end
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
30
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
31
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
32
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
33 % KOLLA NAMN PÅ ALLA FUNKTIONER
c56437d097de I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
34 % TESTA ALLA FUNKTIONER