Mercurial > repos > public > sbplib

annotate +sbp/D2Standard.m @ 1037:2d7ba44340d0 feature/burgers1d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Pass scheme specific parameters as cell array. This will enabale constructDiffOps to be more general. In addition, allow for schemes returning function handles as diffOps, which is currently how non-linear schemes such as Burgers1d are implemented.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 18 Jan 2019 09:02:02 +0100
parents bc78157c89cb
children eca4ca84cf0a
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
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
3fdfad20037e Remove whitespace.
Jonatan Werpers <jonatan@werpers.com>
parents: 269
diff changeset
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
479194747f9c Bug fix
Martin Almquist <martin.almquist@it.uu.se>
parents: 268
diff changeset
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
479194747f9c Bug fix
Martin Almquist <martin.almquist@it.uu.se>
parents: 268
diff changeset
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
479194747f9c Bug fix
Martin Almquist <martin.almquist@it.uu.se>
parents: 268
diff changeset
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
479194747f9c Bug fix
Martin Almquist <martin.almquist@it.uu.se>
parents: 268
diff changeset
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
479194747f9c Bug fix
Martin Almquist <martin.almquist@it.uu.se>
parents: 268
diff changeset
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
479194747f9c Bug fix
Martin Almquist <martin.almquist@it.uu.se>
parents: 268
diff changeset
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