From a845c2b10da6242bd4594049dca4b8d025db1427 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Tue, 9 Jul 2013 16:41:35 +0200
Subject: [PATCH] [Algorit] Initial iterate: 0
If we use the old velocity, it can happen that we hit the kernel of
the singular norm component
---
src/timestepping.cc | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/timestepping.cc b/src/timestepping.cc
index bdadfefc..fe6b09e2 100644
--- a/src/timestepping.cc
+++ b/src/timestepping.cc
@@ -139,7 +139,7 @@ void Newmark<VectorType, MatrixType, FunctionType, dim>::setup(
Arithmetic::addProduct(problem_AB, 2.0 / tau, B);
// v_old makes a good initial iterate; we could use anything, though
- problem_iterate = v_old;
+ problem_iterate = 0.0;
for (size_t i = 0; i < dirichletNodes.size(); ++i)
switch (dirichletNodes[i].count()) {
--
GitLab