Mercurial > repos > public > sbplib

diff +scheme/Elastic2dVariable.m @ 967:368a2773f78b feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix bug in Elastic2dVariable that resulted in using tuning^2
author Martin Almquist <malmquist@stanford.edu>
date Sun, 23 Dec 2018 14:06:26 +0100
parents 4a85cff53107
children a4ad90b37998
line wrap: on
line diff
--- a/+scheme/Elastic2dVariable.m	Sat Dec 22 18:08:51 2018 +0100
+++ b/+scheme/Elastic2dVariable.m	Sun Dec 23 14:06:26 2018 +0100
@@ -483,8 +483,8 @@
 
             [alpha_ii_u, alpha_ij_u] = computeAlpha(theta_R_u, theta_H_u, theta_M_u, lambda_u, mu_u);
             [alpha_ii_v, alpha_ij_v] = computeAlpha(theta_R_v, theta_H_v, theta_M_v, lambda_v, mu_v);
-            sigma_ii = tuning*(alpha_ii_u + alpha_ii_v);
-            sigma_ij = tuning*(alpha_ij_u + alpha_ij_v);
+            sigma_ii = alpha_ii_u + alpha_ii_v;
+            sigma_ij = alpha_ij_u + alpha_ij_v;
 
             d = @kroneckerDelta;  % Kronecker delta
             db = @(i,j) 1-d(i,j); % Logical not of Kronecker delta