Mercurial > repos > public > sbplib

comparison +sbp/D4Compatible.m @ 264:8a625c5a3633 operator_remake

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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
25 end 25 end
26 26
27 27
28 28
29 methods 29 methods
30 function obj = D4Compatible(m,L,order) 30 function obj = D4Compatible(m,lim,order)
31 31
32
33 x_l = lim{1};
34 x_r = lim{2};
35 L = x_r-x_l;
32 obj.h = L/(m-1); 36 obj.h = L/(m-1);
33 obj.x = linspace(0,L,m)'; 37 obj.x = linspace(x_l,x_r,m)';
34 38
35 if order == 2 39 if order == 2
36 [obj.H, obj.HI, obj.D1, obj.D4, obj.e_1, obj.e_m, obj.M4,... 40 [obj.H, obj.HI, obj.D1, obj.D4, obj.e_1, obj.e_m, obj.M4,...
37 obj.Q, obj.S2_1, obj.S2_m, obj.S3_1, obj.S3_m,... 41 obj.Q, obj.S2_1, obj.S2_m, obj.S3_1, obj.S3_m,...
38 obj.S_1, obj.S_m] =... 42 obj.S_1, obj.S_m] =...