From 2aaf4029a44b8eedacabddb8a415e18c4d768961 Mon Sep 17 00:00:00 2001 From: podlesny <podlesny@zedat.fu-berlin.de> Date: Sun, 14 Mar 2021 17:14:28 +0100 Subject: [PATCH] set dirichlet data for dirichlet body correctly --- dune/tectonic/data-structures/program_state.hh | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/dune/tectonic/data-structures/program_state.hh b/dune/tectonic/data-structures/program_state.hh index 6c6f3875..b4092850 100644 --- a/dune/tectonic/data-structures/program_state.hh +++ b/dune/tectonic/data-structures/program_state.hh @@ -159,8 +159,9 @@ template <class VectorTEMPLATE, class ScalarVectorTEMPLATE> class ProgramState { const auto& body = *contactNetwork.body(i); if (body.data()->getYoungModulus() == 0.0) { - for (; dof<body.nVertices(); dof++) { + for (size_t j=0; j<body.nVertices(); j++) { dirichletNodes[dof] = true; + dof++; } } else { dof += body.nVertices(); -- GitLab