Eigensolvers: check and make sure seeding vectors don't have nan or info norms

lattice / quda

QUDA is a library for performing calculations in lattice QCD on GPUs.

https://lattice.github.io/quda

Other

279 stars 94 forks source link

Eigensolvers: check and make sure seeding vectors don't have nan or info norms #1464

Closed weinbe2 closed 2 months ago

weinbe2 commented 2 months ago

By default, the routine that prepared initial guesses for QUDA's eigensolvers only checked to make sure that the initial vector had a non-zero norm. This extends it to checking for nan and inf as well. Reported by @leonhostetler .

Don't merge until I get confirmation from Leon that this works, the initial "fix" I gave him modified interface_quda.cpp instead of checking deeper down in the eigensolver routines.

weinbe2 commented 2 months ago

Tentatively approving subject to confirmation from @leonhostetler that this fixes the issue.

He confirmed via our e-mail chain, so this should be good to go.