This is the mail archive of the
gcc-bugs@gcc.gnu.org
mailing list for the GCC project.
[Bug middle-end/19974] New: Lapack hang in xeigtstc on ia-64
- From: "Thomas dot Koenig at online dot de" <gcc-bugzilla at gcc dot gnu dot org>
- To: gcc-bugs at gcc dot gnu dot org
- Date: 15 Feb 2005 16:18:20 -0000
- Subject: [Bug middle-end/19974] New: Lapack hang in xeigtstc on ia-64
- Reply-to: gcc-bugzilla at gcc dot gnu dot org
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