bug

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++){