From f9ad0c192c2ecc7ca77f45584127ffe18936246a Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Fri, 28 Oct 2011 18:53:55 +0200
Subject: [PATCH] Test a nicer problem
---
src/one-body-sample.cc | 15 ++++++---------
1 file changed, 6 insertions(+), 9 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index fe8bb656..066d9dfe 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -47,10 +47,10 @@ int main() {
typedef Dune::BlockVector<SmallVector> VectorType;
// FIXME: Random values
- double const E = 1e8;
+ double const E = 1e3;
double const nu = 0.3;
- int const refinements = 1;
- size_t const solver_maxIterations = 10;
+ int const refinements = 5;
+ size_t const solver_maxIterations = 100;
double const solver_tolerance = 1e-4;
typedef Dune::YaspGrid<dim> GridType;
@@ -81,13 +81,10 @@ int main() {
// Fill right-hand side with semi-random numbers
for (size_t i = 0; i < f.size(); ++i)
for (size_t j = 0; j < dim; ++j)
- f[i][j] = i + j;
+ f[i][j] = 300;
VectorType u1(grid.size(grid.maxLevel(), dim));
- // Fill initial guess with other semi-random numbers
- for (size_t i = 0; i < f.size(); ++i)
- for (size_t j = 0; j < dim; ++j)
- u1[i][j] = 17 * (j + 1);
+ u1 = 0;
VectorType u2 = u1;
Dune::BitSetVector<VectorType::block_type::dimension> ignoreNodes(
@@ -126,7 +123,7 @@ int main() {
++bounding_nodes;
size_t id = myVertexMapper.map(*it);
// std::cout << "Ignoring id #" << id << std::endl;
- ignoreNodes[id] = true;
+ // ignoreNodes[id] = true;
}
if (extremal) {
++extremal_nodes;
--
GitLab