... | ... | @@ -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");
|
|
|
```
|
... | ... | |