From 72cee6fa89b35fbd30f44da15d59daa402413ed0 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 10 Jul 2013 15:58:13 +0200
Subject: [PATCH] [Algorit] FIX: initial acceleration: do not ignore any DOFs
---
src/one-body-sample.cc | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 55caa315..b6cf7772 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -208,6 +208,7 @@ int main(int argc, char *argv[]) {
// Set up the boundary
Dune::BitSetVector<dims> ignoreNodes(finestSize, false);
+ Dune::BitSetVector<dims> noNodes(finestSize, false);
Dune::BitSetVector<1> neumannNodes(finestSize, false);
Dune::BitSetVector<1> frictionalNodes(finestSize, false);
@@ -370,7 +371,7 @@ int main(int argc, char *argv[]) {
}
BlockGSStep<MatrixType, VectorType> accelerationSolverStep(
massMatrix, a_initial, problem_rhs_initial);
- accelerationSolverStep.ignoreNodes_ = &ignoreNodes;
+ accelerationSolverStep.ignoreNodes_ = &noNodes;
LoopSolver<VectorType> accelerationSolver(
&accelerationSolverStep, 100000, // FIXME
1e-12, // FIXME
--
GitLab