From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 22233 invoked by alias); 29 Nov 2004 10:57:51 -0000 Mailing-List: contact gcc-bugs-help@gcc.gnu.org; run by ezmlm Precedence: bulk List-Archive: List-Post: List-Help: Sender: gcc-bugs-owner@gcc.gnu.org Received: (qmail 21906 invoked by uid 48); 29 Nov 2004 10:57:37 -0000 Date: Mon, 29 Nov 2004 10:57:00 -0000 From: "paul dot richard dot thomas at cea dot fr" To: gcc-bugs@gcc.gnu.org Message-ID: <20041129105735.18714.paul.richard.thomas@cea.fr> Reply-To: gcc-bugzilla@gcc.gnu.org Subject: [Bug fortran/18714] New: Runtime hang in LAPACK routine SLAMC1 - in Quetzal benchmark suite X-Bugzilla-Reason: CC X-SW-Source: 2004-11/txt/msg03574.txt.bz2 List-Id: 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