This is the mail archive of the
gcc-bugs@gcc.gnu.org
mailing list for the GCC project.
g77 (egcs-2.91.33 19980525) ICE in subst_stack_regs_pat
- To: egcs-bugs at cygnus dot com
- Subject: g77 (egcs-2.91.33 19980525) ICE in subst_stack_regs_pat
- From: Klaus-Georg Adams <Klaus-Georg dot Adams at chemie dot uni-karlsruhe dot de>
- Date: Wed, 27 May 1998 10:31:07 +0200
- Cc: fortran at gnu dot org
Hi,
g77 chokes on a file in the LAPACK distribution when optimizing.
-O0 compiles fine, -Os, -O[1-6] crash with:
/usr/local/egcs/bin/g77 -v -O1 -c cgtcon.f
Reading specs from /usr/local/egcs/lib/gcc-lib/i686-pc-linux-gnu/egcs-2.91.33/specs
gcc version egcs-2.91.33 19980525 (gcc2 ss-980502 experimental)
/usr/local/egcs/lib/gcc-lib/i686-pc-linux-gnu/egcs-2.91.33/f771 cgtcon.f -fset-g77-defaults -quiet -dumpbase cgtcon.f -O1 -version -fversion -o /tmp/cca32468.s
GNU F77 version egcs-2.91.33 19980525 (gcc2 ss-980502 experimental) (i686-pc-linux-gnu) compiled by GNU C version egcs-2.91.33 19980525 (gcc2 ss-980502 experimental).
GNU Fortran Front End version 0.5.22-release
/egcs/adams/egcs-cvs/gcc/reg-stack.c:2398: Internal compiler error in function subst_stack_regs_pat
The same thing happens on i686-pc-linux-gnulibc1.
Both compilers where produced with
configure --enable-shared --enable-haifa
and compiled with
make bootstrap BOOT_CFLAGS="-g -O4" CFLAGS="-g -O4"
The offending file is attached.
Let me know if you need more info.
-- kga
-------------------------------------------------------------------------
Klaus-Georg Adams Email: Klaus-Georg.Adams@chemie.uni-karlsruhe.de
Institut f. Anorg. Chemie, Lehrstuhl II Tel: 49(0)721 608 3485
Universität Karlsruhe, D-76128 Karlsruhe
-------------------------------------------------------------------------
SUBROUTINE CGTCON( NORM, N, DL, D, DU, DU2, IPIV, ANORM, RCOND,
$ WORK, INFO )
*
* -- LAPACK routine (version 2.0) --
* Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd.,
* Courant Institute, Argonne National Lab, and Rice University
* September 30, 1994
*
* .. Scalar Arguments ..
CHARACTER NORM
INTEGER INFO, N
REAL ANORM, RCOND
* ..
* .. Array Arguments ..
INTEGER IPIV( * )
COMPLEX D( * ), DL( * ), DU( * ), DU2( * ), WORK( * )
* ..
*
* Purpose
* =======
*
* CGTCON estimates the reciprocal of the condition number of a complex
* tridiagonal matrix A using the LU factorization as computed by
* CGTTRF.
*
* An estimate is obtained for norm(inv(A)), and the reciprocal of the
* condition number is computed as RCOND = 1 / (ANORM * norm(inv(A))).
*
* Arguments
* =========
*
* NORM (input) CHARACTER*1
* Specifies whether the 1-norm condition number or the
* infinity-norm condition number is required:
* = '1' or 'O': 1-norm;
* = 'I': Infinity-norm.
*
* N (input) INTEGER
* The order of the matrix A. N >= 0.
*
* DL (input) COMPLEX array, dimension (N-1)
* The (n-1) multipliers that define the matrix L from the
* LU factorization of A as computed by CGTTRF.
*
* D (input) COMPLEX array, dimension (N)
* The n diagonal elements of the upper triangular matrix U from
* the LU factorization of A.
*
* DU (input) COMPLEX array, dimension (N-1)
* The (n-1) elements of the first superdiagonal of U.
*
* DU2 (input) COMPLEX array, dimension (N-2)
* The (n-2) elements of the second superdiagonal of U.
*
* IPIV (input) INTEGER array, dimension (N)
* The pivot indices; for 1 <= i <= n, row i of the matrix was
* interchanged with row IPIV(i). IPIV(i) will always be either
* i or i+1; IPIV(i) = i indicates a row interchange was not
* required.
*
* ANORM (input) REAL
* If NORM = '1' or 'O', the 1-norm of the original matrix A.
* If NORM = 'I', the infinity-norm of the original matrix A.
*
* RCOND (output) REAL
* The reciprocal of the condition number of the matrix A,
* computed as RCOND = 1/(ANORM * AINVNM), where AINVNM is an
* estimate of the 1-norm of inv(A) computed in this routine.
*
* WORK (workspace) COMPLEX array, dimension (2*N)
*
* INFO (output) INTEGER
* = 0: successful exit
* < 0: if INFO = -i, the i-th argument had an illegal value
*
* =====================================================================
*
* .. Parameters ..
REAL ONE, ZERO
PARAMETER ( ONE = 1.0E+0, ZERO = 0.0E+0 )
* ..
* .. Local Scalars ..
LOGICAL ONENRM
INTEGER I, KASE, KASE1
REAL AINVNM
* ..
* .. External Functions ..
LOGICAL LSAME
EXTERNAL LSAME
* ..
* .. External Subroutines ..
EXTERNAL CGTTRS, CLACON, XERBLA
* ..
* .. Intrinsic Functions ..
INTRINSIC CMPLX
* ..
* .. Executable Statements ..
*
* Test the input arguments.
*
INFO = 0
ONENRM = NORM.EQ.'1' .OR. LSAME( NORM, 'O' )
IF( .NOT.ONENRM .AND. .NOT.LSAME( NORM, 'I' ) ) THEN
INFO = -1
ELSE IF( N.LT.0 ) THEN
INFO = -2
ELSE IF( ANORM.LT.ZERO ) THEN
INFO = -8
END IF
IF( INFO.NE.0 ) THEN
CALL XERBLA( 'CGTCON', -INFO )
RETURN
END IF
*
* Quick return if possible
*
RCOND = ZERO
IF( N.EQ.0 ) THEN
RCOND = ONE
RETURN
ELSE IF( ANORM.EQ.ZERO ) THEN
RETURN
END IF
*
* Check that D(1:N) is non-zero.
*
DO 10 I = 1, N
IF( D( I ).EQ.CMPLX( ZERO ) )
$ RETURN
10 CONTINUE
*
AINVNM = ZERO
IF( ONENRM ) THEN
KASE1 = 1
ELSE
KASE1 = 2
END IF
KASE = 0
20 CONTINUE
CALL CLACON( N, WORK( N+1 ), WORK, AINVNM, KASE )
IF( KASE.NE.0 ) THEN
IF( KASE.EQ.KASE1 ) THEN
*
* Multiply by inv(U)*inv(L).
*
CALL CGTTRS( 'No transpose', N, 1, DL, D, DU, DU2, IPIV,
$ WORK, N, INFO )
ELSE
*
* Multiply by inv(L')*inv(U').
*
CALL CGTTRS( 'Conjugate transpose', N, 1, DL, D, DU, DU2,
$ IPIV, WORK, N, INFO )
END IF
GO TO 20
END IF
*
* Compute the estimate of the reciprocal condition number.
*
IF( AINVNM.NE.ZERO )
$ RCOND = ( ONE / AINVNM ) / ANORM
*
RETURN
*
* End of CGTCON
*
END