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 ...@@ -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; Eigen::UmfPackLU<typename Matrix::Matrix> solver;
DOFVector<double> u(gv, 0); 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)); 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 ...@@ -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`: For the visualization, we write the `DOFVector` in a VTK-format, that can be plotted using `ParaView`:
```cpp ```cpp
Dune::VTKWriter<typename GridBase::LeafGridView> writer(gridBase->leafGridView()); Dune::VTKWriter<GridView> writer(gv);
writer.addVertexData(u, "u"); writer.addVertexData(u, "u");
writer.write("poisson"); writer.write("poisson");
``` ```
... ...
......