diff --git a/dG3D/src/dG3DMaterialLaw.cpp b/dG3D/src/dG3DMaterialLaw.cpp index 1bace907b7f346f7afef2298cf4b8e9751966336..19d1a3d283f99a7df6d5e0fb7849dbfa208bee3a 100644 --- a/dG3D/src/dG3DMaterialLaw.cpp +++ b/dG3D/src/dG3DMaterialLaw.cpp @@ -1723,7 +1723,7 @@ void MultipleDG3DMaterialLaw::stress(IPVariable* ipv, const IPVariable* ipvp, co IPVariable* ipv_i = ipvcur->getIPv(i); const IPVariable* ipvprev_i = ipvprev->getIPv(i); double wi = ipvcur->getWeight(i); - if (wi > 1e-10 && wi<-1e-10) + if (wi > 1e-10 || wi<-1e-10) { dG3DIPVariableBase* ipvdG3D_i = dynamic_cast<dG3DIPVariableBase*>(ipv_i); ipvdG3D_i->getRefToDeformationGradient() = F;