diff --git a/NonLinearSolver/materialLaw/mlawNonLinearTVP.cpp b/NonLinearSolver/materialLaw/mlawNonLinearTVP.cpp
index 88b4e0ff3f66b89321e8f3dc1d36c25a4d1073c6..47241a4449d7010345a9dd9b1b4129a3e3d634dc 100644
--- a/NonLinearSolver/materialLaw/mlawNonLinearTVP.cpp
+++ b/NonLinearSolver/materialLaw/mlawNonLinearTVP.cpp
@@ -1394,7 +1394,7 @@ void mlawNonLinearTVP::plasticCorrector_TVP(const STensor3& F0, const STensor3&
 
   corKir_e = q1->_kirchhoff;
   const STensor3& KS = q1->_kirchhoff;
-  const STensor3& DcorKirDT_TVEcorrection = q1->_DcorKirDT; // Get the corrected value directly from TVE
+  // const STensor3& DcorKirDT_TVEcorrection = q1->_DcorKirDT; // Get the corrected value directly from TVE
   getModifiedMandel(Ce, q0, q1); // update Mandel
 
   const STensor3& MS = q1->_ModMandel;
@@ -2288,14 +2288,15 @@ void mlawNonLinearTVP::predictorCorrector_TVP_nonAssociatedFlow(const STensor3&
         dDgammaDGamma = kk*(A+Gamma*dAdGamma);  // mistake in the paper (VD 2016)
 
         this->getYieldCoefficientDerivatives(q1,q1->_nup,Da);
-             
+        
+        /*     
         // DEBUG FLE
         double a0 = a(0);
         double a1 = a(1);
         double a2 = a(2);
         double da0 = Da(0);
         double da1 = Da(1);
-        double da2 = Da(2);
+        double da2 = Da(2);*/
         
         dfdDgamma = Da(2)*pow(PhiEq,_n) - Da(1)*ptilde -Da(0);
         dfdDgamma += a(2)*_n*pow(PhiEq,(_n-1.))*He - a(1)*Hp;
@@ -3311,7 +3312,7 @@ void mlawNonLinearTVP::predictorCorrector_TVP_nonAssociatedFlow(const STensor3&
 
 
   // mechSourceTVP derivatives
-
+/*
     // Numerical Derivatives 
     double T_plus(0.), Wm_TVE_plus(0.), Wm_TVP_plus(0.), dWm_TVEdT_plus(0.), dWm_TVPdT_plus(0.), DKDTsum_plus(0.), DGDTsum_plus(0.);
     static STensor3 F_plus, Kepr_plus, Mepr_plus, Ke_plus, Me_plus, N_plus, P_plus, DphiPDF_plus, dWm_TVEdF_plus, dWm_TVPdF_plus;
@@ -3342,7 +3343,7 @@ void mlawNonLinearTVP::predictorCorrector_TVP_nonAssociatedFlow(const STensor3&
         dWm_TVEdF_plus(i,j) = (Wm_TVE_plus - Wm_TVE)/_perturbationfactor;
         dWm_TVPdF_plus(i,j) = (Wm_TVP_plus - Wm_TVP)/_perturbationfactor;
       }
-    }
+    }*/
 
   // conversion tensor for mechSourceTVE Terms
   static STensor43 DCeDFe, DEeDFe;