diff --git a/dG3D/src/dG3DMaterialLaw.cpp b/dG3D/src/dG3DMaterialLaw.cpp
index 994fb533b43c88e047c2b556ad5ad47daa693a27..acf40d0067af7eca03ae7195ff5186e0aae23d52 100644
--- a/dG3D/src/dG3DMaterialLaw.cpp
+++ b/dG3D/src/dG3DMaterialLaw.cpp
@@ -2965,7 +2965,10 @@ void HyperViscoElasticDG3DMaterialLaw::stress(IPVariable* ipv, const IPVariable*
 
 HyperViscoElastoPlasticPowerYieldDG3DMaterialLaw::HyperViscoElastoPlasticPowerYieldDG3DMaterialLaw(const int num, const double rho, const double E,const double nu, const double tol ,
                         const bool matrixbyPerturbation, const double pert):
-                        dG3DMaterialLaw(num,rho), _viscoLaw(num,E,nu,rho,tol,matrixbyPerturbation,pert){};
+                        dG3DMaterialLaw(num,rho), _viscoLaw(num,E,nu,rho,tol,matrixbyPerturbation,pert)
+{
+  fillElasticStiffness(E, nu, elasticStiffness);
+};
 HyperViscoElastoPlasticPowerYieldDG3DMaterialLaw::HyperViscoElastoPlasticPowerYieldDG3DMaterialLaw(const HyperViscoElastoPlasticPowerYieldDG3DMaterialLaw& src): dG3DMaterialLaw(src),
                       _viscoLaw(src._viscoLaw){};