[petsc-users] Bug report LMVM matrix class

[Stephan Köhler]

Dear PETSc/Tao team,

it seems to be that there is a bug in the LMVM matrix class:

The function MatMultAdd_LMVM, see, e.g., 
https://petsc.org/release/src/ksp/ksp/utils/lmvm/lmvmimpl.c.html at line 
114, if the vectors Y and Z are the same, then the result is wrong, 
since the first MatMult overwrites also the value in Y.

Best regards
Stephan Köhler

-- 
Stephan Köhler
TU Bergakademie Freiberg
Institut für numerische Mathematik und Optimierung

Akademiestraße 6
09599 Freiberg
Gebäudeteil Mittelbau, Zimmer 2.07

Telefon: +49 (0)3731 39-3173 (Büro)

-------------- next part --------------
A non-text attachment was scrubbed...
Name: OpenPGP_0xC9BF2C20DFE9F713.asc
Type: application/pgp-keys
Size: 758 bytes
Desc: OpenPGP public key
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20230509/d9d3a611/attachment.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: OpenPGP_signature
Type: application/pgp-signature
Size: 236 bytes
Desc: OpenPGP digital signature
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20230509/d9d3a611/attachment.sig>