Mercurial > repos > public > sbplib
annotate +sbp/D2Standard.m @ 1289:2fd2e2337b77 feature/boundary_optimized_grids
Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
| author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
|---|---|
| date | Wed, 01 Jul 2020 15:15:30 +0200 |
| parents | bc78157c89cb |
| children | eca4ca84cf0a |
| rev | line source |
|---|---|
|
263
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
1 classdef D2Standard < sbp.OpSet |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
2 properties |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
3 D1 % SBP operator approximating first derivative |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
4 H % Norm matrix |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
5 HI % H^-1 |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
6 Q % Skew-symmetric matrix |
|
268
4b9310edcdf8
Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents:
265
diff
changeset
|
7 e_l % Left boundary operator |
|
4b9310edcdf8
Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents:
265
diff
changeset
|
8 e_r % Right boundary operator |
|
263
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
9 D2 % SBP operator for second derivative |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
10 M % Norm matrix, second derivative |
|
268
4b9310edcdf8
Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents:
265
diff
changeset
|
11 d1_l % Left boundary first derivative |
|
4b9310edcdf8
Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents:
265
diff
changeset
|
12 d1_r % Right boundary first derivative |
|
263
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
13 m % Number of grid points. |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
14 h % Step size |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
15 x % grid |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
16 borrowing % Struct with borrowing limits for different norm matrices |
| 396 | 17 |
|
263
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
18 end |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
19 |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
20 methods |
|
264
8a625c5a3633
Changed input parameter L (domain length) to lim (cell with domain boundaries)
Martin Almquist <martin.almquist@it.uu.se>
parents:
263
diff
changeset
|
21 function obj = D2Standard(m,lim,order) |
|
263
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
22 |
|
264
8a625c5a3633
Changed input parameter L (domain length) to lim (cell with domain boundaries)
Martin Almquist <martin.almquist@it.uu.se>
parents:
263
diff
changeset
|
23 x_l = lim{1}; |
|
8a625c5a3633
Changed input parameter L (domain length) to lim (cell with domain boundaries)
Martin Almquist <martin.almquist@it.uu.se>
parents:
263
diff
changeset
|
24 x_r = lim{2}; |
|
8a625c5a3633
Changed input parameter L (domain length) to lim (cell with domain boundaries)
Martin Almquist <martin.almquist@it.uu.se>
parents:
263
diff
changeset
|
25 L = x_r-x_l; |
|
263
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
26 obj.h = L/(m-1); |
|
264
8a625c5a3633
Changed input parameter L (domain length) to lim (cell with domain boundaries)
Martin Almquist <martin.almquist@it.uu.se>
parents:
263
diff
changeset
|
27 obj.x = linspace(x_l,x_r,m)'; |
|
263
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
28 |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
29 if order == 2 |
|
268
4b9310edcdf8
Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents:
265
diff
changeset
|
30 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... |
|
4b9310edcdf8
Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents:
265
diff
changeset
|
31 obj.e_r, obj.M, obj.Q, obj.d1_l, obj.d1_r] =... |
| 269 | 32 sbp.implementations.d2_2(m,obj.h); |
|
263
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
33 obj.borrowing.M.S = 0.4000; |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
34 elseif order == 4 |
|
268
4b9310edcdf8
Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents:
265
diff
changeset
|
35 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... |
|
4b9310edcdf8
Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents:
265
diff
changeset
|
36 obj.e_r, obj.M, obj.Q, obj.d1_l, obj.d1_r] =... |
| 269 | 37 sbp.implementations.d2_4(m,obj.h); |
|
263
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
38 obj.borrowing.M.S = 0.2508; |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
39 elseif order == 6 |
|
268
4b9310edcdf8
Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents:
265
diff
changeset
|
40 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... |
|
4b9310edcdf8
Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents:
265
diff
changeset
|
41 obj.e_r, obj.M, obj.Q, obj.d1_l, obj.d1_r] =... |
| 269 | 42 sbp.implementations.d2_6(m,obj.h); |
|
263
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
43 obj.borrowing.M.S = 0.1878; |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
44 elseif order == 8 |
|
268
4b9310edcdf8
Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents:
265
diff
changeset
|
45 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... |
|
4b9310edcdf8
Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents:
265
diff
changeset
|
46 obj.e_r, obj.M, obj.Q, obj.d1_l, obj.d1_r] =... |
| 269 | 47 sbp.implementations.d2_8(m,obj.h); |
|
263
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
48 obj.borrowing.M.S = 0.0015; |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
49 elseif order == 10 |
|
268
4b9310edcdf8
Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents:
265
diff
changeset
|
50 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... |
|
4b9310edcdf8
Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents:
265
diff
changeset
|
51 obj.e_r, obj.M, obj.Q, obj.d1_l, obj.d1_r] =... |
| 269 | 52 sbp.implementations.d2_10(m,obj.h); |
|
263
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
53 obj.borrowing.M.S = 0.0351; |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
54 elseif order == 12 |
|
268
4b9310edcdf8
Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents:
265
diff
changeset
|
55 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... |
|
4b9310edcdf8
Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents:
265
diff
changeset
|
56 obj.e_r, obj.M, obj.Q, obj.d1_l, obj.d1_r] =... |
| 269 | 57 sbp.implementations.d2_12(m,obj.h); |
|
263
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
58 obj.borrowing.M.S = []; % Need to compute! |
|
265
3e0432b3cb68
Added warning for 12th order 2nd der, D2=D1*D1
Martin Almquist <martin.almquist@it.uu.se>
parents:
264
diff
changeset
|
59 warning('D2=D1*D1 (wide stencil).'); |
|
3e0432b3cb68
Added warning for 12th order 2nd der, D2=D1*D1
Martin Almquist <martin.almquist@it.uu.se>
parents:
264
diff
changeset
|
60 warning('The borrowing constant has not been computed.') |
|
263
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
61 else |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
62 error('Invalid operator order %d.',order); |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
63 end |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
64 |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
65 obj.m = m; |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
66 end |
|
402
bc78157c89cb
Add string methods to some sbp classes.
Jonatan Werpers <jonatan@werpers.com>
parents:
397
diff
changeset
|
67 |
|
bc78157c89cb
Add string methods to some sbp classes.
Jonatan Werpers <jonatan@werpers.com>
parents:
397
diff
changeset
|
68 function str = string(obj) |
|
bc78157c89cb
Add string methods to some sbp classes.
Jonatan Werpers <jonatan@werpers.com>
parents:
397
diff
changeset
|
69 str = [class(obj) '_' num2str(obj.order)]; |
|
bc78157c89cb
Add string methods to some sbp classes.
Jonatan Werpers <jonatan@werpers.com>
parents:
397
diff
changeset
|
70 end |
|
bc78157c89cb
Add string methods to some sbp classes.
Jonatan Werpers <jonatan@werpers.com>
parents:
397
diff
changeset
|
71 |
|
263
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
72 end |
|
21a180acbd49
Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff
changeset
|
73 end |