Mercurial > repos > public > sbplib
diff +sbp/D2Standard.m @ 263:21a180acbd49 operator_remake
Renamed standard class to D2Standard etc, Rewrote class properties.
| author | Martin Almquist <martin.almquist@it.uu.se> |
|---|---|
| date | Thu, 08 Sep 2016 17:50:30 +0200 |
| parents | |
| children | 8a625c5a3633 |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/+sbp/D2Standard.m Thu Sep 08 17:50:30 2016 +0200 @@ -0,0 +1,71 @@ +classdef D2Standard < sbp.OpSet + properties + D1 % SBP operator approximating first derivative + H % Norm matrix + HI % H^-1 + Q % Skew-symmetric matrix + e_1 % Left boundary operator + e_m % Right boundary operator + D2 % SBP operator for second derivative + M % Norm matrix, second derivative + S_1 % Left boundary first derivative + S_m % Right boundary first derivative + m % Number of grid points. + h % Step size + x % grid + borrowing % Struct with borrowing limits for different norm matrices + + end + + methods + function obj = D2Standard(m,L,order) + + obj.h = L/(m-1); + obj.x = linspace(0,L,m)'; + + if order == 2 + [obj.H, obj.HI, obj.D1, obj.D2, obj.e_1,... + obj.e_m, obj.M, obj.Q, obj.S_1, obj.S_m] =... + sbp.implementations.ordinary2(m,obj.h); + obj.borrowing.M.S = 0.4000; + elseif order == 4 + [obj.H, obj.HI, obj.D1, obj.D2, obj.e_1,... + obj.e_m, obj.M, obj.Q, obj.S_1, obj.S_m] =... + sbp.implementations.ordinary4(m,obj.h); + obj.borrowing.M.S = 0.2508; + elseif order == 6 + [obj.H, obj.HI, obj.D1, obj.D2, obj.e_1,... + obj.e_m, obj.M, obj.Q, obj.S_1, obj.S_m] =... + sbp.implementations.ordinary6(m,obj.h); + obj.borrowing.M.S = 0.1878; + elseif order == 8 + [obj.H, obj.HI, obj.D1, obj.D2, obj.e_1,... + obj.e_m, obj.M, obj.Q, obj.S_1, obj.S_m] =... + sbp.implementations.ordinary8(m,obj.h); + obj.borrowing.M.S = 0.0015; + elseif order == 10 + [obj.H, obj.HI, obj.D1, obj.D2, obj.e_1,... + obj.e_m, obj.M, obj.Q, obj.S_1, obj.S_m] =... + sbp.implementations.ordinary10(m,obj.h); + obj.borrowing.M.S = 0.0351; + elseif order == 12 + [obj.H, obj.HI, obj.D1, obj.D2, obj.e_1,... + obj.e_m, obj.M, obj.Q, obj.S_1, obj.S_m] =... + sbp.implementations.ordinary12(m,obj.h); + obj.borrowing.M.S = []; % Need to compute! + else + error('Invalid operator order %d.',order); + end + + obj.m = m; + + end + end + + +end + + + + +