diff --git a/NonLinearSolver/materialLaw/mlawHyperelastic.cpp b/NonLinearSolver/materialLaw/mlawHyperelastic.cpp
index 5ea592ab218bd21ef9a1847f2f479ebb0534ecd3..94bb565dcf7eced452ed820356b44add1813b196 100644
--- a/NonLinearSolver/materialLaw/mlawHyperelastic.cpp
+++ b/NonLinearSolver/materialLaw/mlawHyperelastic.cpp
@@ -142,8 +142,8 @@ double mlawHyperViscoElastic::deformationEnergy(const IPHyperViscoElastic& q) co
     
     double Av, dAv,B_d;
     static STensor3 dB_d;
-    evaluatePhiPCorrection(trEe, devEpr, A_v, dAprdEepr, B_d, dB_vddevEe);
-    Psy = getUpdatedBulkModulus(&q)*A_v*(0.5*trEe*trEe)+getUpdatedShearModulus(&q)*B_d*STensorOperation::doubledot(devEe,devEe); //this is not correct we should do an integral
+    evaluatePhiPCorrection(trEe, devEe, Av, dAv, B_d, dB_d);
+    Psy = getUpdatedBulkModulus(&q)*Av*(0.5*trEe*trEe)+getUpdatedShearModulus(&q)*B_d*STensorOperation::doubledot(devEe,devEe); //this is not correct we should do an integral
     
 
     for (int i=0; i<_N; i++){