Tests of the Generalized Nonsymmetric Eigenvalue Problem Driver DGGES LAPACK VERSION 3.0, released June 30, 1999 The following parameter values will be used: M: 2 6 10 12 20 N: 2 6 10 12 20 NB: 1 NBMIN: 1 NX: 1 NS: 2 MAXB: 1 Relative machine underflow is taken to be 0.733442-2465 Relative machine overflow is taken to be 0.136344+2466 Relative machine precision is taken to be 0.142109E-13 Routines pass computational tests if test ratio is less than 10.00 All tests for DGS drivers passed the threshold ( 1560 tests run) ----------------------------------------------------------------------- Tests of the Generalized Nonsymmetric Eigenvalue Problem Driver DGGEV LAPACK VERSION 3.0, released June 30, 1999 The following parameter values will be used: M: 2 6 8 10 15 20 N: 2 6 8 10 15 20 NB: 1 NBMIN: 1 NX: 1 NS: 2 MAXB: 1 Relative machine underflow is taken to be 0.733442-2465 Relative machine overflow is taken to be 0.136344+2466 Relative machine precision is taken to be 0.142109E-13 Routines pass computational tests if test ratio is less than 10.00 DGV -- Real Generalized eigenvalue problem driver Matrix types (see DDRGEV for details): Special Matrices: (J'=transposed Jordan block) 1=(0,0) 2=(I,0) 3=(0,I) 4=(I,I) 5=(J',J') 6=(diag(J',I), diag(I,J')) Diagonal Matrices: ( D=diag(0,1,2,...) ) 7=(D,I) 9=(large*D, small*I) 11=(large*I, small*D) 13=(large*D, large*I) 8=(I,D) 10=(small*D, large*I) 12=(small*I, large*D) 14=(small*D, small*I) 15=(D, reversed D) Matrices Rotated by Random Orthogonal Matrices U, V: 16=Transposed Jordan Blocks 19=geometric alpha, beta=0,1 17=arithm. alpha&beta 20=arithmetic alpha, beta=0,1 18=clustered alpha, beta=0,1 21=random alpha, beta=0,1 Large & Small Matrices: 22=(large, small) 23=(small,large) 24=(small,small) 25=(large,large) 26=random O(1) matrices. Tests performed: 1 = max | ( b A - a B )'*l | / const., 2 = | |VR(i)| - 1 | / ulp, 3 = max | ( b A - a B )*r | / const. 4 = | |VL(i)| - 1 | / ulp, 5 = 0 if W same no matter if r or l computed, 6 = 0 if l same no matter if l computed, 7 = 0 if r same no matter if r computed, Matrix order= 6, type=12, seed=1147,1097,1906, 485, result 5 is 3.518E+13 DGV drivers: 1 out of 1092 tests failed to pass the threshold ----------------------------------------------------------------------- Tests of the Generalized Nonsymmetric Eigenvalue Problem Expert Driver DGGESX LAPACK VERSION 3.0, released June 30, 1999 The following parameter values will be used: N: 6 NB: 1 NBMIN: 1 NX: 1 NS: 2 MAXB: 1 Relative machine underflow is taken to be 0.733442-2465 Relative machine overflow is taken to be 0.136344+2466 Relative machine precision is taken to be 0.142109E-13 Routines pass computational tests if test ratio is less than 10.00 DGX -- Real Expert Generalized Schur form problem driver Matrix types: 1: A is a block diagonal matrix of Jordan blocks and B is the identity matrix, 2: A and B are upper triangular matrices, 3: A and B are as type 2, but each second diagonal block in A_11 and each third diaongal block in A_22 are 2x2 blocks, 4: A and B are block diagonal matrices, 5: (A,B) has potentially close or common eigenvalues. Tests performed: (S is Schur, T is triangular, Q and Z are orthogonal, a is alpha, b is beta, and ' means transpose.) 1 = | A - Q S Z' | / ( |A| n ulp ) 2 = | B - Q T Z' | / ( |B| n ulp ) 3 = | I - QQ' | / ( n ulp ) 4 = | I - ZZ' | / ( n ulp ) 5 = 1/ULP if A is not in Schur form S 6 = difference between (alpha,beta) and diagonals of (S,T) 7 = 1/ULP if SDIM is not the correct number of selected eigenvalues 8 = 1/ULP if DIFEST/DIFTRU > 10*THRESH or DIFTRU/DIFEST > 10*THRESH 9 = 1/ULP if DIFEST <> 0 or DIFTRU > ULP*norm(A,B) when reordering fails 10 = 1/ULP if PLEST/PLTRU > THRESH or PLTRU/PLEST > THRESH ( Test 10 is only for input examples ) Matrix order= 6, type= 1, a=0.1686E-06, order(A_11)= 3, result 8 is 0.2585E+08 Matrix order= 6, type= 1, a=0.1686E-06, order(A_11)= 4, result 8 is 0.1577E+09 Matrix order= 6, type= 1, a=0.1686E-06, order(A_11)= 1, result 8 is 0.7372E+06 Matrix order= 6, type= 1, a=0.1686E-06, order(A_11)= 2, result 8 is 0.1811E+08 Matrix order= 5, type= 1, a=0.1686E-06, order(A_11)= 3, result 8 is 154.62 Matrix order= 6, type= 1, a=0.1686E-06, order(A_11)= 5, result 8 is 0.7193E+07 Matrix order= 3, type= 5, a=0.5932E+07, order(A_11)= 1, result 8 is 0.1437E+06 Matrix order= 5, type= 5, a=0.5932E+07, order(A_11)= 1, result 8 is 0.1936E+06 Matrix order= 3, type= 5, a=0.5932E+07, order(A_11)= 2, result 8 is 0.1500E+06 DGX drivers: 9 out of 2250 tests failed to pass the threshold ----------------------------------------------------------------------- Tests of the Generalized Nonsymmetric Eigenvalue Problem Expert Driver DGGESX LAPACK VERSION 3.0, released June 30, 1999 The following parameter values will be used: N: 0 NB: 1 NBMIN: 1 NX: 1 NS: 2 MAXB: 1 Relative machine underflow is taken to be 0.733442-2465 Relative machine overflow is taken to be 0.136344+2466 Relative machine precision is taken to be 0.142109E-13 Routines pass computational tests if test ratio is less than 10.00 All tests for DGX drivers passed the threshold ( 20 tests run) ----------------------------------------------------------------------- Tests of the Generalized Nonsymmetric Eigenvalue Problem Expert Driver DGGEVX LAPACK VERSION 3.0, released June 30, 1999 The following parameter values will be used: N: 5 NB: 1 NBMIN: 1 NX: 1 NS: 2 MAXB: 1 Relative machine underflow is taken to be 0.733442-2465 Relative machine overflow is taken to be 0.136344+2466 Relative machine precision is taken to be 0.142109E-13 Routines pass computational tests if test ratio is less than 10.00 DXV -- Real Expert Eigenvalue/vector problem driver Matrix types: TYPE 1: Da is diagonal, Db is identity, A = Y^(-H) Da X^(-1), B = Y^(-H) Db X^(-1) YH and X are left and right eigenvectors. TYPE 2: Da is quasi-diagonal, Db is identity, A = Y^(-H) Da X^(-1), B = Y^(-H) Db X^(-1) YH and X are left and right eigenvectors. Tests performed: a is alpha, b is beta, l is a left eigenvector, r is a right eigenvector and ' means transpose. 1 = max | ( b A - a B )' l | / const. 2 = max | ( b A - a B ) r | / const. 3 = max ( Sest/Stru, Stru/Sest ) over all eigenvalues 4 = max( DIFest/DIFtru, DIFtru/DIFest ) over the 1st and 5th eigenvectors Type= 2, IWA= 1, IWB= 1, IWX= 5, IWY= 1, result 4 is 1119.91 Type= 2, IWA= 1, IWB= 1, IWX= 5, IWY= 2, result 4 is 1168.72 Type= 2, IWA= 1, IWB= 1, IWX= 5, IWY= 3, result 4 is 1570.81 Type= 2, IWA= 1, IWB= 1, IWX= 5, IWY= 4, result 4 is 1501.66 Type= 2, IWA= 1, IWB= 1, IWX= 5, IWY= 5, result 4 is 490.19 Type= 2, IWA= 1, IWB= 2, IWX= 5, IWY= 5, result 4 is 673.31 Type= 2, IWA= 1, IWB= 3, IWX= 5, IWY= 5, result 4 is 394.38 Type= 2, IWA= 1, IWB= 4, IWX= 5, IWY= 5, result 4 is 141.48 Type= 2, IWA= 1, IWB= 5, IWX= 1, IWY= 1, result 4 is 1218.04 Type= 2, IWA= 1, IWB= 5, IWX= 1, IWY= 2, result 4 is 1135.61 Type= 2, IWA= 1, IWB= 5, IWX= 1, IWY= 3, result 4 is 630.87 Type= 2, IWA= 1, IWB= 5, IWX= 2, IWY= 1, result 4 is 1135.61 Type= 2, IWA= 1, IWB= 5, IWX= 2, IWY= 2, result 4 is 1102.95 Type= 2, IWA= 1, IWB= 5, IWX= 2, IWY= 3, result 4 is 629.06 Type= 2, IWA= 1, IWB= 5, IWX= 3, IWY= 1, result 4 is 630.87 Type= 2, IWA= 1, IWB= 5, IWX= 3, IWY= 2, result 4 is 629.06 Type= 2, IWA= 1, IWB= 5, IWX= 3, IWY= 3, result 4 is 504.82 Type= 2, IWA= 2, IWB= 1, IWX= 5, IWY= 5, result 4 is 1350.15 Type= 2, IWA= 2, IWB= 2, IWX= 5, IWY= 5, result 4 is 1083.76 Type= 2, IWA= 2, IWB= 3, IWX= 5, IWY= 5, result 4 is 485.53 Type= 2, IWA= 2, IWB= 4, IWX= 5, IWY= 5, result 4 is 143.11 Type= 2, IWA= 2, IWB= 5, IWX= 1, IWY= 1, result 4 is 1218.09 Type= 2, IWA= 2, IWB= 5, IWX= 1, IWY= 2, result 4 is 1135.66 Type= 2, IWA= 2, IWB= 5, IWX= 1, IWY= 3, result 4 is 630.90 Type= 2, IWA= 2, IWB= 5, IWX= 2, IWY= 1, result 4 is 1135.66 Type= 2, IWA= 2, IWB= 5, IWX= 2, IWY= 2, result 4 is 1103.00 Type= 2, IWA= 2, IWB= 5, IWX= 2, IWY= 3, result 4 is 629.08 Type= 2, IWA= 2, IWB= 5, IWX= 3, IWY= 1, result 4 is 630.90 Type= 2, IWA= 2, IWB= 5, IWX= 3, IWY= 2, result 4 is 629.08 Type= 2, IWA= 2, IWB= 5, IWX= 3, IWY= 3, result 4 is 504.84 Type= 2, IWA= 3, IWB= 1, IWX= 5, IWY= 5, result 4 is 1421.21 Type= 2, IWA= 3, IWB= 2, IWX= 5, IWY= 5, result 4 is 1406.28 Type= 2, IWA= 3, IWB= 3, IWX= 5, IWY= 5, result 4 is 944.50 Type= 2, IWA= 3, IWB= 4, IWX= 5, IWY= 5, result 4 is 156.07 Type= 2, IWA= 3, IWB= 5, IWX= 1, IWY= 1, result 4 is 1218.54 Type= 2, IWA= 3, IWB= 5, IWX= 1, IWY= 2, result 4 is 1136.08 Type= 2, IWA= 3, IWB= 5, IWX= 1, IWY= 3, result 4 is 631.13 Type= 2, IWA= 3, IWB= 5, IWX= 2, IWY= 1, result 4 is 1136.08 Type= 2, IWA= 3, IWB= 5, IWX= 2, IWY= 2, result 4 is 1103.40 Type= 2, IWA= 3, IWB= 5, IWX= 2, IWY= 3, result 4 is 629.32 Type= 2, IWA= 3, IWB= 5, IWX= 3, IWY= 1, result 4 is 631.13 Type= 2, IWA= 3, IWB= 5, IWX= 3, IWY= 2, result 4 is 629.32 Type= 2, IWA= 3, IWB= 5, IWX= 3, IWY= 3, result 4 is 505.03 Type= 2, IWA= 4, IWB= 1, IWX= 5, IWY= 5, result 4 is 189.92 Type= 2, IWA= 4, IWB= 2, IWX= 5, IWY= 5, result 4 is 191.27 Type= 2, IWA= 4, IWB= 3, IWX= 5, IWY= 5, result 4 is 201.65 Type= 2, IWA= 4, IWB= 4, IWX= 5, IWY= 5, result 4 is 162.77 Type= 2, IWA= 4, IWB= 5, IWX= 1, IWY= 1, result 4 is 1223.05 Type= 2, IWA= 4, IWB= 5, IWX= 1, IWY= 2, result 4 is 1140.28 Type= 2, IWA= 4, IWB= 5, IWX= 1, IWY= 3, result 4 is 633.46 Type= 2, IWA= 4, IWB= 5, IWX= 2, IWY= 1, result 4 is 1140.28 Type= 2, IWA= 4, IWB= 5, IWX= 2, IWY= 2, result 4 is 1107.49 Type= 2, IWA= 4, IWB= 5, IWX= 2, IWY= 3, result 4 is 631.64 Type= 2, IWA= 4, IWB= 5, IWX= 3, IWY= 1, result 4 is 633.46 Type= 2, IWA= 4, IWB= 5, IWX= 3, IWY= 2, result 4 is 631.64 Type= 2, IWA= 4, IWB= 5, IWX= 3, IWY= 3, result 4 is 506.90 Type= 2, IWA= 5, IWB= 1, IWX= 1, IWY= 1, result 4 is 2.966E+06 Type= 2, IWA= 5, IWB= 1, IWX= 1, IWY= 2, result 4 is 2.765E+06 Type= 2, IWA= 5, IWB= 1, IWX= 1, IWY= 3, result 4 is 1.536E+06 Type= 2, IWA= 5, IWB= 1, IWX= 1, IWY= 4, result 4 is 2.088E+05 Type= 2, IWA= 5, IWB= 1, IWX= 1, IWY= 5, result 4 is 861.43 Type= 2, IWA= 5, IWB= 1, IWX= 2, IWY= 1, result 4 is 2.765E+06 Type= 2, IWA= 5, IWB= 1, IWX= 2, IWY= 2, result 4 is 2.685E+06 Type= 2, IWA= 5, IWB= 1, IWX= 2, IWY= 3, result 4 is 1.531E+06 Type= 2, IWA= 5, IWB= 1, IWX= 2, IWY= 4, result 4 is 2.088E+05 Type= 2, IWA= 5, IWB= 1, IWX= 2, IWY= 5, result 4 is 861.43 Type= 2, IWA= 5, IWB= 1, IWX= 3, IWY= 1, result 4 is 1.536E+06 Type= 2, IWA= 5, IWB= 1, IWX= 3, IWY= 2, result 4 is 1.531E+06 Type= 2, IWA= 5, IWB= 1, IWX= 3, IWY= 3, result 4 is 1.229E+06 Type= 2, IWA= 5, IWB= 1, IWX= 3, IWY= 4, result 4 is 2.077E+05 Type= 2, IWA= 5, IWB= 1, IWX= 3, IWY= 5, result 4 is 861.43 Type= 2, IWA= 5, IWB= 1, IWX= 4, IWY= 1, result 4 is 2.088E+05 Type= 2, IWA= 5, IWB= 1, IWX= 4, IWY= 2, result 4 is 2.088E+05 Type= 2, IWA= 5, IWB= 1, IWX= 4, IWY= 3, result 4 is 2.077E+05 Type= 2, IWA= 5, IWB= 1, IWX= 4, IWY= 4, result 4 is 1.480E+05 Type= 2, IWA= 5, IWB= 1, IWX= 4, IWY= 5, result 4 is 861.43 Type= 2, IWA= 5, IWB= 1, IWX= 5, IWY= 1, result 4 is 861.43 Type= 2, IWA= 5, IWB= 1, IWX= 5, IWY= 2, result 4 is 861.43 Type= 2, IWA= 5, IWB= 1, IWX= 5, IWY= 3, result 4 is 861.43 Type= 2, IWA= 5, IWB= 1, IWX= 5, IWY= 4, result 4 is 861.43 Type= 2, IWA= 5, IWB= 1, IWX= 5, IWY= 5, result 4 is 609.25 Type= 2, IWA= 5, IWB= 2, IWX= 1, IWY= 1, result 4 is 2.697E+06 Type= 2, IWA= 5, IWB= 2, IWX= 1, IWY= 2, result 4 is 2.514E+06 Type= 2, IWA= 5, IWB= 2, IWX= 1, IWY= 3, result 4 is 1.397E+06 Type= 2, IWA= 5, IWB= 2, IWX= 1, IWY= 4, result 4 is 1.899E+05 Type= 2, IWA= 5, IWB= 2, IWX= 1, IWY= 5, result 4 is 783.47 Type= 2, IWA= 5, IWB= 2, IWX= 2, IWY= 1, result 4 is 2.514E+06 Type= 2, IWA= 5, IWB= 2, IWX= 2, IWY= 2, result 4 is 2.442E+06 Type= 2, IWA= 5, IWB= 2, IWX= 2, IWY= 3, result 4 is 1.393E+06 Type= 2, IWA= 5, IWB= 2, IWX= 2, IWY= 4, result 4 is 1.899E+05 Type= 2, IWA= 5, IWB= 2, IWX= 2, IWY= 5, result 4 is 783.47 Type= 2, IWA= 5, IWB= 2, IWX= 3, IWY= 1, result 4 is 1.397E+06 Type= 2, IWA= 5, IWB= 2, IWX= 3, IWY= 2, result 4 is 1.393E+06 Type= 2, IWA= 5, IWB= 2, IWX= 3, IWY= 3, result 4 is 1.118E+06 Type= 2, IWA= 5, IWB= 2, IWX= 3, IWY= 4, result 4 is 1.889E+05 Type= 2, IWA= 5, IWB= 2, IWX= 3, IWY= 5, result 4 is 783.47 Type= 2, IWA= 5, IWB= 2, IWX= 4, IWY= 1, result 4 is 1.899E+05 Type= 2, IWA= 5, IWB= 2, IWX= 4, IWY= 2, result 4 is 1.899E+05 Type= 2, IWA= 5, IWB= 2, IWX= 4, IWY= 3, result 4 is 1.889E+05 Type= 2, IWA= 5, IWB= 2, IWX= 4, IWY= 4, result 4 is 1.346E+05 Type= 2, IWA= 5, IWB= 2, IWX= 4, IWY= 5, result 4 is 783.47 Type= 2, IWA= 5, IWB= 2, IWX= 5, IWY= 1, result 4 is 783.47 Type= 2, IWA= 5, IWB= 2, IWX= 5, IWY= 2, result 4 is 783.47 Type= 2, IWA= 5, IWB= 2, IWX= 5, IWY= 3, result 4 is 783.47 Type= 2, IWA= 5, IWB= 2, IWX= 5, IWY= 4, result 4 is 783.47 Type= 2, IWA= 5, IWB= 2, IWX= 5, IWY= 5, result 4 is 554.11 Type= 2, IWA= 5, IWB= 3, IWX= 1, IWY= 1, result 4 is 1.484E+06 Type= 2, IWA= 5, IWB= 3, IWX= 1, IWY= 2, result 4 is 1.383E+06 Type= 2, IWA= 5, IWB= 3, IWX= 1, IWY= 3, result 4 is 7.686E+05 Type= 2, IWA= 5, IWB= 3, IWX= 1, IWY= 4, result 4 is 1.045E+05 Type= 2, IWA= 5, IWB= 3, IWX= 1, IWY= 5, result 4 is 431.07 Type= 2, IWA= 5, IWB= 3, IWX= 2, IWY= 1, result 4 is 1.383E+06 Type= 2, IWA= 5, IWB= 3, IWX= 2, IWY= 2, result 4 is 1.344E+06 Type= 2, IWA= 5, IWB= 3, IWX= 2, IWY= 3, result 4 is 7.664E+05 Type= 2, IWA= 5, IWB= 3, IWX= 2, IWY= 4, result 4 is 1.045E+05 Type= 2, IWA= 5, IWB= 3, IWX= 2, IWY= 5, result 4 is 431.07 Type= 2, IWA= 5, IWB= 3, IWX= 3, IWY= 1, result 4 is 7.686E+05 Type= 2, IWA= 5, IWB= 3, IWX= 3, IWY= 2, result 4 is 7.664E+05 Type= 2, IWA= 5, IWB= 3, IWX= 3, IWY= 3, result 4 is 6.151E+05 Type= 2, IWA= 5, IWB= 3, IWX= 3, IWY= 4, result 4 is 1.040E+05 Type= 2, IWA= 5, IWB= 3, IWX= 3, IWY= 5, result 4 is 431.07 Type= 2, IWA= 5, IWB= 3, IWX= 4, IWY= 1, result 4 is 1.045E+05 Type= 2, IWA= 5, IWB= 3, IWX= 4, IWY= 2, result 4 is 1.045E+05 Type= 2, IWA= 5, IWB= 3, IWX= 4, IWY= 3, result 4 is 1.040E+05 Type= 2, IWA= 5, IWB= 3, IWX= 4, IWY= 4, result 4 is 7.407E+04 Type= 2, IWA= 5, IWB= 3, IWX= 4, IWY= 5, result 4 is 431.07 Type= 2, IWA= 5, IWB= 3, IWX= 5, IWY= 1, result 4 is 431.07 Type= 2, IWA= 5, IWB= 3, IWX= 5, IWY= 2, result 4 is 431.07 Type= 2, IWA= 5, IWB= 3, IWX= 5, IWY= 3, result 4 is 431.07 Type= 2, IWA= 5, IWB= 3, IWX= 5, IWY= 4, result 4 is 431.07 Type= 2, IWA= 5, IWB= 3, IWX= 5, IWY= 5, result 4 is 304.87 Type= 2, IWA= 5, IWB= 4, IWX= 1, IWY= 1, result 4 is 2.708E+05 Type= 2, IWA= 5, IWB= 4, IWX= 1, IWY= 2, result 4 is 2.525E+05 Type= 2, IWA= 5, IWB= 4, IWX= 1, IWY= 3, result 4 is 1.403E+05 Type= 2, IWA= 5, IWB= 4, IWX= 1, IWY= 4, result 4 is 1.906E+04 Type= 2, IWA= 5, IWB= 4, IWX= 2, IWY= 1, result 4 is 2.525E+05 Type= 2, IWA= 5, IWB= 4, IWX= 2, IWY= 2, result 4 is 2.452E+05 Type= 2, IWA= 5, IWB= 4, IWX= 2, IWY= 3, result 4 is 1.399E+05 Type= 2, IWA= 5, IWB= 4, IWX= 2, IWY= 4, result 4 is 1.906E+04 Type= 2, IWA= 5, IWB= 4, IWX= 3, IWY= 1, result 4 is 1.403E+05 Type= 2, IWA= 5, IWB= 4, IWX= 3, IWY= 2, result 4 is 1.399E+05 Type= 2, IWA= 5, IWB= 4, IWX= 3, IWY= 3, result 4 is 1.122E+05 Type= 2, IWA= 5, IWB= 4, IWX= 3, IWY= 4, result 4 is 1.897E+04 Type= 2, IWA= 5, IWB= 4, IWX= 4, IWY= 1, result 4 is 1.906E+04 Type= 2, IWA= 5, IWB= 4, IWX= 4, IWY= 2, result 4 is 1.906E+04 Type= 2, IWA= 5, IWB= 4, IWX= 4, IWY= 3, result 4 is 1.897E+04 Type= 2, IWA= 5, IWB= 4, IWX= 4, IWY= 4, result 4 is 1.352E+04 Type= 2, IWA= 5, IWB= 5, IWX= 1, IWY= 1, result 4 is 3443.44 Type= 2, IWA= 5, IWB= 5, IWX= 1, IWY= 2, result 4 is 3210.02 Type= 2, IWA= 5, IWB= 5, IWX= 1, IWY= 3, result 4 is 1783.27 Type= 2, IWA= 5, IWB= 5, IWX= 1, IWY= 4, result 4 is 242.39 Type= 2, IWA= 5, IWB= 5, IWX= 2, IWY= 1, result 4 is 3210.02 Type= 2, IWA= 5, IWB= 5, IWX= 2, IWY= 2, result 4 is 3117.73 Type= 2, IWA= 5, IWB= 5, IWX= 2, IWY= 3, result 4 is 1778.17 Type= 2, IWA= 5, IWB= 5, IWX= 2, IWY= 4, result 4 is 242.38 Type= 2, IWA= 5, IWB= 5, IWX= 3, IWY= 1, result 4 is 1783.27 Type= 2, IWA= 5, IWB= 5, IWX= 3, IWY= 2, result 4 is 1778.17 Type= 2, IWA= 5, IWB= 5, IWX= 3, IWY= 3, result 4 is 1427.08 Type= 2, IWA= 5, IWB= 5, IWX= 3, IWY= 4, result 4 is 241.21 Type= 2, IWA= 5, IWB= 5, IWX= 4, IWY= 1, result 4 is 242.39 Type= 2, IWA= 5, IWB= 5, IWX= 4, IWY= 2, result 4 is 242.38 Type= 2, IWA= 5, IWB= 5, IWX= 4, IWY= 3, result 4 is 241.21 Type= 2, IWA= 5, IWB= 5, IWX= 4, IWY= 4, result 4 is 171.84 DXV drivers: 163 out of 5000 tests failed to pass the threshold ----------------------------------------------------------------------- Tests of the Generalized Nonsymmetric Eigenvalue Problem Expert Driver DGGEVX LAPACK VERSION 3.0, released June 30, 1999 The following parameter values will be used: N: 0 NB: 1 NBMIN: 1 NX: 1 NS: 2 MAXB: 1 Relative machine underflow is taken to be 0.733442-2465 Relative machine overflow is taken to be 0.136344+2466 Relative machine precision is taken to be 0.142109E-13 Routines pass computational tests if test ratio is less than 10.00 All tests for DXV drivers passed the threshold ( 8 tests run) ----------------------------------------------------------------------- End of tests Total time used = 9.74 seconds