Update home authored by Praetorius, Simon's avatar Praetorius, Simon
......@@ -93,7 +93,8 @@ The resulting linear system `A*u=f` can now be solved using any linear-solver yo
Eigen::UmfPackLU<typename Matrix::Matrix> solver;
DOFVector<double> u(gv, 0);
u = solver.compute(A).solve(f);
solver.compute(A);
u = solver.solve(f);
msg("|A*u-f| = ", residuum(A, u, f));
```
......@@ -102,7 +103,7 @@ The resulting linear system `A*u=f` can now be solved using any linear-solver yo
For the visualization, we write the `DOFVector` in a VTK-format, that can be plotted using `ParaView`:
```cpp
Dune::VTKWriter<typename GridBase::LeafGridView> writer(gridBase->leafGridView());
Dune::VTKWriter<GridView> writer(gv);
writer.addVertexData(u, "u");
writer.write("poisson");
```
......
......