diff --git a/dune/solvers/solvers/umfpacksolver.hh b/dune/solvers/solvers/umfpacksolver.hh index 90c43c3fc4c11c868fc02084c4e7e8164922a27d..d98a7dac89e38b31d2ea5139e4d20351028b1c4d 100644 --- a/dune/solvers/solvers/umfpacksolver.hh +++ b/dune/solvers/solvers/umfpacksolver.hh @@ -90,10 +90,6 @@ public: // Modify Dirichlet rows /////////////////////////////////////////// - for (size_t j=0; j<modifiedStiffness.N(); j++) - for (int k=0; k<blocksize; k++) - modifiedStiffness[j][j][k][k] += 0.0; - for (size_t j=0; j<modifiedStiffness.N(); j++) { auto cIt = modifiedStiffness[j].begin();