From 870f6e7ebdc6311435bc56bf25695d0f9f7aceba Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 4 Jan 2012 16:42:17 +0100
Subject: [PATCH] Start from zero again
Since truncation works properly now, this hack is no longer necessary;
it actually does more harm than good
---
src/one-body-sample.cc | 23 -----------------------
1 file changed, 23 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 40425882..12429526 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -342,29 +342,6 @@ int main(int argc, char *argv[]) {
MyConvexProblemType myConvexProblem(stiffnessMatrix,
*myGlobalNonlinearity, b4, u4_diff);
- // FIXME: This is sort of nasty
- // Choose an initial iterate that is admissable
- // and non-zero in as many places as possible
- if (run == 1)
- for (size_t i = 0; i < ignoreNodes.size(); ++i) {
- size_t const count = ignoreNodes[i].count();
-
- if (count == dim)
- continue;
-
- if (count == 0) {
- u4_diff[i] = 1e-4;
- continue;
- }
-
- assert(count == 1);
- assert(ignoreNodes[i][1]); // Only the y coordinate is allowed to be
- // fixed
- for (int j = 0; j < dim; ++j)
- if (j != 1)
- u4_diff[i][j] = 1e-4;
- }
-
MyBlockProblemType myBlockProblem(parset, myConvexProblem);
multigridStep.setProblem(u4_diff, myBlockProblem);
--
GitLab