Mercurial > repos > public > sbplib

annotate +sbp/upwind2.m @ 132:6ec2248b83c4

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Refactored +anim.animate to be more clear.
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 05 Feb 2016 10:49:16 +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