Mercurial > repos > public > sbplib
changeset 966:4a85cff53107 feature/poroelastic
Fix bug in Elastic2dVariable that caused crash for m_x neq m_y
| author | Martin Almquist <malmquist@stanford.edu> |
|---|---|
| date | Sat, 22 Dec 2018 18:08:51 +0100 |
| parents | db3411264b96 |
| children | 368a2773f78b |
| files | +scheme/Elastic2dVariable.m |
| diffstat | 1 files changed, 4 insertions(+), 4 deletions(-) [+] |
line wrap: on
line diff
diff -r db3411264b96 -r 4a85cff53107 +scheme/Elastic2dVariable.m --- a/+scheme/Elastic2dVariable.m Sat Dec 22 16:41:17 2018 +0100 +++ b/+scheme/Elastic2dVariable.m Sat Dec 22 18:08:51 2018 +0100 @@ -264,13 +264,13 @@ tau_r{j}{i} = sparse(n_r, dim*m_tot); for k = 1:dim T_l{j}{i,k} = ... - -d(i,j)*LAMBDA_l*(d(i,k)*d1_l{k}' + db(i,k)*e_l{j}'*D1{k})... - -d(j,k)*MU_l*(d(i,j)*d1_l{i}' + db(i,j)*e_l{j}'*D1{i})... + -d(i,j)*LAMBDA_l*(d(i,k)*d1_l{j}' + db(i,k)*e_l{j}'*D1{k})... + -d(j,k)*MU_l*(d(i,j)*d1_l{j}' + db(i,j)*e_l{j}'*D1{i})... -d(i,k)*MU_l*d1_l{j}'; T_r{j}{i,k} = ... - d(i,j)*LAMBDA_r*(d(i,k)*d1_r{k}' + db(i,k)*e_r{j}'*D1{k})... - +d(j,k)*MU_r*(d(i,j)*d1_r{i}' + db(i,j)*e_r{j}'*D1{i})... + d(i,j)*LAMBDA_r*(d(i,k)*d1_r{j}' + db(i,k)*e_r{j}'*D1{k})... + +d(j,k)*MU_r*(d(i,j)*d1_r{j}' + db(i,j)*e_r{j}'*D1{i})... +d(i,k)*MU_r*d1_r{j}'; tau_l{j}{i} = tau_l{j}{i} + T_l{j}{i,k}*E{k}';