Fix complex conjugate in MGS dot product

ginkgo-project / ginkgo

Numerical linear algebra software package

https://ginkgo-project.github.io/

BSD 3-Clause "New" or "Revised" License

415 stars 90 forks source link

Fix complex conjugate in MGS dot product #1678

Closed nbeams closed 2 months ago

nbeams commented 2 months ago

This fixes the discrepancy between MGS and CGS in GMRES convergence for complex numbers, since compute_conj_dot takes the conjugate of the calling vector, rather than the argument.