This is the mail archive of the
gcc@gcc.gnu.org
mailing list for the GCC project.
RE: Advance warning: g77-2.96 and LAPACK 3.0.
- To: "'Toon Moene'" <toon at moene dot indiv dot nluug dot nl>, gcc at gcc dot gnu dot org
- Subject: RE: Advance warning: g77-2.96 and LAPACK 3.0.
- From: "Billinghurst, David (CRTS)" <David dot Billinghurst at riotinto dot com>
- Date: Mon, 4 Sep 2000 03:26:28 -0000
Conversely, I have seen progressive improvement in LAPACK test results over
the last year on Irix6.5
* no longer need to emulate complex arithmetic
* -mabi=64 works (although there is one ICE)
* a couple of bugs identified then fixed
* can run tests with enable-checking
Results for latest snapshot egcs-20000828 on mips-sgi-irix6.5 show a few
failures, but results are comparable to SGI f77 7.2.1 so they may be
IEEE-754 vs CPU issues
With compiler flags "-O3 -funroll-loops" and default -n32 ABI I get
csep.out: CST: 1 out of 4662 tests failed to pass the threshold
csep.out: CST drivers: 1 out of 11664 tests failed to pass the
threshold
csep.out: CST: 1 out of 4662 tests failed to pass the threshold
csep.out: CST drivers: 1 out of 11664 tests failed to pass the
threshold
csep.out: CST drivers: 1 out of 11664 tests failed to pass the
threshold
ctest.out: CPB: 12 out of 3458 tests failed to pass the threshold
ctest.out: CPB drivers: 3 out of 4750 tests failed to pass the
threshold
ctest.out: CLS drivers: 1 out of 65268 tests failed to pass the
threshold
dgd.out: DXV drivers: 200 out of 5000 tests failed to pass the
threshold
sgd.out: SXV drivers: 37 out of 5000 tests failed to pass the
threshold
ssep.out: SST: 1 out of 4662 tests failed to pass the threshold
ssep.out: SST drivers: 1 out of 14256 tests failed to pass the
threshold
stest.out: SPB drivers: 1 out of 4750 tests failed to pass the
threshold
stest.out: SLS drivers: 3 out of 65268 tests failed to pass the
threshold
zgd.out: ZXV drivers: 24 out of 5000 tests failed to pass the
threshold
The failures of the _XV drivers in _gd.out are known LAPACK problems.
{C,S}ST failures in {c,s}sep.out are for matrix type = 9. Other systems
have
occasional problems with these tests, too.
Failures in {c,s}test.out are near misses. Similar failures are seem with
MIPSPro f77 v7.2.1.
With "-mabi=64 -O3 -funroll-loops" the results are similar, although there
is an ICE on one file - it is in GNATS.
::::::::::::::: Gory details of failures :::::::::::::::::
::::::::::::
csep.out
::::::::::::
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
60.73
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)
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=2247,3966,3935,3605, result 47 is
9.034E+05
CST drivers: 1 out of 11664 tests failed to pass the threshold
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 36 is
1460.44
CST: 1 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 47 is
334.75
CST drivers: 1 out of 11664 tests failed to pass the threshold
SEP: NB = 10, NBMIN = 2, NX = 1
All tests for CST passed the threshold ( 4662 tests run)
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=1139,3868,2100,1985, result 101 is
8.458E+05
CST drivers: 1 out of 11664 tests failed to pass the threshold
:::::::::::
ctest.out
:::::::::::
CPB: Hermitian positive definite band matrices
Matrix types:
1. Random, CNDNUM = 2 5. Random, CNDNUM = sqrt(0.1/EPS)
*2. First row and column zero 6. Random, CNDNUM = 0.1/EPS
*3. Last row and column zero 7. Scaled near underflow
*4. Middle row and column zero 8. Scaled near overflow
(* - tests error exits from CPBTRF, no test ratios are computed)
Test ratios:
1: norm( U' * U - A ) / ( N * norm(A) * EPS ), or
norm( L * L' - A ) / ( N * norm(A) * EPS )
2: norm( B - A * X ) / ( norm(A) * norm(X) * EPS )
3: norm( X - XACT ) / ( norm(XACT) * CNDNUM * EPS )
4: norm( X - XACT ) / ( norm(XACT) * CNDNUM * EPS ), refined
5: norm( X - XACT ) / ( norm(XACT) * (error bound) )
6: (backward error) / EPS
7: RCOND * CNDNUM - 1.0
Messages:
UPLO='U', N= 50, KD= 37, NRHS= 1, type 7, test( 2) = 36.192
UPLO='U', N= 50, KD= 37, NRHS= 1, type 7, test( 3) = 34.406
UPLO='U', N= 50, KD= 37, NRHS= 2, type 7, test( 2) = 42.724
UPLO='U', N= 50, KD= 37, NRHS= 2, type 7, test( 3) = 38.736
UPLO='U', N= 50, KD= 37, NRHS= 15, type 7, test( 2) = 44.307
UPLO='U', N= 50, KD= 37, NRHS= 15, type 7, test( 3) = 44.029
UPLO='L', N= 50, KD= 37, NRHS= 1, type 7, test( 2) = 45.681
UPLO='L', N= 50, KD= 37, NRHS= 1, type 7, test( 3) = 38.842
UPLO='L', N= 50, KD= 37, NRHS= 2, type 7, test( 2) = 48.998
UPLO='L', N= 50, KD= 37, NRHS= 2, type 7, test( 3) = 41.583
UPLO='L', N= 50, KD= 37, NRHS= 15, type 7, test( 2) = 50.468
UPLO='L', N= 50, KD= 37, NRHS= 15, type 7, test( 3) = 41.896
...
CPB drivers: Hermitian positive definite band matrices
Matrix types:
1. Random, CNDNUM = 2 5. Random, CNDNUM = sqrt(0.1/EPS)
*2. First row and column zero 6. Random, CNDNUM = 0.1/EPS
*3. Last row and column zero 7. Scaled near underflow
*4. Middle row and column zero 8. Scaled near overflow
(* - tests error exits from CPBTRF, no test ratios are computed)
Test ratios:
1: norm( U' * U - A ) / ( N * norm(A) * EPS ), or
norm( L * L' - A ) / ( N * norm(A) * EPS )
2: norm( B - A * X ) / ( norm(A) * norm(X) * EPS )
3: norm( X - XACT ) / ( norm(XACT) * CNDNUM * EPS )
4: norm( X - XACT ) / ( norm(XACT) * (error bound) )
5: (backward error) / EPS
6: RCOND * CNDNUM - 1.0
Messages:
CPBSV , UPLO='U', N = 50, KD = 37, type 7, test(2)= 45.423
CPBSV , UPLO='U', N = 50, KD = 37, type 7, test(3)= 38.374
CPBSV , UPLO='L', N = 50, KD = 37, type 7, test(3)= 36.370
CPB drivers: 3 out of 4750 tests failed to pass the threshold
...
CLS: Least squares driver routines
Matrix types (1-3: full rank, 4-6: rank deficient):
1 and 4. Normal scaling
2 and 5. Scaled near overflow
3 and 6. Scaled near underflow
Test ratios:
(1-2: CGELS, 3-6: CGELSX, 7-10: CGELSY, 11-14: CGELSS, 15-18: CGELSD)
1: norm( B - A * X ) / ( max(M,N) * norm(A) * norm(X) * EPS )
2: norm( (A*X-B)' *A ) / ( max(M,N,NRHS) * norm(A) * norm(B) * EPS )
if TRANS='N' and M.GE.N or TRANS='T' and M.LT.N, otherwise
check if X is in the row space of A or A' (overdetermined case)
3: norm(svd(A)-svd(R)) / ( min(M,N) * norm(svd(R)) * EPS )
4: norm( B - A * X ) / ( max(M,N) * norm(A) * norm(X) * EPS )
5: norm( (A*X-B)' *A ) / ( max(M,N,NRHS) * norm(A) * norm(B) * EPS )
6: Check if X is in the row space of A or A'
7-10: same as 3-6 11-14: same as 3-6 15-18: same as 3-6
Messages:
TRANS='C', M= 50, N= 1, NRHS= 15, NB= 1, type 3, test( 1)= 40.213
CLS drivers: 1 out of 65268 tests failed to pass the threshold
:::::::::::
ssep.out
:::::::::::
SST -- Real Symmetric eigenvalue problem
Matrix types (see SCHKST 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 Symmetric 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 SCHKST for details.
N= 20, seed=2989,1119,3793,1781, type 9, test(36)= 195.
SST: 1 out of 4662 tests failed to pass the threshold
...
SST -- Real Symmetric 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 Symmetric 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 sdrvst.f
Matrix order= 20, type= 9, seed=3966,3411,3597,2265, result 72 is
61.73
SST drivers: 1 out of 14256 tests failed to pass the threshold
:::::::::::
stest.out
:::::::::::
SPB drivers: Symmetric positive definite band matrices
Matrix types:
1. Random, CNDNUM = 2 5. Random, CNDNUM = sqrt(0.1/EPS)
*2. First row and column zero 6. Random, CNDNUM = 0.1/EPS
*3. Last row and column zero 7. Scaled near underflow
*4. Middle row and column zero 8. Scaled near overflow
(* - tests error exits from SPBTRF, no test ratios are computed)
Test ratios:
1: norm( U' * U - A ) / ( N * norm(A) * EPS ), or
norm( L * L' - A ) / ( N * norm(A) * EPS )
2: norm( B - A * X ) / ( norm(A) * norm(X) * EPS )
3: norm( X - XACT ) / ( norm(XACT) * CNDNUM * EPS )
4: norm( X - XACT ) / ( norm(XACT) * (error bound) )
5: (backward error) / EPS
6: RCOND * CNDNUM - 1.0
Messages:
SPBSV , UPLO='L', N = 50, KD = 37, type 7, test(3)= 30.333
SPB drivers: 1 out of 4750 tests failed to pass the threshold
...
SLS: Least squares driver routines
Matrix types (1-3: full rank, 4-6: rank deficient):
1 and 4. Normal scaling
2 and 5. Scaled near overflow
3 and 6. Scaled near underflow
Test ratios:
(1-2: SGELS, 3-6: SGELSX, 7-10: SGELSY, 11-14: SGELSS, 15-18: SGELSD)
1: norm( B - A * X ) / ( max(M,N) * norm(A) * norm(X) * EPS )
2: norm( (A*X-B)' *A ) / ( max(M,N,NRHS) * norm(A) * norm(B) * EPS )
if TRANS='N' and M.GE.N or TRANS='T' and M.LT.N, otherwise
check if X is in the row space of A or A' (overdetermined case)
3: norm(svd(A)-svd(R)) / ( min(M,N) * norm(svd(R)) * EPS )
4: norm( B - A * X ) / ( max(M,N) * norm(A) * norm(X) * EPS )
5: norm( (A*X-B)' *A ) / ( max(M,N,NRHS) * norm(A) * norm(B) * EPS )
6: Check if X is in the row space of A or A'
7-10: same as 3-6 11-14: same as 3-6 15-18: same as 3-6
Messages:
TRANS='T', M= 50, N= 1, NRHS= 2, NB= 3, type 3, test( 1)= 139.96
TRANS='T', M= 50, N= 1, NRHS= 15, NB= 3, type 3, test( 1)= 1235.6
TRANS='T', M= 50, N= 1, NRHS= 15, NB= 3, type 3, test( 1)= 91.860
SLS drivers: 3 out of 65268 tests failed to pass the threshold
>