Mercurial > repos > public > sbplib

comparison +sbp/D2VariablePeriodic.m @ 686:5ccf6aaf6d6b feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add D2VariablePeriodic orders 4 and 6.
author Martin Almquist <malmquist@stanford.edu>
date Thu, 08 Feb 2018 16:44:46 -0800
parents 7368affc8f78
children bf2554f1825d
comparison
equal deleted inserted replaced
685:b035902869a8 686:5ccf6aaf6d6b
27 obj.x = x(1:end-1); 27 obj.x = x(1:end-1);
28 28
29 switch order 29 switch order
30 30
31 case 6 31 case 6
32 error('Not impl') 32 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,...
33 33 obj.e_r, obj.d1_l, obj.d1_r] = ...
34 [obj.H, obj.HI, obj.D1, obj.D2, ... 34 sbp.implementations.d2_variable_periodic_6(m,obj.h);
35 ~, obj.e_l, obj.e_r, ~, ~, ~, ~, ~,...
36 obj.d1_l, obj.d1_r] = ...
37 sbp.implementations.d4_variable_periodic_6(m, obj.h);
38 obj.borrowing.M.d1 = 0.1878; 35 obj.borrowing.M.d1 = 0.1878;
39 obj.borrowing.R.delta_D = 0.3696; 36 obj.borrowing.R.delta_D = 0.3696;
40 % Borrowing e^T*D1 - d1 from R 37 % Borrowing e^T*D1 - d1 from R
41 38
42 case 4 39 case 4
43 error('Not impl')
44
45 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... 40 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,...
46 obj.e_r, obj.d1_l, obj.d1_r] = ... 41 obj.e_r, obj.d1_l, obj.d1_r] = ...
47 sbp.implementations.d2_variable_periodic_4(m,obj.h); 42 sbp.implementations.d2_variable_periodic_4(m,obj.h);
48 obj.borrowing.M.d1 = 0.2505765857; 43 obj.borrowing.M.d1 = 0.2505765857;
49 44
61 56
62 otherwise 57 otherwise
63 error('Invalid operator order %d.',order); 58 error('Invalid operator order %d.',order);
64 end 59 end
65 obj.borrowing.H11 = obj.H(1,1)/obj.h; % First element in H/h, 60 obj.borrowing.H11 = obj.H(1,1)/obj.h; % First element in H/h,
61
66 obj.m = m; 62 obj.m = m;
67 obj.M = []; 63 obj.M = [];
68 end 64 end
69 function str = string(obj) 65 function str = string(obj)
70 str = [class(obj) '_' num2str(obj.order)]; 66 str = [class(obj) '_' num2str(obj.order)];