diff --git a/dG3D/src/dG3DMaterialLaw.cpp b/dG3D/src/dG3DMaterialLaw.cpp index eba380b8d03f0d00a5c045f3ba58e838b8cc53bc..cc0f448383c0b17e3d2899dd7a52746a9d9af9b7 100644 --- a/dG3D/src/dG3DMaterialLaw.cpp +++ b/dG3D/src/dG3DMaterialLaw.cpp @@ -6211,14 +6211,14 @@ void StochTMDMNDG3DMaterialLaw::stress(IPVariable*ipv, const IPVariable*ipvprev, for(int s=0; s<3; s++) dUrot_dF(m,n,k,l) += I4_symm(m,n,p,q)*dUrot_dU(p,q,r,s)*dUdF(r,s,k,l); } - + /* for(int m=0; m<3; m++) for(int n=0; n<3; n++) for(int k=0; k<3; k++) for(int l=0; l<3; l++) for(int p=0; p<3; p++) - dFrot_dF(m,n,k,l) += R_def(m,p)*dUrot_dF(p,n,k,l); - // dFrot_dF = dUrot_dF; // R_def_rot = I; + dFrot_dF(m,n,k,l) += R_def(m,p)*dUrot_dF(p,n,k,l);*/ + dFrot_dF = dUrot_dF; // R_def_rot = I; for(int m=0; m<3; m++) for(int n=0; n<3; n++)