diff --git a/nonlinelast.cc b/nonlinelast.cc index c29a256046b24704a8fd880a5900ac85bb62c0b5..795b27da31efc20b174e086c944fbff27ab5e0f8 100644 --- a/nonlinelast.cc +++ b/nonlinelast.cc @@ -1,6 +1,5 @@ #include <config.h> -#include <dune/fufem/utilities/adolcnamespaceinjections.hh> #include <dune/common/parametertreeparser.hh> #include <dune/common/bitsetvector.hh> @@ -29,12 +28,11 @@ #include <dune/solvers/transferoperators/mandelobsrestrictor.hh> #include <dune/solvers/transferoperators/truncatedcompressedmgtransfer.hh> #include <dune/solvers/iterationsteps/trustregiongsstep.hh> +#include <dune/solvers/solvers/trustregionsolver.hh> #include <dune/elasticity/common/nonlinearelasticityproblem.hh> -#include <dune/elasticity/common/trustregionsolver.hh> #define LAURSEN #include <dune/elasticity/materials/neohookeanmaterial.hh> -#include <dune/elasticity/materials/neohookematerial.hh> // The grid dimension const int dim = 3; @@ -143,11 +141,10 @@ int main (int argc, char *argv[]) try // typedef P1NodalBasis<GridType::LeafGridView> P1Basis; P1Basis p1Basis(grid->leafGridView()); - typedef NeoHookeMaterial<P1Basis> MaterialType; + typedef NeoHookeanMaterial<P1Basis> MaterialType; MaterialType material(p1Basis, parameterSet.get<double>("E"), - parameterSet.get<double>("nu"), - parameterSet.get<bool>("vectorMode")); + parameterSet.get<double>("nu")); // /////////////////////////////////////////////////// @@ -287,10 +284,9 @@ int main (int argc, char *argv[]) try typedef P2NodalBasis<GridType::LeafGridView> P2Basis; P2Basis p2Basis(grid->leafGridView()); - typedef NeoHookeMaterial<P2Basis> MaterialTypeP2; + typedef NeoHookeanMaterial<P2Basis> MaterialTypeP2; MaterialTypeP2 p2Material(p2Basis, parameterSet.get<field_type>("E"), - parameterSet.get<field_type>("nu"), - parameterSet.get<bool>("vectorMode")); + parameterSet.get<field_type>("nu")); // P2 Forces VectorType p2ExtForces(p2Basis.size());