Eigen/src/Eigenvalues/EigenSolver.h should pass Scalar typename to cdiv explicitly
Submitted by Thad Hughes
Assigned to Nobody
Link to original bugzilla bug (#636)
Description
For example, I think you should write this:
std::complex<Scalar> cc = cdiv<Scalar>(-ra,-sa,w,q);
Instead of this:
std::complex<Scalar> cc = cdiv(-ra,-sa,w,q);
There is one place in the file where it is done correctly, other places where it is not.
The reason is that if Scalar is a type that uses expression templates -ra is not necessarily the same type as w, so the type inference when invoking cdiv does not work correctly.
Edited by Eigen Bugzilla