Mercurial > repos > public > sbplib

comparison +sbp/D2VariablePeriodic.m @ 1337:bf2554f1825d feature/D2_boundary_opt

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add periodic D1 and D2 operators for orders 8,10,12
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 13 May 2022 13:28:10 +0200
parents 5ccf6aaf6d6b
children
comparison
equal deleted inserted replaced
1336:0666629aa183 1337:bf2554f1825d
25 obj.h = L/m; 25 obj.h = L/m;
26 x = linspace(x_l,x_r,m+1)'; 26 x = linspace(x_l,x_r,m+1)';
27 obj.x = x(1:end-1); 27 obj.x = x(1:end-1);
28 28
29 switch order 29 switch order
30
31 case 12
32 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,...
33 obj.e_r, obj.d1_l, obj.d1_r] = ...
34 sbp.implementations.d2_variable_periodic_12(m,obj.h);
35 obj.borrowing.M.d1 = obj.H(1, 1) / obj.h; % First element in H/h is borrowing also for M
36 obj.borrowing.R.delta_D = inf;
37
38 case 10
39 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,...
40 obj.e_r, obj.d1_l, obj.d1_r] = ...
41 sbp.implementations.d2_variable_periodic_10(m,obj.h);
42 obj.borrowing.M.d1 = obj.H(1, 1) / obj.h; % First element in H/h is borrowing also for M
43 obj.borrowing.R.delta_D = inf;
44
45 case 8
46 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,...
47 obj.e_r, obj.d1_l, obj.d1_r] = ...
48 sbp.implementations.d2_variable_periodic_8(m,obj.h);
49 obj.borrowing.M.d1 = obj.H(1, 1) / obj.h; % First element in H/h is borrowing also for M
50 obj.borrowing.R.delta_D = inf;
30 51
31 case 6 52 case 6
32 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... 53 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,...
33 obj.e_r, obj.d1_l, obj.d1_r] = ... 54 obj.e_r, obj.d1_l, obj.d1_r] = ...
34 sbp.implementations.d2_variable_periodic_6(m,obj.h); 55 sbp.implementations.d2_variable_periodic_6(m,obj.h);