Mercurial > repos > public > sbplib
diff +sbp/D4Lonely.m @ 324:c0cbffcf6513 feature/beams
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 |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/+sbp/D4Lonely.m Mon Sep 26 09:51:45 2016 +0200 @@ -0,0 +1,104 @@ +classdef D4Lonely < sbp.OpSet + properties + H % Norm matrix + HI % H^-1 + e_l, e_r % Left and right boundary operator + d1_l, d1_r % Left and right boundary first derivative + d2_l, d2_r % Left and right boundary second derivative + d3_l, d3_r % Left and right boundary third derivative + D4 % SBP operator for fourth derivative + M4 % Norm matrix, fourth derivative + m % Number of grid points. + h % Step size + x % grid + borrowing % Struct with borrowing limits for different norm matrices + end + + methods + function obj = D4Lonely(m, lim, order, opt) + default_arg('opt', '') + + x_l = lim{1}; + x_r = lim{2}; + L = x_r-x_l; + obj.h = L/(m-1); + obj.x = linspace(x_l, x_r,m)'; + + if order == 2 + [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... + sbp.implementations.d4_variable_2(m, obj.h); + obj.borrowing.N.S2 = 1.2500; + obj.borrowing.N.S3 = 0.4000; + + elseif order == 4 + switch opt + case 'min_boundary_points' + [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... + sbp.implementations.d4_variable_4_min_boundary_points(m, obj.h); + % obj.borrowing.N.S2 = 0.5055; + % obj.borrowing.N.S3 = 0.9290; + otherwise + [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... + sbp.implementations.d4_variable_4(m, obj.h); + obj.borrowing.N.S2 = 0.5055; + obj.borrowing.N.S3 = 0.9290; + end + + elseif order == 6 + switch opt + case '2' + [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... + sbp.implementations.d4_variable_6_2(m, obj.h); + % obj.borrowing.N.S2 = 0.3259; + % obj.borrowing.N.S3 = 0.1580; + case '3' + [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... + sbp.implementations.d4_variable_6_3(m, obj.h); + % obj.borrowing.N.S2 = 0.3259; + % obj.borrowing.N.S3 = 0.1580; + case 'min_boundary_points' + [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... + sbp.implementations.d4_variable_6_min_boundary_points(m, obj.h); + % obj.borrowing.N.S2 = 0.3259; + % obj.borrowing.N.S3 = 0.1580; + otherwise + [H, HI, D1, D2, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... + sbp.implementations.d4_variable_6(m, obj.h); + obj.borrowing.N.S2 = 0.3259; + obj.borrowing.N.S3 = 0.1580; + end + + elseif order == 8 + switch opt + case 'min_boundary_points' + [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... + sbp.implementations.d4_variable_8_min_boundary_points(m, obj.h); + % obj.borrowing.N.S2 = 0.3259; + % obj.borrowing.N.S3 = 0.1580; + otherwise + [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = ... + sbp.implementations.d4_variable_8_higher_boundary_order(m, obj.h); + % obj.borrowing.N.S2 = 0.3259; + % obj.borrowing.N.S3 = 0.1580; + end + else + error('Invalid operator order.'); + end + + obj.m = m; + + obj.H = H; + obj.HI = HI; + obj.D4 = D4; + obj.M4 = M4; + obj.e_l = e_l; + obj.e_r = e_r; + obj.d1_l = d1_l; + obj.d1_r = d1_r; + obj.d2_l = d2_l; + obj.d2_r = d2_r; + obj.d3_l = d3_l; + obj.d3_r = d3_r; + end + end +end \ No newline at end of file