Skip to content
Snippets Groups Projects
Commit cbb51725 authored by Elias Pipping's avatar Elias Pipping Committed by Elias Pipping
Browse files

Typedef

parent e4d2b8c3
Branches releases/2.7
No related tags found
No related merge requests found
......@@ -143,25 +143,22 @@ void assemble_nonlinearity(
Dune::shared_ptr<Dune::GlobalNonlinearity<dim> const> &myGlobalNonlinearity,
Dune::shared_ptr<Dune::BlockVector<Dune::FieldVector<double, 1>>>
nodalIntegrals) {
typedef Dune::BlockVector<Dune::FieldVector<double, 1>> SingletonVectorType;
// {{{ Assemble terms for the nonlinearity
auto mu =
Dune::make_shared<Dune::BlockVector<Dune::FieldVector<double, 1>>>(size);
auto mu = Dune::make_shared<SingletonVectorType>(size);
*mu = parset.get<double>("boundary.friction.mu");
auto normalStress =
Dune::make_shared<Dune::BlockVector<Dune::FieldVector<double, 1>>>(size);
auto normalStress = Dune::make_shared<SingletonVectorType>(size);
*normalStress = parset.get<double>("boundary.friction.normalstress");
std::string const friction_model =
parset.get<std::string>("boundary.friction.model");
if (friction_model == std::string("Ruina")) {
auto a = Dune::make_shared<Dune::BlockVector<Dune::FieldVector<double, 1>>>(
size);
auto a = Dune::make_shared<SingletonVectorType>(size);
*a = parset.get<double>("boundary.friction.ruina.a");
auto eta =
Dune::make_shared<Dune::BlockVector<Dune::FieldVector<double, 1>>>(
size);
auto eta = Dune::make_shared<SingletonVectorType>(size);
*eta = parset.get<double>("boundary.friction.eta");
myGlobalNonlinearity =
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment