Mercurial > repos > public > sbplib

comparison +sbp/D4Lonely.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 +sbp/D4Variable.m@9230c056a574
children 72468bc9b63f
comparison
equal deleted inserted replaced
323:7579c2abbf9f 324:c0cbffcf6513
1 classdef D4Lonely < sbp.OpSet
2 properties
3 H % Norm matrix
4 HI % H^-1
5 e_l, e_r % Left and right boundary operator
6 d1_l, d1_r % Left and right boundary first derivative
7 d2_l, d2_r % Left and right boundary second derivative
8 d3_l, d3_r % Left and right boundary third derivative
9 D4 % SBP operator for fourth derivative
10 M4 % Norm matrix, fourth derivative
11 m % Number of grid points.
12 h % Step size
13 x % grid
14 borrowing % Struct with borrowing limits for different norm matrices
15 end
16
17 methods
18 function obj = D4Lonely(m, lim, order, opt)
19 default_arg('opt', '')
20
21 x_l = lim{1};
22 x_r = lim{2};
23 L = x_r-x_l;
24 obj.h = L/(m-1);
25 obj.x = linspace(x_l, x_r,m)';
26
27 if order == 2
28 [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
29 sbp.implementations.d4_variable_2(m, obj.h);
30 obj.borrowing.N.S2 = 1.2500;
31 obj.borrowing.N.S3 = 0.4000;
32
33 elseif order == 4
34 switch opt
35 case 'min_boundary_points'
36 [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
37 sbp.implementations.d4_variable_4_min_boundary_points(m, obj.h);
38 % obj.borrowing.N.S2 = 0.5055;
39 % obj.borrowing.N.S3 = 0.9290;
40 otherwise
41 [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
42 sbp.implementations.d4_variable_4(m, obj.h);
43 obj.borrowing.N.S2 = 0.5055;
44 obj.borrowing.N.S3 = 0.9290;
45 end
46
47 elseif order == 6
48 switch opt
49 case '2'
50 [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
51 sbp.implementations.d4_variable_6_2(m, obj.h);
52 % obj.borrowing.N.S2 = 0.3259;
53 % obj.borrowing.N.S3 = 0.1580;
54 case '3'
55 [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
56 sbp.implementations.d4_variable_6_3(m, obj.h);
57 % obj.borrowing.N.S2 = 0.3259;
58 % obj.borrowing.N.S3 = 0.1580;
59 case 'min_boundary_points'
60 [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
61 sbp.implementations.d4_variable_6_min_boundary_points(m, obj.h);
62 % obj.borrowing.N.S2 = 0.3259;
63 % obj.borrowing.N.S3 = 0.1580;
64 otherwise
65 [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
66 sbp.implementations.d4_variable_6(m, obj.h);
67 obj.borrowing.N.S2 = 0.3259;
68 obj.borrowing.N.S3 = 0.1580;
69 end
70
71 elseif order == 8
72 switch opt
73 case 'min_boundary_points'
74 [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
75 sbp.implementations.d4_variable_8_min_boundary_points(m, obj.h);
76 % obj.borrowing.N.S2 = 0.3259;
77 % obj.borrowing.N.S3 = 0.1580;
78 otherwise
79 [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ...
80 sbp.implementations.d4_variable_8_higher_boundary_order(m, obj.h);
81 % obj.borrowing.N.S2 = 0.3259;
82 % obj.borrowing.N.S3 = 0.1580;
83 end
84 else
85 error('Invalid operator order.');
86 end
87
88 obj.m = m;
89
90 obj.H = H;
91 obj.HI = HI;
92 obj.D4 = D4;
93 obj.M4 = M4;
94 obj.e_l = e_l;
95 obj.e_r = e_r;
96 obj.d1_l = d1_l;
97 obj.d1_r = d1_r;
98 obj.d2_l = d2_l;
99 obj.d2_r = d2_r;
100 obj.d3_l = d3_l;
101 obj.d3_r = d3_r;
102 end
103 end
104 end