I narrowed the problem down to the magma_stedx call within the zheevd_gpu.cpp file.
More precicely, I modified the magma_zheevd function to use lapackf77_zstemr instead of magma_zstedx.
This modification suffices to converge to the right solution.
What is confusing me, hower, is that I'm not able to build a selfcontained example of this bug (i.e. a separate binary which reads all the input arguments from a file).
My selfcontained example allways converges for exactly the same input arguments and by using the original MAGMA version. However,
within our application the original MAGMA version is never converging.
I'm also linking my selfcontained binary against some (not all) of the libraries which we are linking against in our application as well.
Can you think of any other problems that I'm running into?
Do you think that this is a MAGMA issue or an issue with the Bisection Algorithm?
Last edited by Paul2822
on Tue Sep 24, 2013 2:01 pm, edited 1 time in total.