diff --git a/dG3D/src/dG3DMaterialLaw.cpp b/dG3D/src/dG3DMaterialLaw.cpp
index fcb797284d9ca244997433be8646a363d43bf259..c568c2a25494c8d60d2b1e753ab125493a9384e0 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){};