23 std::unique_ptr<SimplifiedElasticityModel> simplified_elasticity =
24 std::make_unique<RigidElasticityModel>();
25 if (
auto const simplified_elasticity_switch =
29 DBUG(
"Using simplified_elasticity for the Richards flow equation");
30 if (*simplified_elasticity_switch ==
"uniaxial")
32 DBUG(
"assuming local uniaxial deformation only.");
33 simplified_elasticity = std::make_unique<UniaxialElasticityModel>();
35 else if (*simplified_elasticity_switch ==
"hydrostatic")
37 DBUG(
"assuming constant hydrostatic stress locally.");
38 simplified_elasticity =
39 std::make_unique<HydrostaticElasticityModel>();
41 else if (*simplified_elasticity_switch ==
"user_defined")
43 DBUG(
"using user defined elasticity model.");
44 simplified_elasticity =
45 std::make_unique<UserDefinedElasticityModel>();
47 else if (*simplified_elasticity_switch ==
"rigid")
49 DBUG(
"using user defined elasticity model.");
50 simplified_elasticity = std::make_unique<RigidElasticityModel>();
53 return simplified_elasticity;