Mercurial > repos > public > sbplib
comparison +sbp/D2BlockNorm.m @ 252:07fa0d6a05bb operator_remake
Renamned class files and added nonequidistant operators.
| author | Martin Almquist <martin.almquist@it.uu.se> |
|---|---|
| date | Wed, 07 Sep 2016 13:40:41 +0200 |
| parents | |
| children | 6009f2712d13 |
comparison
equal
deleted
inserted
replaced
| 251:6a5e94bb5e13 | 252:07fa0d6a05bb |
|---|---|
| 1 classdef D2BlockNorm < sbp.OpSet | |
| 2 properties | |
| 3 norms % Struct containing norm matrices such as H,Q, M | |
| 4 boundary % Struct contanging vectors for boundry point approximations | |
| 5 derivatives % Struct containging differentiation operators | |
| 6 borrowing % Struct with borrowing limits for different norm matrices | |
| 7 m % Number of grid points. | |
| 8 h % Step size | |
| 9 end | |
| 10 | |
| 11 | |
| 12 | |
| 13 methods | |
| 14 function obj = D2BlockNorm(m,h,order) | |
| 15 | |
| 16 if order == 4 | |
| 17 [H, HI, D1, D2, e_1, e_m, M,Q S_1, S_m] = sbp.blocknorm4(m,h); | |
| 18 elseif order == 6 | |
| 19 [H, HI, D1, D2, e_1, e_m, M,Q S_1, S_m] = sbp.blocknorm6(m,h); | |
| 20 elseif order == 8 | |
| 21 [H, HI, D1, D2, e_1, e_m, M,Q S_1, S_m] = sbp.blocknorm8(m,h); | |
| 22 elseif order == 10 | |
| 23 [H, HI, D1, D2, e_1, e_m, M,Q S_1, S_m] = sbp.blocknorm10(m,h); | |
| 24 else | |
| 25 error('Invalid operator order %d.',order); | |
| 26 end | |
| 27 | |
| 28 obj.h = h; | |
| 29 obj.m = m; | |
| 30 | |
| 31 obj.norms.H = H; | |
| 32 obj.norms.HI = HI; | |
| 33 obj.norms.Q = Q; | |
| 34 obj.norms.M = M; | |
| 35 | |
| 36 obj.boundary.e_1 = e_1; | |
| 37 obj.boundary.S_1 = S_1; | |
| 38 | |
| 39 obj.boundary.e_m = e_m; | |
| 40 obj.boundary.S_m = S_m; | |
| 41 | |
| 42 obj.derivatives.D1 = D1; | |
| 43 obj.derivatives.D2 = D2; | |
| 44 end | |
| 45 end | |
| 46 | |
| 47 methods (Static) | |
| 48 function lambda = smallestGrid(obj) | |
| 49 error('Not implmented') | |
| 50 end | |
| 51 end | |
| 52 | |
| 53 | |
| 54 | |
| 55 end |