Mercurial > repos > public > sbplib

annotate +sbp/D4Variable.m @ 324:c0cbffcf6513 feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Created new operator class for lonely D4 operators. Removed some output parameters of implementations.
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 26 Sep 2016 09:51:45 +0200
parents 9230c056a574
children 31d6698c1edf
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
324
c0cbffcf6513 Created new operator class for lonely D4 operators. Removed some output parameters of implementations.
Jonatan Werpers <jonatan@werpers.com>
parents: 312
diff changeset
1 classdef D4Variable < sbp.OpSet
252
07fa0d6a05bb Renamned class files and added nonequidistant operators.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff changeset
2 properties
308
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
3 D1 % SBP operator approximating first derivative
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
4 H % Norm matrix
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
5 HI % H^-1
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
6 Q % Skew-symmetric matrix
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
7 e_l % Left boundary operator
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
8 e_r % Right boundary operator
268
4b9310edcdf8 Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents: 264
diff changeset
9 d1_l % Left boundary first derivative
4b9310edcdf8 Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents: 264
diff changeset
10 d1_r % Right boundary first derivative
4b9310edcdf8 Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents: 264
diff changeset
11 d2_l % Left boundary second derivative
4b9310edcdf8 Renamned boundary operators!
Martin Almquist <martin.almquist@it.uu.se>
parents: 264
diff changeset
12 d2_r % Right boundary second derivative
324
c0cbffcf6513 Created new operator class for lonely D4 operators. Removed some output parameters of implementations.
Jonatan Werpers <jonatan@werpers.com>
parents: 312
diff changeset
13 d3_l % Left boundary third derivative
c0cbffcf6513 Created new operator class for lonely D4 operators. Removed some output parameters of implementations.
Jonatan Werpers <jonatan@werpers.com>
parents: 312
diff changeset
14 d3_r % Right boundary third derivative
c0cbffcf6513 Created new operator class for lonely D4 operators. Removed some output parameters of implementations.
Jonatan Werpers <jonatan@werpers.com>
parents: 312
diff changeset
15 D2 % SBP operator for second derivative
c0cbffcf6513 Created new operator class for lonely D4 operators. Removed some output parameters of implementations.
Jonatan Werpers <jonatan@werpers.com>
parents: 312
diff changeset
16 M % Norm matrix, second derivative
308
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
17 D4 % SBP operator for fourth derivative
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
18 M4 % Norm matrix, fourth derivative
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
19 m % Number of grid points.
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
20 h % Step size
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
21 x % grid
263
21a180acbd49 Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents: 262
diff changeset
22 borrowing % Struct with borrowing limits for different norm matrices
252
07fa0d6a05bb Renamned class files and added nonequidistant operators.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff changeset
23 end
07fa0d6a05bb Renamned class files and added nonequidistant operators.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff changeset
24
07fa0d6a05bb Renamned class files and added nonequidistant operators.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff changeset
25 methods
324
c0cbffcf6513 Created new operator class for lonely D4 operators. Removed some output parameters of implementations.
Jonatan Werpers <jonatan@werpers.com>
parents: 312
diff changeset
26 function obj = D4Variable(m, lim, order)
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 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
28 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
29 L = x_r-x_l;
263
21a180acbd49 Renamed standard class to D2Standard etc, Rewrote class properties.
Martin Almquist <martin.almquist@it.uu.se>
parents: 262
diff changeset
30 obj.h = L/(m-1);
307
fefb2f9884f7 Merge with default.
Jonatan Werpers <jonatan@werpers.com>
parents: 268
diff changeset
31 obj.x = linspace(x_l, x_r,m)';
252
07fa0d6a05bb Renamned class files and added nonequidistant operators.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff changeset
32
07fa0d6a05bb Renamned class files and added nonequidistant operators.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff changeset
33 if order == 2
312
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
34 [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
310
ffa5d557942b Moved operator implementations and fixed some naming.
Jonatan Werpers <jonatan@werpers.com>
parents: 309
diff changeset
35 sbp.implementations.d4_variable_2(m, obj.h);
252
07fa0d6a05bb Renamned class files and added nonequidistant operators.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff changeset
36 obj.borrowing.N.S2 = 1.2500;
07fa0d6a05bb Renamned class files and added nonequidistant operators.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff changeset
37 obj.borrowing.N.S3 = 0.4000;
07fa0d6a05bb Renamned class files and added nonequidistant operators.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff changeset
38 elseif order == 4
324
c0cbffcf6513 Created new operator class for lonely D4 operators. Removed some output parameters of implementations.
Jonatan Werpers <jonatan@werpers.com>
parents: 312
diff changeset
39 [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
c0cbffcf6513 Created new operator class for lonely D4 operators. Removed some output parameters of implementations.
Jonatan Werpers <jonatan@werpers.com>
parents: 312
diff changeset
40 sbp.implementations.d4_variable_4(m, obj.h);
c0cbffcf6513 Created new operator class for lonely D4 operators. Removed some output parameters of implementations.
Jonatan Werpers <jonatan@werpers.com>
parents: 312
diff changeset
41 obj.borrowing.N.S2 = 0.5055;
c0cbffcf6513 Created new operator class for lonely D4 operators. Removed some output parameters of implementations.
Jonatan Werpers <jonatan@werpers.com>
parents: 312
diff changeset
42 obj.borrowing.N.S3 = 0.9290;
252
07fa0d6a05bb Renamned class files and added nonequidistant operators.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff changeset
43 elseif order == 6
324
c0cbffcf6513 Created new operator class for lonely D4 operators. Removed some output parameters of implementations.
Jonatan Werpers <jonatan@werpers.com>
parents: 312
diff changeset
44 [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
c0cbffcf6513 Created new operator class for lonely D4 operators. Removed some output parameters of implementations.
Jonatan Werpers <jonatan@werpers.com>
parents: 312
diff changeset
45 sbp.implementations.d4_variable_6(m, obj.h);
c0cbffcf6513 Created new operator class for lonely D4 operators. Removed some output parameters of implementations.
Jonatan Werpers <jonatan@werpers.com>
parents: 312
diff changeset
46 obj.borrowing.N.S2 = 0.3259;
c0cbffcf6513 Created new operator class for lonely D4 operators. Removed some output parameters of implementations.
Jonatan Werpers <jonatan@werpers.com>
parents: 312
diff changeset
47 obj.borrowing.N.S3 = 0.1580;
252
07fa0d6a05bb Renamned class files and added nonequidistant operators.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff changeset
48 else
07fa0d6a05bb Renamned class files and added nonequidistant operators.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff changeset
49 error('Invalid operator order.');
07fa0d6a05bb Renamned class files and added nonequidistant operators.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff changeset
50 end
07fa0d6a05bb Renamned class files and added nonequidistant operators.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff changeset
51
07fa0d6a05bb Renamned class files and added nonequidistant operators.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff changeset
52 obj.m = m;
307
fefb2f9884f7 Merge with default.
Jonatan Werpers <jonatan@werpers.com>
parents: 268
diff changeset
53
308
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
54 obj.H = H;
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
55 obj.HI = HI;
324
c0cbffcf6513 Created new operator class for lonely D4 operators. Removed some output parameters of implementations.
Jonatan Werpers <jonatan@werpers.com>
parents: 312
diff changeset
56 obj.D1 = D1;
308
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
57 obj.D2 = D2;
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
58 obj.D4 = D4;
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
59 obj.M4 = M4;
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
60 obj.e_l = e_l;
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
61 obj.e_r = e_r;
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
62 obj.d1_l = d1_l;
312
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
63 obj.d1_r = d1_r;
308
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
64 obj.d2_l = d2_l;
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
65 obj.d2_r = d2_r;
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
66 obj.d3_l = d3_l;
067fd13ba320 Improved form of D4CompatibleVariable.
Jonatan Werpers <jonatan@werpers.com>
parents: 307
diff changeset
67 obj.d3_r = d3_r;
252
07fa0d6a05bb Renamned class files and added nonequidistant operators.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff changeset
68 end
07fa0d6a05bb Renamned class files and added nonequidistant operators.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff changeset
69 end
07fa0d6a05bb Renamned class files and added nonequidistant operators.
Martin Almquist <martin.almquist@it.uu.se>
parents:
diff changeset
70 end