Mercurial > repos > public > sbplib
annotate +sbp/D1StaggeredUpwind.m @ 757:179f234f6cbf feature/d1_staggered
Merge with poroelastic to get blockEvalOn.
author | Martin Almquist <malmquist@stanford.edu> |
---|---|
date | Sun, 17 Jun 2018 14:44:05 -0700 |
parents | 62147d8ce062 |
children | a3f2c1781612 |
rev | line source |
---|---|
643
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
1 classdef D1StaggeredUpwind < sbp.OpSet |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
2 % Compatible staggered and upwind operators by Ken Mattsson and Ossian O'reilly |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
3 properties |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
4 % x_primal: "primal" grid with m points. Equidistant. Called Plus grid in Ossian's paper. |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
5 % x_dual: "dual" grid with m+1 points. Called Minus grid in Ossian's paper. |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
6 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
7 % D1_primal takes FROM dual grid TO primal grid |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
8 % D1_dual takes FROM primal grid TO dual grid |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
9 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
10 D1_primal % SBP operator approximating first derivative |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
11 D1_dual % SBP operator approximating first derivative |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
12 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
13 Dplus_primal % Upwind operator on primal grid |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
14 Dminus_primal % Upwind operator on primal grid |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
15 Dplus_dual % Upwind operator on dual grid |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
16 Dminus_dual % Upwind operator on dual grid |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
17 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
18 H_primal % Norm matrix |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
19 H_dual % Norm matrix |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
20 H_primalI % H^-1 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
21 H_dualI % H^-1 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
22 e_primal_l % Left boundary operator |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
23 e_dual_l % Left boundary operator |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
24 e_primal_r % Right boundary operator |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
25 e_dual_r % Right boundary operator |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
26 m % Number of grid points. |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
27 m_primal % Number of grid points. |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
28 m_dual % Number of grid points. |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
29 h % Step size |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
30 x_primal % grid |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
31 x_dual % grid |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
32 x |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
33 borrowing % Struct with borrowing limits for different norm matrices |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
34 end |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
35 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
36 methods |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
37 function obj = D1StaggeredUpwind(m,lim,order) |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
38 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
39 x_l = lim{1}; |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
40 x_r = lim{2}; |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
41 L = x_r-x_l; |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
42 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
43 m_primal = m; |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
44 m_dual = m+1; |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
45 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
46 switch order |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
47 case 2 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
48 [obj.x_primal, obj.x_dual, obj.H_primal, obj.H_dual,... |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
49 obj.H_primalI, obj.H_dualI,... |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
50 obj.D1_primal, obj.D1_dual, obj.Dplus_primal, obj.Dminus_primal,... |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
51 obj.Dplus_dual, obj.Dminus_dual] = sbp.implementations.d1_staggered_upwind_2(m, L); |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
52 case 4 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
53 [obj.x_primal, obj.x_dual, obj.H_primal, obj.H_dual,... |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
54 obj.H_primalI, obj.H_dualI,... |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
55 obj.D1_primal, obj.D1_dual, obj.Dplus_primal, obj.Dminus_primal,... |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
56 obj.Dplus_dual, obj.Dminus_dual] = sbp.implementations.d1_staggered_upwind_4(m, L); |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
57 case 6 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
58 [obj.x_primal, obj.x_dual, obj.H_primal, obj.H_dual,... |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
59 obj.H_primalI, obj.H_dualI,... |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
60 obj.D1_primal, obj.D1_dual, obj.Dplus_primal, obj.Dminus_primal,... |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
61 obj.Dplus_dual, obj.Dminus_dual] = sbp.implementations.d1_staggered_upwind_6(m, L); |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
62 case 8 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
63 [obj.x_primal, obj.x_dual, obj.H_primal, obj.H_dual,... |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
64 obj.H_primalI, obj.H_dualI,... |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
65 obj.D1_primal, obj.D1_dual, obj.Dplus_primal, obj.Dminus_primal,... |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
66 obj.Dplus_dual, obj.Dminus_dual] = sbp.implementations.d1_staggered_upwind_8(m, L); |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
67 otherwise |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
68 error('Invalid operator order %d.',order); |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
69 end |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
70 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
71 obj.m = m; |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
72 obj.m_primal = m_primal; |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
73 obj.m_dual = m_dual; |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
74 obj.h = L/(m-1); |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
75 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
76 obj.e_primal_l = sparse(m_primal,1); |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
77 obj.e_primal_r = sparse(m_primal,1); |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
78 obj.e_primal_l(1) = 1; |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
79 obj.e_primal_r(m_primal) = 1; |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
80 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
81 obj.e_dual_l = sparse(m_dual,1); |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
82 obj.e_dual_r = sparse(m_dual,1); |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
83 obj.e_dual_l(1) = 1; |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
84 obj.e_dual_r(m_dual) = 1; |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
85 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
86 obj.borrowing = []; |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
87 obj.x = []; |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
88 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
89 end |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
90 |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
91 function str = string(obj) |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
92 str = [class(obj) '_' num2str(obj.order)]; |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
93 end |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
94 end |
62147d8ce062
Add staggered-upwind operators by Ken and Ossian.
Martin Almquist <malmquist@stanford.edu>
parents:
diff
changeset
|
95 end |