changeset 677:eeaf9a00e304 feature/poroelastic

Correct factor of 2 in Dirichlet penalty parameter.
author Martin Almquist <malmquist@stanford.edu>
date Fri, 26 Jan 2018 15:17:49 -0800
parents 9926efb39330
children 06676c40e77f
files +scheme/elasticShearVariable.m
diffstat 1 files changed, 1 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/elasticShearVariable.m	Thu Jan 18 14:36:59 2018 -0800
+++ b/+scheme/elasticShearVariable.m	Fri Jan 26 15:17:49 2018 -0800
@@ -237,7 +237,7 @@
                     for i = 1:obj.dim
                         H_gamma_i = obj.H_boundary{i};
                         sigma_ij = tuning*delta(i,j)*2/(gamma*h(j)) +...
-                                   tuning*delta_b(i,j)*(2/(H11*h(j)) + 2/(H11*h(j)*phi));
+                                   tuning*delta_b(i,j)*(2/(H11*h(j)) + 1/(H11*h(j)*phi));
 
                         ci = E{i}*inv(RHO)*nj*Hi*...
                                 ( (e{j}*H_gamma*e{j}'*A*e{j}*d{j}')'*E{i}' + ...