LAPACK95 diagonalization problem

sknepnek — Fri, 02 Dec 2005 04:10:46 GMT

I have installed MKL 8.0 (evaluation version) on my Linux box running Mandrake 10.0 with ifort and icc versions 8.0

I compiled LAPACK95 interface by simply using

make PLAT=lnx32 lib

ifort complined about file source/cgelsd1.f90 line 30
USE MKL77_LAPACK, ONLY: ILAENV => MKL77_ILAENV

A simple change to
USE MKL77_LAPACK, ONLY: MKL77_ILAENV
plus change of line #76
SMLSIZ = ILAENV(9,'CGELSD',' ', 0,0,0,0)
to
SMLSIZ = MKL77_ILAENV(9,'CGELSD',' ', 0,0,0,0)

fixed the probelm, i.e., compilation went smoothy.

However, if I compile examples

make lib32

although the compilation goes witout any problems, the results are wrong.

For example, if you look into syev.res the output looks like:

SYEV Example Program Results
Matrix A:
-5.00000 -3.00000 1.00000 5.00000 3.00000
0.00000 2.00000 -1.00000 4.00000 -1.00000
0.00000 0.00000 6.00000 3.00000 -1.00000
0.00000 0.00000 0.00000 1.00000 -1.00000
0.00000 0.00000 0.00000 0.00000 -7.00000
CALL SYEV( A, W )
W on exit :
-11.24619
-6.24975
1.18524
5.00130
8.30939

Maple gives: 1, 2, -5, 6, -7

Simple inspection showes that -5 has to be an eigenvalue.

Thanks.

topic LAPACK95 diagonalization problem in Intel® oneAPI Math Kernel Library

LAPACK95 diagonalization problem