|Re: [eigen] problem with MKL bindings for dense Eigen|
[ Thread Index |
| More lists.tuxfamily.org/eigen Archives
Am 06.04.2018 um 22:49 schrieb Gael Guennebaud:
> Just a guess, but maybe it's just that MKL gives you different values for the lowest eigenvalues (e.g., 1e-16 vs 1e-30 vs exact 0 while the largest is 1) but your code is lacking some thresholds to properly deal with them (e.g.,
> to ignore them). So tracking the largest and the smallest eigenvalues might a be a good start.
It's all about tracking the lowest eigen values and that's done properly, at least
I intended to do it properly.
That the MKL BLAS version is giving slightly different results is fine.
E.g. if FMA are handled differently it will already influence the result.
My concern is about the mess I get from the MKL diagonalization.
I use MKL a lot in my main code (that one is still without Eigen) and there
everything is fine.