change prestrain_material_geometry in compute_mu_Gamma to parametrized_Laminate

inputs/computeMuGamma.parset

@@ -29,7 +29,7 @@
 
 nElements = 32 32
 
-gamma=0.5
+gamma=100
 
 
 #############################################
@@ -40,15 +40,16 @@ write_materialFunctions = true   # VTK mu-functions , lambda-functions
 
 beta=2.0
 
-mu1=10.0
-
+mu1=1.0
 
+lambda1=0.0
 #lambda1=10.0
 
 
 ####  material_implementation("analytical_Example") ? 
 
-material_prestrain_imp= "analytical_Example"
+material_prestrain_imp= "parametrized_Laminate"
+#material_prestrain_imp= "analytical_Example"
 #material_prestrain_imp ="isotropic_bilayer"
 
 
@@ -63,9 +64,9 @@ write_prestrainFunctions = true  # VTK norm of B ,
 
 
 rho1=1.0
-alpha=9.764705882352942
+alpha=2.0
 
-theta=0.15000000000000002
+theta=0.541859296482412
 
 #theta = 0.3    # volume fraction   #default = 1.0/4.0
 #theta = 0.25   # volume fraction
@@ -111,12 +112,3 @@ Solvertype = 1
 
 
 
-)*mu2)
-
-
-
-
-
-
-
-

src/Compute_MuGamma.cc

@@ -744,7 +744,7 @@ int main(int argc, char *argv[])
     using Domain = GridView::Codim<0>::Geometry::GlobalCoordinate;
     
     // ----------- INPUT PARAMETERS -----------------------------
-    std::string imp = parameterSet.get<std::string>("material_prestrain_imp", "analytical_Example");
+    std::string imp = parameterSet.get<std::string>("material_prestrain_imp", "parametrized_Laminate");
     log << "material_prestrain used: "<< imp  << std::endl;
     double gamma   = parameterSet.get<double>("gamma", 1.0);   
     double theta   = parameterSet.get<double>("theta", 1.0/4.0);