Test a nicer problem

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;