From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 5542 invoked by alias); 15 Feb 2005 16:18:28 -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 5473 invoked by uid 48); 15 Feb 2005 16:18:20 -0000 Date: Tue, 15 Feb 2005 21:20:00 -0000 From: "Thomas dot Koenig at online dot de" To: gcc-bugs@gcc.gnu.org Message-ID: <20050215161817.19974.Thomas.Koenig@online.de> Reply-To: gcc-bugzilla@gcc.gnu.org Subject: [Bug middle-end/19974] New: Lapack hang in xeigtstc on ia-64 X-Bugzilla-Reason: CC X-SW-Source: 2005-02/txt/msg01619.txt.bz2 List-Id: Checking Lapack with flag_complex_method = 2 on ia64-unknown-linux-gnu with -O0, I got a hang in xeigtstc when processing sep.in, in the routine slarrb.f. The compiler was the 20050213 snapshot. This may not have anything to do with the setting of flag_complex_method, because the routine in question does not have any complex variables. The code where it hangs looks like this: DO 70 J = 1, N - 1 TMP = D( J ) + S S = S*( LD( J ) / TMP )*L( J ) - RIGHT IF( TMP.LT.ZERO ) $ CNT = CNT + 1 70 CONTINUE Here's a sample gdb session: $ gdb xeigtstc GNU gdb Red Hat Linux (6.1post-1.20040607.52rh) Copyright 2004 Free Software Foundation, Inc. GDB is free software, covered by the GNU General Public License, and you are welcome to change it and/or distribute copies of it under certain conditions. Type "show copying" to see the conditions. There is absolutely no warranty for GDB. Type "show warranty" for details. This GDB was configured as "ia64-redhat-linux-gnu"...Using host libthread_db library "/lib/tls/libthread_db.so.1". (gdb) r Starting program: /home/zfkts/zfkts/LAPACK-gfortran/TESTING/xeigtstc Program received signal SIGINT, Interrupt. 0x200000000037bbc1 in read () from /lib/tls/libc.so.6.1 (gdb) r < sep.in The program being debugged has been started already. Start it from the beginning? (y or n) y Starting program: /home/zfkts/zfkts/LAPACK-gfortran/TESTING/xeigtstc < sep.in Tests of the Hermitian Eigenvalue Problem routines LAPACK VERSION 3.0, released June 30, 1999 The following parameter values will be used: M: 0 1 2 3 5 20 N: 0 1 2 3 5 20 NB: 1 3 3 3 10 NBMIN: 2 2 2 2 2 NX: 1 0 5 9 1 Relative machine underflow is taken to be 0.117549E-37 Relative machine overflow is taken to be 0.340282E+39 Relative machine precision is taken to be 0.596046E-07 Routines pass computational tests if test ratio is less than 50.00 CST routines passed the tests of the error exits (114 tests done) SEP: NB = 1, NBMIN = 2, NX = 1 All tests for CST passed the threshold ( 4662 tests run) All tests for CST drivers passed the threshold ( 11664 tests run) SEP: NB = 3, NBMIN = 2, NX = 0 CST -- Complex Hermitian eigenvalue problem Matrix types (see CCHKST for details): Special Matrices: 1=Zero matrix. 5=Diagonal: clustered entries. 2=Identity matrix. 6=Diagonal: large, evenly spaced. 3=Diagonal: evenly spaced entries. 7=Diagonal: small, evenly spaced. 4=Diagonal: geometr. spaced entries. Dense Hermitian Matrices: 8=Evenly spaced eigenvals. 12=Small, evenly spaced eigenvals. 9=Geometrically spaced eigenvals. 13=Matrix with random O(1) entries. 10=Clustered eigenvalues. 14=Matrix with large random entries. 11=Large, evenly spaced eigenvals. 15=Matrix with small random entries. 16=Positive definite, evenly spaced eigenvalues 17=Positive definite, geometrically spaced eigenvlaues 18=Positive definite, clustered eigenvalues 19=Positive definite, small evenly spaced eigenvalues 20=Positive definite, large evenly spaced eigenvalues 21=Diagonally dominant tridiagonal, geometrically spaced eigenvalues Test performed: see CCHKST for details. Matrix order= 5, type= 9, seed= 894,3587, 994, 217, result 36 is 62.23 CST: 1 out of 4662 tests failed to pass the threshold All tests for CST drivers passed the threshold ( 11664 tests run) SEP: NB = 3, NBMIN = 2, NX = 5 All tests for CST passed the threshold ( 4662 tests run) All tests for CST drivers passed the threshold ( 11664 tests run) SEP: NB = 3, NBMIN = 2, NX = 9 CST -- Complex Hermitian eigenvalue problem Matrix types (see CCHKST for details): Special Matrices: 1=Zero matrix. 5=Diagonal: clustered entries. 2=Identity matrix. 6=Diagonal: large, evenly spaced. 3=Diagonal: evenly spaced entries. 7=Diagonal: small, evenly spaced. 4=Diagonal: geometr. spaced entries. Dense Hermitian Matrices: 8=Evenly spaced eigenvals. 12=Small, evenly spaced eigenvals. 9=Geometrically spaced eigenvals. 13=Matrix with random O(1) entries. 10=Clustered eigenvalues. 14=Matrix with large random entries. 11=Large, evenly spaced eigenvals. 15=Matrix with small random entries. 16=Positive definite, evenly spaced eigenvalues 17=Positive definite, geometrically spaced eigenvlaues 18=Positive definite, clustered eigenvalues 19=Positive definite, small evenly spaced eigenvalues 20=Positive definite, large evenly spaced eigenvalues 21=Diagonally dominant tridiagonal, geometrically spaced eigenvalues Test performed: see CCHKST for details. Matrix order= 20, type= 9, seed=1052,3651,3662,3633, result 35 is 69.27 Matrix order= 20, type= 9, seed=1052,3651,3662,3633, result 36 is 2182.68 CST: 2 out of 4662 tests failed to pass the threshold CST -- Complex Hermitian eigenvalue problem Matrix types (see xDRVST for details): Special Matrices: 1=Zero matrix. 5=Diagonal: clustered entries. 2=Identity matrix. 6=Diagonal: large, evenly spaced. 3=Diagonal: evenly spaced entries. 7=Diagonal: small, evenly spaced. 4=Diagonal: geometr. spaced entries. Dense Hermitian Matrices: 8=Evenly spaced eigenvals. 12=Small, evenly spaced eigenvals. 9=Geometrically spaced eigenvals. 13=Matrix with random O(1) entries. 10=Clustered eigenvalues. 14=Matrix with large random entries. 11=Large, evenly spaced eigenvals. 15=Matrix with small random entries. Tests performed: See cdrvst.f Matrix order= 20, type= 9, seed=1494,3156,1807,2209, result 101 is 353.18 (wait for a long time...) Program received signal SIGINT, Interrupt. slarrb_ (n=@0x4, d=0x4, l=0x4, ld=0x4, lld=0x4, ifirst=@0x4, ilast=@0x4, sigma=@0x4, reltol=Cannot access memory at address 0x1 ) at slarrb.f:191 191 TMP = D( N ) + S Current language: auto; currently fortran (gdb) c Continuing. Program received signal SIGINT, Interrupt. slarrb_ (n=@0x7fc000007f800000, d=0x7fc000007f800000, l=0x7fc000007f800000, ld=0x7fc000007f800000, lld=0x7fc000007f800000, ifirst=@0x7fc000007f800000, ilast=@0x7fc000007f800000, sigma=@0x7fc000007f800000, reltol=Cannot access memory at address 0x0 ) at slarrb.f:195 195 DELTA = TWO*DELTA (gdb) c Continuing. Program received signal SIGINT, Interrupt. slarrb_ (n=Cannot access memory at address 0x1 ) at slarrb.f:185 185 DO 70 J = 1, N - 1 (gdb) c Continuing. Program received signal SIGINT, Interrupt. slarrb_ (n=@0x47fc00000, d=0x47fc00000, l=0x47fc00000, ld=0x47fc00000, lld=0x47fc00000, ifirst=@0x47fc00000, ilast=@0x47fc00000, sigma=@0x47fc00000, reltol=Cannot access memory at address 0x0 ) at slarrb.f:183 183 S = -RIGHT (gdb) c Continuing. Program received signal SIGINT, Interrupt. 0x400000000047fd71 in slarrb_ (n=@0x47fc00000, d=0x47fc00000, l=0x47fc00000, ld=0x47fc00000, lld=0x47fc00000, ifirst=@0x47fc00000, ilast=@0x47fc00000, sigma=@0x47fc00000, reltol=Cannot access memory at address 0x0 ) at slarrb.f:183 183 S = -RIGHT (gdb) c Continuing. Program received signal SIGINT, Interrupt. 0x400000000047fe91 in slarrb_ (n=Cannot access memory at address 0x3 ) at slarrb.f:187 187 S = S*( LD( J ) / TMP )*L( J ) - RIGHT (gdb) c Continuing. Program received signal SIGINT, Interrupt. 0x4000000000480111 in slarrb_ (n=@0x7fc000007fc00000, d=0x7fc000007fc00000, l=0x7fc000007fc00000, ld=0x7fc000007fc00000, lld=0x7fc000007fc00000, ifirst=@0x7fc000007fc00000, ilast=@0x7fc000007fc00000, sigma=@0x7fc000007fc00000, reltol=Cannot access memory at address 0x3 ) at slarrb.f:191 191 TMP = D( N ) + S (gdb) c Continuing. Program received signal SIGINT, Interrupt. 0x400000000047fe70 in slarrb_ (n=@0x7fc000007fc00000, d=0x7fc000007fc00000, l=0x7fc000007fc00000, ld=0x7fc000007fc00000, lld=0x7fc000007fc00000, ifirst=@0x7fc000007fc00000, ilast=@0x7fc000007fc00000, sigma=@0x7fc000007fc00000, reltol=Cannot access memory at address 0x1 ) at slarrb.f:186 186 TMP = D( J ) + S (gdb) c Continuing. Program received signal SIGINT, Interrupt. slarrb_ (n=@0x60000000000d5cf8, d=0x60000000000d5cf8, l=0x60000000000d5cf8, ld=0x60000000000d5cf8, lld=0x60000000000d5cf8, ifirst=@0x60000000000d5cf8, ilast=@0x60000000000d5cf8, sigma=@0x60000000000d5cf8, reltol=Cannot access memory at address 0x1 ) at slarrb.f:187 187 S = S*( LD( J ) / TMP )*L( J ) - RIGHT -- Summary: Lapack hang in xeigtstc on ia-64 Product: gcc Version: unknown Status: UNCONFIRMED Severity: normal Priority: P2 Component: middle-end AssignedTo: unassigned at gcc dot gnu dot org ReportedBy: Thomas dot Koenig at online dot de CC: gcc-bugs at gcc dot gnu dot org GCC target triplet: ia64-unknown-linux-gnu http://gcc.gnu.org/bugzilla/show_bug.cgi?id=19974