#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.