[MINOR PATCH] Aliter to the previous

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);
     }