Changes
Page history
Update home
authored
May 14, 2017
by
Praetorius, Simon
Show whitespace changes
Inline
Side-by-side
home.md
View page @
1693e29d
...
@@ -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"
);
```
```
...
...
...
...