[Users] Cactus-built LAPACK, BLAS giving incorrect results

Kelly, Bernard J. (GSFC-660.0)[UNIVERSITY OF MARYLAND BALTIMORE COUNTY] bernard.j.kelly at nasa.gov
Fri May 27 11:09:06 CDT 2011

Hi all. I recently needed to use a LAPACK routine in a Cactus thorn. Since my local machine's installation of LAPACK & BLAS was old, I decided to force an installation by including the thorns "ExternalLibraries/[BLAS,LAPACK]". This seemed to work, with the resulting .a files ending up in configs/<my config name>/scratch/external/LAPACK, etc. Compiling and running my cactus executable, however, I found the LAPACK gave completely wrong results.

For a much simpler test case that (see attached F90 standalone program), consider the 2x2 system A x = b, where A = (3, 2\\2, 6), and b = (5, 8).

Correct result for "x" is (1,1).

Cactus-forced LAPACK (using routine "dgesv") yields (-2,1).

I later installed LAPACK (v. 3.3.0) + BLAS by hand on my machine. Linking against these yields the correct answer.

I'm on a 32-bit Linux machine, using the Intel 10 compilers.

I'd appreciate any insights on this. Thanks,

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/users/attachments/20110527/5290b805/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: lapack_test.F90
Type: application/octet-stream
Size: 677 bytes
Desc: lapack_test.F90
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20110527/5290b805/attachment.obj 

More information about the Users mailing list