diff --git a/NonLinearSolver/materialLaw/mlawHyperelastic.cpp b/NonLinearSolver/materialLaw/mlawHyperelastic.cpp
index 08a3b53bd47aa0e9175f87a8f6d1f525cd5d04a9..5ca1ad71c2aa90e291f8961952e5f5cd8cb2e063 100644
--- a/NonLinearSolver/materialLaw/mlawHyperelastic.cpp
+++ b/NonLinearSolver/materialLaw/mlawHyperelastic.cpp
@@ -892,7 +892,7 @@ void mlawPowerYieldHyper::evaluatePhiPCorrection(const IPHyperViscoElastic *q1,
   
   for (int i=0; i<3; i++){
         for (int j=0; j<3; j++){
-          J_dev(i,j)=2.*getUpdatedShearModulus(q1)*getDevCorrection()*pow(exp(getThetaDevCorrection()*devE.dotprod())-1.,getPiDevCorrection())*devE(i,j)-K_dev(i,j) ;
+          J_dev(i,j)=2.*getUpdatedShearModulus(q1)*getDevCorrection()*pow(exp(getThetaDevCorrection()*Temp_Jdev)-1.,getPiDevCorrection())*devE(i,j)-K_dev(i,j) ;
         }
   }
   
@@ -1059,7 +1059,7 @@ void mlawPowerYieldHyper::evaluatePhiPCorrection(const IPHyperViscoElastic *q1,
   			for (int l=0; l<3; l++){
   				for (int m=0; m<3; m++){
   					for (int n=0; n<3; n++){
-  						H(i,j,k,l) += (1.-6.*Ge*Gamma/u)*G(i,j,n,m)*_Idev(m,n,k,l);  				
+  						H(i,j,k,l) += (1.-6.*Ge*Gamma/u)*G(i,j,m,n)*_Idev(n,m,k,l);  				
   					}
   				}
   			}
@@ -1072,7 +1072,6 @@ void mlawPowerYieldHyper::evaluatePhiPCorrection(const IPHyperViscoElastic *q1,
   	for (int j=0; j<3; j++){
   		for (int k=0; k<3; k++){
   			for (int l=0; l<3; l++){
-  				//S_dev(i,j,k,l) = (1.+6.*getUpdatedShearModulus(q1)*getDevCorrection()*pow(D_dev-1.,getPiDevCorrection())*Gamma/u)*_I4(i,j,k,l)+6.*Gamma*D_dev*getThetaDevCorrection()/u*G(i,j,j,i);
   				S_dev(i,j,k,l) = (1.+6.*getUpdatedShearModulus(q1)*getDevCorrection()*pow(D_dev-1.,getPiDevCorrection())*Gamma/u)*_I4(i,j,k,l)+6.*Gamma*D_dev*getThetaDevCorrection()/u*G(i,j,k,l);
   				U_dev(i,j,k,l) = 2.*D_dev*getThetaDevCorrection()*H(i,j,k,l)+2.*getUpdatedShearModulus(q1)*getDevCorrection()*pow(D_dev-1.,getPiDevCorrection())*(1.-6.*Ge*Gamma/u)*_Idev(i,j,k,l);
   			}
@@ -1333,7 +1332,7 @@ void mlawPowerYieldHyper::predictorCorrector_nonAssociatedFlow(const STensor3& F
        
         evaluatePhiPCorrection(q1, _b, v, trEepr, Gamma, ptildepr_Vol, ptildepr, Ke, u, devEepr, K_dev_pr, devPhipr, Ge, q1->getRefTo_ptilde_Vol(), q1->getRefTo_X_vol(), q1->getRefTo_dptilde_VoldGamma(), q1->getRefTo_dptilde_VoldCepr(), q1->getRefTo_K_dev(), q1->getRefTo_D_dev(), q1->getRefTo_dK_devdGamma(), q1->getRefTo_dK_devdCepr()); 
         
-        //Anadir: STensor3 devEepr, double u , Ge, K_dev_pr. Rehacer derivadas, la Gt no deberia estar ahi
+
         
         ptilde_Vol = q1->getRefTo_ptilde_Vol();
         Dptilde_vol = ptilde_Vol-ptildepr_Vol;
@@ -1520,7 +1519,6 @@ void mlawPowerYieldHyper::predictorCorrector_nonAssociatedFlow(const STensor3& F
           dfDCepr(i,j) =  -(a(1)/v)*(dptildeOldpr_VoldCepr(i,j)+dptilde_VoldCepr(i,j)); 
           for (int k=0; k<3; k++){
             for (int l=0; l<3; l++){
-              //dAdCepr(i,j) += (9./(A*u*u))*devPhipr(k,l)*(DdevKprDCepr(k,l,i,j)+dK_devdCepr(k,l,i,j));
               dAdCepr(i,j) += (9./(A*u*u))*devPhipr_Dev_cor(k,l)*(DdevKprDCepr(k,l,i,j)+dK_devdCepr(k,l,i,j));
               dfDCepr(i,j) += fact*devPhipr_Dev_cor(k,l)*(DdevKprDCepr(k,l,i,j)+dK_devdCepr(k,l,i,j)); 
             }