Mercurial > repos > public > sbplib

comparison +sbp/D2VariableCompatibleHollow.m @ 1305:b5907140c069 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Clean up hollow D2 generation by using precomputed coefficients stored in mat files
author Martin Almquist <malmquist@stanford.edu>
date Mon, 13 Jul 2020 14:21:41 -0700
parents 49e3870335ef
children
comparison
equal deleted inserted replaced
1304:a38e80fdbf60 1305:b5907140c069
57 e_r = obj.e_r; 57 e_r = obj.e_r;
58 e_l = obj.e_l; 58 e_l = obj.e_l;
59 59
60 % D2 = Hinv * (-M + br*er*d1r^T - bl*el*d1l^T); 60 % D2 = Hinv * (-M + br*er*d1r^T - bl*el*d1l^T);
61 % Replace d1' by e'*D1 in D2. 61 % Replace d1' by e'*D1 in D2.
62 D2_compatible = @(b) D2(b) - obj.HI*(b(m)*e_r*d1_r' - b(m)*e_r*e_r'*D1) ... 62 % D2_compatible = @(b) D2(b) - obj.HI*(b(m)*e_r*d1_r' - b(m)*e_r*e_r'*D1) ...
63 + obj.HI*(b(1)*e_l*d1_l' - b(1)*e_l*e_l'*D1); 63 % + obj.HI*(b(1)*e_l*d1_l' - b(1)*e_l*e_l'*D1);
64 64
65 obj.D2 = D2_compatible; 65 obj.D2 = D2;
66 obj.d1_l = (e_l'*D1)'; 66 obj.d1_l = (e_l'*D1)';
67 obj.d1_r = (e_r'*D1)'; 67 obj.d1_r = (e_r'*D1)';
68 68
69 end 69 end
70 function str = string(obj) 70 function str = string(obj)