Using Lapack++ 1.5.1
Affects: LaEigSolve, LaLUInverseIP, others?
Attempting to call LaEigSolve on a LaGenMatComplex which contains a NaN causes Lapack++ to hang if the NaN is not in the (1,1) position, or to throw an exception if NaN is in (1,1). It would make sense to me to throw an exception if there is a NaN anywhere in the matrix.
Using LUFactorizeIP and LaLUInverseIP to find the inverse of the same matrix returns an inverse of all NaN; I would expect that this should also throw an exception?
It's easiest just to demonstrate. bug.cpp and bug2.cpp show the two behaviors.