Mercurial > repos > public > sbplib

annotate +sbp/D1StaggeredUpwind.m @ 1344:b4e5e45bd239 feature/D2_boundary_opt

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove round off zeros from D2Nonequidistant operators
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Sat, 15 Oct 2022 15:48:20 +0200
parents a3f2c1781612
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
1259
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
1 classdef D1StaggeredUpwind < sbp.OpSet
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
2 % Compatible staggered and upwind operators by Ken Mattsson and Ossian O'reilly
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
3 properties
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
4 % x_primal: "primal" grid with m points. Equidistant. Called Plus grid in Ossian's paper.
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
5 % x_dual: "dual" grid with m+1 points. Called Minus grid in Ossian's paper.
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
6
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
7 % D1_primal takes FROM dual grid TO primal grid
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
8 % D1_dual takes FROM primal grid TO dual grid
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
9
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
10 D1_primal % SBP operator approximating first derivative
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
11 D1_dual % SBP operator approximating first derivative
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
12
1261
a3f2c1781612 Fix bug with incorrect grids returned by D1StaggeredUpwind
Martin Almquist <malmquist@stanford.edu>
parents: 1259
diff changeset
13 Dplus_primal % Upwind operator on primal grid
a3f2c1781612 Fix bug with incorrect grids returned by D1StaggeredUpwind
Martin Almquist <malmquist@stanford.edu>
parents: 1259
diff changeset
14 Dminus_primal % Upwind operator on primal grid
a3f2c1781612 Fix bug with incorrect grids returned by D1StaggeredUpwind
Martin Almquist <malmquist@stanford.edu>
parents: 1259
diff changeset
15 Dplus_dual % Upwind operator on dual grid
a3f2c1781612 Fix bug with incorrect grids returned by D1StaggeredUpwind
Martin Almquist <malmquist@stanford.edu>
parents: 1259
diff changeset
16 Dminus_dual % Upwind operator on dual grid
1259
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
17
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
18 H_primal % Norm matrix
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
19 H_dual % Norm matrix
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
20 H_primalI % H^-1
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
21 H_dualI % H^-1
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
22 e_primal_l % Left boundary operator
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
23 e_dual_l % Left boundary operator
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
24 e_primal_r % Right boundary operator
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
25 e_dual_r % Right boundary operator
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
26 m % Number of grid points.
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
27 m_primal % Number of grid points.
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
28 m_dual % Number of grid points.
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
29 h % Step size
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
30 x_primal % grid
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
31 x_dual % grid
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
32 x
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
33 borrowing % Struct with borrowing limits for different norm matrices
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
34 end
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
35
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
36 methods
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
37 function obj = D1StaggeredUpwind(m,lim,order)
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
38
1261
a3f2c1781612 Fix bug with incorrect grids returned by D1StaggeredUpwind
Martin Almquist <malmquist@stanford.edu>
parents: 1259
diff changeset
39 xl = lim{1};
a3f2c1781612 Fix bug with incorrect grids returned by D1StaggeredUpwind
Martin Almquist <malmquist@stanford.edu>
parents: 1259
diff changeset
40 xr = lim{2};
a3f2c1781612 Fix bug with incorrect grids returned by D1StaggeredUpwind
Martin Almquist <malmquist@stanford.edu>
parents: 1259
diff changeset
41 L = xr-xl;
a3f2c1781612 Fix bug with incorrect grids returned by D1StaggeredUpwind
Martin Almquist <malmquist@stanford.edu>
parents: 1259
diff changeset
42 h = L/(m-1);
1259
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
43
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
44 m_primal = m;
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
45 m_dual = m+1;
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
46
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
47 switch order
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
48 case 2
1261
a3f2c1781612 Fix bug with incorrect grids returned by D1StaggeredUpwind
Martin Almquist <malmquist@stanford.edu>
parents: 1259
diff changeset
49 [~, ~, obj.H_primal, obj.H_dual,...
1259
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
50 obj.H_primalI, obj.H_dualI,...
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
51 obj.D1_primal, obj.D1_dual, obj.Dplus_primal, obj.Dminus_primal,...
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
52 obj.Dplus_dual, obj.Dminus_dual] = sbp.implementations.d1_staggered_upwind_2(m, L);
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
53 case 4
1261
a3f2c1781612 Fix bug with incorrect grids returned by D1StaggeredUpwind
Martin Almquist <malmquist@stanford.edu>
parents: 1259
diff changeset
54 [~, ~, obj.H_primal, obj.H_dual,...
1259
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
55 obj.H_primalI, obj.H_dualI,...
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
56 obj.D1_primal, obj.D1_dual, obj.Dplus_primal, obj.Dminus_primal,...
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
57 obj.Dplus_dual, obj.Dminus_dual] = sbp.implementations.d1_staggered_upwind_4(m, L);
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
58 case 6
1261
a3f2c1781612 Fix bug with incorrect grids returned by D1StaggeredUpwind
Martin Almquist <malmquist@stanford.edu>
parents: 1259
diff changeset
59 [~, ~, obj.H_primal, obj.H_dual,...
1259
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
60 obj.H_primalI, obj.H_dualI,...
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
61 obj.D1_primal, obj.D1_dual, obj.Dplus_primal, obj.Dminus_primal,...
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
62 obj.Dplus_dual, obj.Dminus_dual] = sbp.implementations.d1_staggered_upwind_6(m, L);
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
63 case 8
1261
a3f2c1781612 Fix bug with incorrect grids returned by D1StaggeredUpwind
Martin Almquist <malmquist@stanford.edu>
parents: 1259
diff changeset
64 [~, ~, obj.H_primal, obj.H_dual,...
1259
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
65 obj.H_primalI, obj.H_dualI,...
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
66 obj.D1_primal, obj.D1_dual, obj.Dplus_primal, obj.Dminus_primal,...
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
67 obj.Dplus_dual, obj.Dminus_dual] = sbp.implementations.d1_staggered_upwind_8(m, L);
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
68 otherwise
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
69 error('Invalid operator order %d.',order);
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
70 end
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
71
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
72 obj.m = m;
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
73 obj.m_primal = m_primal;
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
74 obj.m_dual = m_dual;
1261
a3f2c1781612 Fix bug with incorrect grids returned by D1StaggeredUpwind
Martin Almquist <malmquist@stanford.edu>
parents: 1259
diff changeset
75 obj.h = h;
a3f2c1781612 Fix bug with incorrect grids returned by D1StaggeredUpwind
Martin Almquist <malmquist@stanford.edu>
parents: 1259
diff changeset
76
a3f2c1781612 Fix bug with incorrect grids returned by D1StaggeredUpwind
Martin Almquist <malmquist@stanford.edu>
parents: 1259
diff changeset
77 obj.x_primal = linspace(xl, xr, m)';
a3f2c1781612 Fix bug with incorrect grids returned by D1StaggeredUpwind
Martin Almquist <malmquist@stanford.edu>
parents: 1259
diff changeset
78 obj.x_dual = [xl, linspace(xl+h/2, xr-h/2, m-1), xr]';
1259
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
79
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
80 obj.e_primal_l = sparse(m_primal,1);
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
81 obj.e_primal_r = sparse(m_primal,1);
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
82 obj.e_primal_l(1) = 1;
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
83 obj.e_primal_r(m_primal) = 1;
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
84
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
85 obj.e_dual_l = sparse(m_dual,1);
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
86 obj.e_dual_r = sparse(m_dual,1);
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
87 obj.e_dual_l(1) = 1;
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
88 obj.e_dual_r(m_dual) = 1;
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
89
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
90 obj.borrowing = [];
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
91 obj.x = [];
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
92
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
93 end
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
94
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
95 function str = string(obj)
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
96 str = [class(obj) '_' num2str(obj.order)];
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
97 end
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
98 end
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
99 end