diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index fe8bb656be1bf0078574f1c576e440500f59c1d8..066d9dfe726187cfaf5c9cd59a16882955f5ea9f 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;