Mercurial > repos > public > sbplib

comparison +sbp/D2Standard.m @ 265:3e0432b3cb68 operator_remake

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added warning for 12th order 2nd der, D2=D1*D1
author Martin Almquist <martin.almquist@it.uu.se>
date Fri, 09 Sep 2016 10:41:31 +0200
parents 8a625c5a3633
children 4b9310edcdf8
comparison
equal deleted inserted replaced
264:8a625c5a3633 265:3e0432b3cb68
54 elseif order == 12 54 elseif order == 12
55 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_1,... 55 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_1,...
56 obj.e_m, obj.M, obj.Q, obj.S_1, obj.S_m] =... 56 obj.e_m, obj.M, obj.Q, obj.S_1, obj.S_m] =...
57 sbp.implementations.ordinary12(m,obj.h); 57 sbp.implementations.ordinary12(m,obj.h);
58 obj.borrowing.M.S = []; % Need to compute! 58 obj.borrowing.M.S = []; % Need to compute!
59 warning('D2=D1*D1 (wide stencil).');
60 warning('The borrowing constant has not been computed.')
59 else 61 else
60 error('Invalid operator order %d.',order); 62 error('Invalid operator order %d.',order);
61 end 63 end
62 64
63 obj.m = m; 65 obj.m = m;