diff --git a/dG3D/src/dG3DMaterialLaw.cpp b/dG3D/src/dG3DMaterialLaw.cpp
index 8cf971c47c2765eed55b0c24a7f6620195b28012..eba380b8d03f0d00a5c045f3ba58e838b8cc53bc 100644
--- a/dG3D/src/dG3DMaterialLaw.cpp
+++ b/dG3D/src/dG3DMaterialLaw.cpp
@@ -5751,12 +5751,12 @@ void StochTMDMNDG3DMaterialLaw::stress(IPVariable*ipv, const IPVariable*ipvprev,
     static STensor3 RT; STensorOperation::zero(RT);
     STensorOperation::transposeSTensor3(R,RT);
 
-    // STensorOperation::multSTensor3InPlace2nd(R,U);
-    // STensorOperation::multSTensor3InPlace1st(U,RT);
+    STensorOperation::multSTensor3InPlace2nd(R,U);
+    STensorOperation::multSTensor3InPlace1st(U,RT);
     // STensorOperation::multSTensor3(R_def,U,F);
-    // F = U; // R_def in the rotated frame = I;
+    F = U; // R_def in the rotated frame = I;
 
-    STensorOperation::multSTensor3InPlace1st(F,RT);
+    // STensorOperation::multSTensor3InPlace1st(F,RT);
     STensorOperation::multSTensor3SVector3InPlace2nd(R,H);
   }
 
@@ -6283,9 +6283,9 @@ void StochTMDMNDG3DMaterialLaw::stress(IPVariable*ipv, const IPVariable*ipvprev,
         L = L_rot; dqdF = dqdF_rot; dwdf = dwdf_rot; dmechSourcedf = dmechSourcedf_rot; 
       }
     
-      // STensorOperation::multSTensor3(F_unrot,S,P);
+      STensorOperation::multSTensor3(F_unrot,S,P);
 
-      STensorOperation::multSTensor3(Prot,R,P);
+      // STensorOperation::multSTensor3(Prot,R,P);
       STensorOperation::multSTensor3SVector3InPlace2nd(RT,Q);
 
     }