diff --git a/dune/tectonic/frictionpotential.hh b/dune/tectonic/frictionpotential.hh index de2d1dc53fc5d924700e166f4de189478241a370..89f29e5196ece66e9e961097b3e20262dc82e090 100644 --- a/dune/tectonic/frictionpotential.hh +++ b/dune/tectonic/frictionpotential.hh @@ -102,7 +102,7 @@ class RegularisedRateState : public FrictionPotential { double Vmin; }; -class TrivialFunction : public FrictionPotential { +class ZeroFunction : public FrictionPotential { public: double evaluate(double) const { return 0; } diff --git a/dune/tectonic/globalratestatefriction.hh b/dune/tectonic/globalratestatefriction.hh index a43619c1b59fd38b07242cc2cf0f4fba24ea5735..3edfe081d3ae5e624d6eda435ad63e9ef8fe4418 100644 --- a/dune/tectonic/globalratestatefriction.hh +++ b/dune/tectonic/globalratestatefriction.hh @@ -29,8 +29,8 @@ class GlobalRateStateFriction : public GlobalFriction<Matrix, Vector> { ScalarVector const &weights, ScalarVector const &weightedNormalStress) : restrictions(weightedNormalStress.size()) { - auto trivialNonlinearity = - std::make_shared<Friction>(std::make_shared<TrivialFunction>()); + auto zeroNonlinearity = + std::make_shared<Friction>(std::make_shared<ZeroFunction>()); Dune::MultipleCodimMultipleGeomTypeMapper< GridView, Dune::MCMGVertexLayout> const vertexMapper(gridView); @@ -39,7 +39,7 @@ class GlobalRateStateFriction : public GlobalFriction<Matrix, Vector> { auto const i = vertexMapper.index(*it); auto const coordinate = it->geometry().corner(0); if (not frictionalBoundary.containsVertex(i)) { - restrictions[i] = trivialNonlinearity; + restrictions[i] = zeroNonlinearity; continue; } auto const fp = std::make_shared<ScalarFriction>(