Mercurial > repos > public > sbplib
diff +sbp/D2VariableCompatibleHollow.m @ 1305:b5907140c069 feature/poroelastic
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 |
line wrap: on
line diff
--- a/+sbp/D2VariableCompatibleHollow.m Sat Jul 11 07:25:06 2020 -0700 +++ b/+sbp/D2VariableCompatibleHollow.m Mon Jul 13 14:21:41 2020 -0700 @@ -59,10 +59,10 @@ % D2 = Hinv * (-M + br*er*d1r^T - bl*el*d1l^T); % Replace d1' by e'*D1 in D2. - D2_compatible = @(b) D2(b) - obj.HI*(b(m)*e_r*d1_r' - b(m)*e_r*e_r'*D1) ... - + obj.HI*(b(1)*e_l*d1_l' - b(1)*e_l*e_l'*D1); + % D2_compatible = @(b) D2(b) - obj.HI*(b(m)*e_r*d1_r' - b(m)*e_r*e_r'*D1) ... + % + obj.HI*(b(1)*e_l*d1_l' - b(1)*e_l*e_l'*D1); - obj.D2 = D2_compatible; + obj.D2 = D2; obj.d1_l = (e_l'*D1)'; obj.d1_r = (e_r'*D1)';