diff --git a/NonLinearSolver/materialLaw/mlawNonLinearTVENonLinearTVP.cpp b/NonLinearSolver/materialLaw/mlawNonLinearTVENonLinearTVP.cpp
index 45e52b6c6c2ef30009e2576d741ff085897154d6..a4f98f54443c93fbd61dc5356ae8e3baaf8d86ec 100644
--- a/NonLinearSolver/materialLaw/mlawNonLinearTVENonLinearTVP.cpp
+++ b/NonLinearSolver/materialLaw/mlawNonLinearTVENonLinearTVP.cpp
@@ -583,8 +583,8 @@ void mlawNonLinearTVENonLinearTVP::predictorCorrector_TVP_nonAssociatedFlow_nonL
                         DdevPhidGamma_RHS(i,j) += ( 0.5*2*_b*Ke*1/q1->_Jacobian*(ptilde+Gamma*dPhiPdGamma)*_I(i,j) - 3*pow(kk,1)*Hb/Cxdev*devPhi(i,j) ); // pow(kk,2) DEBUG
                         for (int k=0; k<3; k++)
                             for (int l=0; l<3; l++){
-                                DdevPhidDgamma_RHS(i,j) += 0.5* 2.*_b/3.*Gamma*(-Ge_TrEeTensor(i,j,k,l)*1/q1->_Jacobian + Dho(i,j,k,l))*_I(k,l)*Hp*1/q1->_Jacobian;
-                                DdevPhidGamma_RHS(i,j) += 0.5* ( -3.*Ge_Tensor(i,j,k,l)*1/q1->_Jacobian*devPhi(k,l) + 2.*_b/3.*Gamma*(-Ge_TrEeTensor(i,j,k,l)*1/q1->_Jacobian + Dho(i,j,k,l))*_I(k,l)*dPhiPdGamma );
+                                DdevPhidDgamma_RHS(i,j) += 0.5* 2.*_b/3.*Gamma*(-Ge_TrEeTensor(i,j,k,l)*1/q1->_Jacobian + Dho(i,j,k,l)*1/q1->_Jacobian)*_I(k,l)*Hp;
+                                DdevPhidGamma_RHS(i,j) += 0.5* ( -3.*Ge_Tensor(i,j,k,l)*1/q1->_Jacobian*devPhi(k,l) + 2.*_b/3.*Gamma*(-Ge_TrEeTensor(i,j,k,l)*1/q1->_Jacobian + Dho(i,j,k,l)*1/q1->_Jacobian)*_I(k,l)*dPhiPdGamma );
                                 DdevPhidGamma_RHS(i,j) += 0.5*(-Ge_TrEeTensor(i,j,k,l)*1/q1->_Jacobian + Dho(i,j,k,l)*1/q1->_Jacobian)*N(k,l); // DEBUG
                             }
                     }