diff --git a/dG3D/src/dG3DMaterialLaw.cpp b/dG3D/src/dG3DMaterialLaw.cpp
index bca3df91422aa1cdb7b5fb0634917d6c85a3e6a5..8cf971c47c2765eed55b0c24a7f6620195b28012 100644
--- a/dG3D/src/dG3DMaterialLaw.cpp
+++ b/dG3D/src/dG3DMaterialLaw.cpp
@@ -5751,10 +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::multSTensor3(R_def,U,F);
+    // STensorOperation::multSTensor3InPlace2nd(R,U);
+    // STensorOperation::multSTensor3InPlace1st(U,RT);
+    // STensorOperation::multSTensor3(R_def,U,F);
     // F = U; // R_def in the rotated frame = I;
+
+    STensorOperation::multSTensor3InPlace1st(F,RT);
     STensorOperation::multSTensor3SVector3InPlace2nd(R,H);
   }
 
@@ -6281,7 +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::multSTensor3SVector3InPlace2nd(RT,Q);
 
     }