#397 Preconditioner< X, Y >::apply documentation-implementation discrepancy
Metadata
| Property | Value |
|---|---|
| Reported by | Uli Sack (usack@math.fu-berlin.de) |
| Reported at | May 20, 2008 10:14 |
| Type | Bug Report |
| Version | Git (pre2.4) [autotools] |
| Operating System | Linux |
| Last edited by | Markus Blatt (markus@dr-blatt.de) |
| Last edited at | Jul 25, 2008 12:25 |
| Closed by | Markus Blatt (markus@dr-blatt.de) |
| Closed at | Jul 25, 2008 12:25 |
| Closed in version | Unknown |
| Resolution | Works for me |
| Comment |
Description
The documentation of
Dune::Preconditioner< X, Y >::apply (X& v, const Y&d)
specifies that on exit v=M^{-1}d. This, however, is not implemented as such - at least in the Jacobi Preconditioner SeqJac. If v already contains old data it is not overwritten but used in the dbjac method.
I would therefore suggest either adjusting documentation or setting v=0.0; in the beginning of apply(v,d) to match the documentation/specification.