Mercurial > repos > public > sbplib

comparison +sbp/D2Variable.m @ 668:43ea848b6aa1 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add borrowing constants to D2variable
author Martin Almquist <malmquist@stanford.edu>
date Fri, 22 Dec 2017 14:14:24 +0100
parents 359861563866
children 52a9dedb9171
comparison
equal deleted inserted replaced
667:ed853945ee99 668:43ea848b6aa1
29 case 4 29 case 4
30 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... 30 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,...
31 obj.e_r, obj.d1_l, obj.d1_r] = ... 31 obj.e_r, obj.d1_l, obj.d1_r] = ...
32 sbp.implementations.d2_variable_4(m,obj.h); 32 sbp.implementations.d2_variable_4(m,obj.h);
33 obj.borrowing.M.d1 = 0.2505765857; 33 obj.borrowing.M.d1 = 0.2505765857;
34
35 obj.borrowing.R.delta_D = 0.577587500088313;
36 % Borrowing e^T*D1 - d1 from R
34 case 2 37 case 2
35 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... 38 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,...
36 obj.e_r, obj.d1_l, obj.d1_r] = ... 39 obj.e_r, obj.d1_l, obj.d1_r] = ...
37 sbp.implementations.d2_variable_2(m,obj.h); 40 sbp.implementations.d2_variable_2(m,obj.h);
38 obj.borrowing.M.d1 = 0.3636363636; 41 obj.borrowing.M.d1 = 0.3636363636;
39 % Borrowing const taken from Virta 2014 42 % Borrowing const taken from Virta 2014
43
44 obj.borrowing.R.delta_D = 1.000000538455350;
45 % Borrowing e^T*D1 - d1 from R
40 46
41 otherwise 47 otherwise
42 error('Invalid operator order %d.',order); 48 error('Invalid operator order %d.',order);
43 end 49 end
44 50 obj.borrowing.H11 = obj.H(1,1)/obj.h; % First element in H/h,
45 obj.m = m; 51 obj.m = m;
46 obj.M = []; 52 obj.M = [];
47 end 53 end
48 function str = string(obj) 54 function str = string(obj)
49 str = [class(obj) '_' num2str(obj.order)]; 55 str = [class(obj) '_' num2str(obj.order)];