The full source is available by anonymous ftp from:
ftp://ftp.cs.unm.edu/pub/smith-quetzal/codes.tar.Z
Directory lu (or I can provide directory lu, with gfc makefile)
The programme hangs in this routine, unless the print * , "" indicated below is
included.
SUBROUTINE SLAMC1( BETA, T, RND, IEEE1 )
*
* -- LAPACK auxiliary routine (version 1.0) --
* Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd.,
* Courant Institute, Argonne National Lab, and Rice University
* February 29, 1992
etc., etc. until
*
*+ WHILE( C.EQ.ONE )LOOP
10 CONTINUE
IF( C.EQ.ONE ) THEN
A = 2*A
C = SLAMC3( A, ONE )
C = SLAMC3( C, -A )
GO TO 10
END IF
*+ END WHILE
*
* Now compute b = 2.0**m with the smallest positive integer m
* such that
*
* fl( a + b ) .gt. a.
*
B = 1
C = SLAMC3( A, B )
print *,"" !!!!!The code hangs without this line
!!!!!Putting C = C does not cure the problem
*
*+ WHILE( C.EQ.A )LOOP
20 CONTINUE
IF( C.EQ.A ) THEN
B = 2*B
C = SLAMC3( A, B )
print *, C , C.EQ.A
GO TO 20
END IF
*+ END WHILE
*
and so on......
--
Summary: Runtime hang in LAPACK routine SLAMC1 - in Quetzal
benchmark suite
Product: gcc
Version: 4.0.0
Status: UNCONFIRMED
Severity: normal
Priority: P2
Component: fortran
AssignedTo: unassigned at gcc dot gnu dot org
ReportedBy: paul dot richard dot thomas at cea dot fr
CC: gcc-bugs at gcc dot gnu dot org,paulthomas2 at wanadoo
dot fr
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=18714
