Merge branch 'issue/dirichletBC' into 'develop'

replaced symmetric DBC with unsymmetric one, since this seems to be wrong. Needs…

See merge request spraetor/dune-amdis!23

examples/init/stokes.dat.2d

 dimension of world: 2
 
-stokesMesh->macro file name: ./macro/macro.stand.2d
 stokesMesh->global refinements: 0
+stokesMesh->dimension: 1.0 1.0
+stokesMesh->num cells: 4 4
 
 stokes->mesh: stokesMesh
 
@@ -14,7 +15,7 @@ stokes->solver->info:	2
 
 stokes->output[0]->filename: stokes_u.2d
 stokes->output[0]->name: u
-stokes->output[0]->subsampling: 0
+stokes->output[0]->subsampling: 2
 stokes->output[0]->output directory: ./output
 stokes->output[0]->ParaView format: 1
 stokes->output[0]->ParaView mode: 1
@@ -25,7 +26,7 @@ stokes->output[1]->output directory: ./output
 stokes->output[1]->ParaView format: 1
 stokes->output[1]->ParaView mode: 1
 
-stokes->viscosity: 0.1
+stokes->viscosity: 1.0
 stokes->density: 1.0
 stokes->boundary velocity: 10.0
 

src/amdis/DirichletBC.hpp

@@ -78,8 +78,8 @@ namespace AMDiS
       finishImpl(matrix, rhs, solution, rowBasis, colBasis, ValueCategory_t<Range>{});
 
       // subtract columns of dirichlet nodes from rhs
-      for (auto const& triplet : columns)
-        rhs[triplet.row] -= triplet.value * solution[triplet.col];
+      // for (auto const& triplet : columns)
+      //   rhs[triplet.row] -= triplet.value * solution[triplet.col];
     }
 
   protected:

src/amdis/linear_algebra/mtl/DOFMatrix.hpp

@@ -181,6 +181,7 @@ namespace AMDiS
       auto value = mtl::mat::value_map(*matrix);
 
       // iterate over the matrix
+#if 0
       for (auto r : mtl::major_of(*matrix)) {   // rows or columns
         for (auto i : mtl::nz_of(r)) {          // non-zeros within
           if (dirichletNodes[row(i)]) {
@@ -193,6 +194,16 @@ namespace AMDiS
           }
         }
       }
+#endif
+
+      for (auto r : mtl::rows_of(*matrix)) {   // rows or columns
+        if (dirichletNodes[r.value()]) {
+          for (auto i : mtl::nz_of(r)) {          // non-zeros within
+            // set identity row
+            value(i, (row(i) == col(i) ? value_type(1) : value_type(0)) );
+          }
+        }
+      }
 
       return columns;
     }