Mercurial > repos > public > sbplib
comparison +sbp/D2BlockNorm.m @ 264:8a625c5a3633 operator_remake
Changed input parameter L (domain length) to lim (cell with domain boundaries)
author | Martin Almquist <martin.almquist@it.uu.se> |
---|---|
date | Fri, 09 Sep 2016 10:39:12 +0200 |
parents | 21a180acbd49 |
children | 4b9310edcdf8 |
comparison
equal
deleted
inserted
replaced
263:21a180acbd49 | 264:8a625c5a3633 |
---|---|
17 end | 17 end |
18 | 18 |
19 | 19 |
20 | 20 |
21 methods | 21 methods |
22 function obj = D2BlockNorm(m,L,order) | 22 function obj = D2BlockNorm(m,lim,order) |
23 | 23 |
24 x_l = lim{1}; | |
25 x_r = lim{2}; | |
26 L = x_r-x_l; | |
24 obj.h = L/(m-1); | 27 obj.h = L/(m-1); |
25 obj.x = linspace(0,L,m)'; | 28 obj.x = linspace(x_l,x_r,m)'; |
26 | 29 |
27 if order == 4 | 30 if order == 4 |
28 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_1,... | 31 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_1,... |
29 obj.e_m, obj.M, obj.Q, obj.S_1, obj.S_m] = ... | 32 obj.e_m, obj.M, obj.Q, obj.S_1, obj.S_m] = ... |
30 sbp.implementations.d2_blocknorm_4(m,obj.h); | 33 sbp.implementations.d2_blocknorm_4(m,obj.h); |