Problems installing LAPACK in ALPHA-DEC

[wd20lp]

Dear All,

We are trying to compile LAPACK in a DEC ALPHA machine and
we are getting some errors messages (please, see below). Does
anybody have suggestions on how to overcome this problem?

Cheers,

Luciano

---------------------------------------------------------
.....
f77 -04 -fpe0 -c xlaenv.f
f77 -04 -fpe0 -c chkxer.f
f77 -o ../xeigtsts schkee.o sbdt01.o sbdt02.o sbdt03.o schkbb.o schkbd.o
schkbk.o schkbl.o schkec.o schkgg.o schkgk.o schkgl.o schkhs.o schksb.o
schkst.o sckglm.o sckgqr.o sckgsv.o scklse.o sdrges.o sdrgev.o sdrgsx.o
sdrgvx.o sdrvbd.o sdrves.o sdrvev.o sdrvgg.o sdrvsg.o sdrvst.o sdrvsx.o
sdrvvx.o serrbd.o serrec.o serred.o serrgg.o serrhs.o serrst.o sget02.o
sget10.o sget22.o sget23.o sget24.o sget31.o sget32.o sget33.o sget34.o
sget35.o sget36.o sget37.o sget38.o sget39.o sget51.o sget52.o sget53.o
sget54.o sglmts.o sgqrts.o sgrqts.o sgsvts.o shst01.o slarfy.o slarhs.o
slatm4.o slctes.o slctsx.o slsets.o sort01.o sort03.o ssbt21.o ssgt01.o
sslect.o sspt21.o sstt21.o sstt22.o ssyt21.o ssyt22.o slafts.o slahd2.o
slasum.o slatb9.o sstech.o sstect.o ssvdch.o ssvdct.o ssxt1.o alahdg.o
alasum.o alasvm.o alareq.o ilaenv.o xerbla.o xlaenv.o chkxer.o
.../../tmglib_alpha.a ../../lapack_alpha.a -ldxml
NEP: Testing Nonsymmetric Eigenvalue Problem routines
../xeigtsts < nep.in > snep.out 2>&1
SEP: Testing Symmetric Eigenvalue Problem routines
../xeigtsts < sep.in > ssep.out 2>&1
sh: 85269 Abort - core dumped
*** Exit 134
Stop.
*** Exit 1
Stop.

[Julie]

Luciano,

Please take a look at http://icl.cs.utk.edu/lapack-forum/viewtopic.php?t=78

Hope it helps.

Julie

[dmday]

I observed similar problems using a two x86 four-core processors, Goto's threaded BLAS library (1.19) with lapack3.1.1. At first I used gfortran version 4.1.1, and I worked really hard to get just stest to pass. Then I upgraded gfortran to version 4.2.2., and the tests ran.

I digress to mention the steps needed for stest to pass with gfortran 4.1.1. I hacked LAPACK in several ways. In schkee.f I to avoid a segmentation fault I changed
$ RESULT( 500 ), TAUA( NMAX ), TAUB( NMAX ),
to
$ RESULT( 5000 ), TAUA( NMAX ), TAUB( NMAX ),
There were several problems with Infinity and NaN's. I added a redundant test for NaNs to slange.f. slascl.f and slartg.f were being called with infinite arguments, and were stuck in infinite loops. So again I added code to detect such (illegal) input and handle it accordingly.
And here is the make.inc (minus comments).

Code: Select all

SHELL = /bin/sh
PLAT = _LINUX
FORTRAN  = gfortran
OPTS     = -funroll-all-loops -frounding-math -fsignaling-nans -fno-range-check -O3
DRVOPTS  = $(OPTS)
NOOPT    =
LOADER   = gfortran
LOADOPTS = -lpthread
TIMER    = INT_CPU_TIME
ARCH     = ar
ARCHFLAGS= cr
RANLIB   = ranlib
BLASLIB      = /home/dmday/GotoBLAS/libgoto.a -lpthread
LAPACKLIB    = lapack$(PLAT).a
TMGLIB       = tmglib$(PLAT).a
EIGSRCLIB    = eigsrc$(PLAT).a
LINSRCLIB    = linsrc$(PLAT).a


--David Day