I'm trying to use the CLAPACK package, especially zgesvd to do sv decomposition, but found that the zgesvd_ function hangs. When debugging the library, I found that it was actually the dlamch_ (dlamch1 or dlamch2 actually) function that stucks. When searching for solutions, I found that the same problem occured in the Fortran version on certain platforms (among x86), and that compiling with the -ffloat-store would solve the problem. That did indeed solve the problem.
Would it be an idea to compile the library with -ffloat-store as default?
My platform is a dual core Intel x86 machine, gcc --version: "gcc (GCC) 4.1.2 20070626 (Red Hat 4.1.2-14)" and CLAPACK-3.1.1 with the internally provided blas library.