TESTING/EIG/schkee.F(3)
π NAME π
TESTING/EIG/schkee.F β bestie, this is a page straight out of the LAPACK manual, and it's all about schkee, the test program for the real (single precision) eigenvalue problem routines. Periodt.
β¨ SYNOPSIS β¨
Functions/Subroutines
program schkeeβ the whole vibe, right hereSCHKEE
πΈ Function/Subroutine Documentation πΈ
program schkee π«Ά
SCHKEE
π Purpose π
SCHKEE tests the REAL LAPACK subroutines for the matrix eigenvalue problem. The test paths in this version are, and I need you to pay attention because there's a LOT of them:
- NEP (Nonsymmetric Eigenvalue Problem): Test
SGEHRD,SORGHR,SHSEQR,STREVC,SHSEIN, andSORMHR - SEP (Symmetric Eigenvalue Problem): Test
SSYTRD,SORGTR,SSTEQR,SSTERF,SSTEIN,SSTEDC, and driversSSYEV(X),SSBEV(X),SSPEV(X),SSTEV(X),SSYEVD,SSBEVD,SSPEVD,SSTEVD - SVD (Singular Value Decomposition): Test
SGEBRD,SORGBR,SBDSQR,SBDSDCand the driversSGESVD,SGESDD - SEV (Nonsymmetric Eigenvalue/eigenvector Driver): Test
SGEEV - SES (Nonsymmetric Schur form Driver): Test
SGEES - SVX (Nonsymmetric Eigenvalue/eigenvector Expert Driver): Test
SGEEVX - SSX (Nonsymmetric Schur form Expert Driver): Test
SGEESX - SGG (Generalized Nonsymmetric Eigenvalue Problem): Test
SGGHD3,SGGBAL,SGGBAK,SHGEQZ, andSTGEVC - SGS (Generalized Nonsymmetric Schur form Driver): Test
SGGES - SGV (Generalized Nonsymmetric Eigenvalue/eigenvector Driver): Test
SGGEV - SGX (Generalized Nonsymmetric Schur form Expert Driver): Test
SGGESX - SXV (Generalized Nonsymmetric Eigenvalue/eigenvector Expert Driver): Test
SGGEVX - SSG (Symmetric Generalized Eigenvalue Problem): Test
SSYGST,SSYGV,SSYGVD,SSYGVX,SSPGST,SSPGV,SSPGVD,SSPGVX,SSBGST,SSBGV,SSBGVD, andSSBGVX - SSB (Symmetric Band Eigenvalue Problem): Test
SSBTRD - SBB (Band Singular Value Decomposition): Test
SGBBRD - SEC (Eigencondition estimation): Test
SLALN2,SLASY2,SLAEQU,SLAEXC,STRSYL,STREXC,STRSNA,STRSEN, andSLAQTR - SBL (Balancing a general matrix) Test
SGEBAL - SBK (Back transformation on a balanced matrix) Test
SGEBAK - SGL (Balancing a matrix pair) Test
SGGBAL - SGK (Back transformation on a matrix pair) Test
SGGBAK - GLM (Generalized Linear Regression Model): Tests
SGGGLM - GQR (Generalized QR and RQ factorizations): Tests
SGGQRFandSGGRQF - GSV (Generalized Singular Value Decomposition): Tests
SGGSVD,SGGSVP,STGSJA,SLAGS2,SLAPLL, andSLAPMT - CSD (CS decomposition): Tests
SORCSD - LSE (Constrained Linear Least Squares): Tests
SGGLSE
Each test path has a different set of inputs, but the data sets for the driver routines xEV, xES, xVX, and xSX can be concatenated in a single input file. The first line of input should contain one of the 3-character path names in columns 1-3. The number of remaining lines depends on what is found on the first line.
The number of matrix types used in testing is often controllable from the input file. The number of matrix types for each path, and the test routine that describes them, is as follows β okay bestie, here's the whole list, no skips:
SHSorNEPβ 21 types βSCHKHSSSTorSEPβ 21 types βSCHKST(routines); 18 types βSDRVST(drivers)SBDorSVDβ 16 types βSCHKBD(routines); 5 types βSDRVBD(drivers)SEVβ 21 types βSDRVEVSESβ 21 types βSDRVESSVXβ 21 types βSDRVVXSSXβ 21 types βSDRVSXSGGβ 26 types βSCHKGG(routines)SGSβ 26 types βSDRGESSGXβ 5 types βSDRGSXSGVβ 26 types βSDRGEVSXVβ 2 types βSDRGVXSSGβ 21 types βSDRVSGSSBβ 15 types βSCHKSBSBBβ 15 types βSCHKBBSECβ-βSCHKECSBLβ-βSCHKBLSBKβ-βSCHKBKSGLβ-βSCHKGLSGKβ-βSCHKGKGLMβ 8 types βSCKGLMGQRβ 8 types βSCKGQRGSVβ 8 types βSCKGSVCSDβ 3 types βSCKCSDLSEβ 8 types βSCKLSE
π¦ NEP input file π¦
- line 2:
NN, INTEGER β Number of values of N. - line 3:
NVAL, INTEGER array, dimension (NN) β The values for the matrix dimension N. - line 4:
NPARMS, INTEGER β Number of values of the parametersNB,NBMIN,NX,NS, andMAXB. - line 5:
NBVAL, INTEGER array, dimension (NPARMS) β The values for the blocksizeNB. - line 6:
NBMIN, INTEGER array, dimension (NPARMS) β The values for the minimum blocksizeNBMIN. - line 7:
NXVAL, INTEGER array, dimension (NPARMS) β The values for the crossover pointNX. - line 8:
INMIN, INTEGER array, dimension (NPARMS) βLAHQRvsTTQREcrossover point,>= 11 - line 9:
INWIN, INTEGER array, dimension (NPARMS) β recommended deflation window size - line 10:
INIBL, INTEGER array, dimension (NPARMS) β nibble crossover point - line 11:
ISHFTS, INTEGER array, dimension (NPARMS) β number of simultaneous shifts) - line 12:
IACC22, INTEGER array, dimension (NPARMS) β select structured matrix multiply: 0, 1 or 2) - line 13:
THRESHβ Threshold value for the test ratios. Information will be printed about each test for which the test ratio is greater than or equal to the threshold. To have all of the test ratios printed, useTHRESH = 0.0. - line 14:
NEWSD, INTEGER β A code indicating how to set the random number seed.= 0: Set the seed to a default value before each run= 1: Initialize the seed to a default value only before the first run= 2: Like 1, but use the seed values on the next line
If line 14 was 2:
- line 15: INTEGER array, dimension (4) β Four integer values for the random number seed.
- lines 15-EOF: The remaining lines occur in sets of 1 or 2 and allow the user to specify the matrix types. Each line contains a 3-character path name in columns 1-3, and the number of matrix types must be the first nonblank item in columns 4-80. If the number of matrix types is at least 1 but is less than the maximum number of possible types, a second line will be read to get the numbers of the matrix types to be used. For example,
NEP 21
requests all of the matrix types for the nonsymmetric eigenvalue problem, while
NEP 4
9 10 11 12
requests only matrices of type 9, 10, 11, and 12.
The valid 3-character path names are 'NEP' or 'SHS' for the nonsymmetric eigenvalue routines.
π SEP or SSG input file π
- line 2:
NN, INTEGER β Number of values of N. - line 3:
NVAL, INTEGER array, dimension (NN) β The values for the matrix dimension N. - line 4:
NPARMS, INTEGER β Number of values of the parametersNB,NBMIN, andNX. - line 5:
NBVAL, INTEGER array, dimension (NPARMS) β The values for the blocksizeNB. - line 6:
NBMIN, INTEGER array, dimension (NPARMS) β The values for the minimum blocksizeNBMIN. - line 7:
NXVAL, INTEGER array, dimension (NPARMS) β The values for the crossover pointNX. - line 8:
THRESHβ Threshold value for the test ratios. Information will be printed about each test for which the test ratio is greater than or equal to the threshold. - line 9:
TSTCHK, LOGICAL β Flag indicating whether or not to test the LAPACK routines. - line 10:
TSTDRV, LOGICAL β Flag indicating whether or not to test the driver routines. - line 11:
TSTERR, LOGICAL β Flag indicating whether or not to test the error exits for the LAPACK routines and driver routines. - line 12:
NEWSD, INTEGER β A code indicating how to set the random number seed.= 0: Set the seed to a default value before each run= 1: Initialize the seed to a default value only before the first run= 2: Like 1, but use the seed values on the next line
If line 12 was 2:
- line 13: INTEGER array, dimension (4) β Four integer values for the random number seed.
- lines 13-EOF: Lines specifying matrix types, as for NEP. The 3-character path names are
'SEP'or'SST'for the symmetric eigenvalue routines and driver routines, and'SSG'for the routines for the symmetric generalized eigenvalue problem.
π§ SVD input file π§
- line 2:
NN, INTEGER β Number of values of M and N. - line 3:
MVAL, INTEGER array, dimension (NN) β The values for the matrix row dimension M. - line 4:
NVAL, INTEGER array, dimension (NN) β The values for the matrix column dimension N. - line 5:
NPARMS, INTEGER β Number of values of the parameterNB,NBMIN,NX, andNRHS. - line 6:
NBVAL, INTEGER array, dimension (NPARMS) β The values for the blocksizeNB. - line 7:
NBMIN, INTEGER array, dimension (NPARMS) β The values for the minimum blocksizeNBMIN. - line 8:
NXVAL, INTEGER array, dimension (NPARMS) β The values for the crossover pointNX. - line 9:
NSVAL, INTEGER array, dimension (NPARMS) β The values for the number of right hand sidesNRHS. - line 10:
THRESHβ Threshold value for the test ratios. Information will be printed about each test for which the test ratio is greater than or equal to the threshold. - line 11:
TSTCHK, LOGICAL β Flag indicating whether or not to test the LAPACK routines. line 12: TSTDRV, LOGICALβ a flag saying whether we're testing the driver routines or not. that's it, that's the whole job πline 13: TSTERR, LOGICALβ a flag for whether or not we test the error exits for the LAPACK routines and the driver routines.line 14: NEWSD, INTEGERβ a code that tells you how to set the random number seed π²= 0: set the seed to a default value before each run= 1: initialize the seed to a default value only before the first run= 2: like1, but use the seed values on the next line
If line 14 was 2:
line 15: INTEGER array, dimension (4)β four integer values for the random number seed.lines 15-EOFβ lines specifying matrix types, as for NEP. The 3-character path names are'SVD'or'SBD', and they cover both the SVD routines and the SVD driver routines. π¦
SEV and SES data files: β¨
line 1β'SEV'or'SES'sitting in columns 1 to 3.line 2: NSIZES, INTEGERβ number of sizes of matrices to use. Should be at least 0 and at most 20. IfNSIZES = 0, no testing is done β but heads up, the remaining 3 lines are still read, so don't go deleting them bestie πline 3: NN, INTEGER array, dimension(NSIZES)β dimensions of the matrices to be tested.line 4: NB, NBMIN, NX, NS, NBCOL, INTEGERsβ these integer parameters decide how blocking gets done (peepILAENVfor the details π)NB: block sizeNBMIN: minimum block sizeNX: minimum dimension for blockingNS: number of shifts inxHSEQRNBCOL: minimum column dimension for blocking
line 5: THRESH, REALβ the test threshold that computed residuals get compared against. Should generally land in the range from10.to20.If it's0., every single test case's data gets printed. yes, all of it. every crumb. πline 6: TSTERR, LOGICALβ flag for whether or not to test the error exits.line 7: NEWSD, INTEGERβ the seed-setting code again:= 0: set the seed to a default value before each run= 1: initialize the seed to a default value only before the first run= 2: like1, but use the seed values on the next line
If line 7 was 2:
line 8: INTEGER array, dimension (4)β four integer values for the random number seed.lines 9 and followingβ lines specifying matrix types, as for NEP. The 3-character path name is'SEV'to testSGEEV, or'SES'to testSGEES. pick your fighter π
The SVX data has two parts. Part one is identical to SEV, and part two is all about test matrices with precomputed solutions. iconic, honestly.
line 1β'SVX'in columns 1-3.line 2: NSIZES, INTEGERβ ifNSIZES = 0, no testing of the randomly generated examples happens, but any precomputed examples do get tested.line 3: NN, INTEGER array, dimension(NSIZES)line 4: NB, NBMIN, NX, NS, NBCOL, INTEGERsline 5: THRESH, REALline 6: TSTERR, LOGICALline 7: NEWSD, INTEGER
If line 7 was 2:
line 8: INTEGER array, dimension (4)lines 9 and followingβ the first line holds'SVX'in columns 1-3, followed by the number of matrix types, possibly with a second line to specify certain matrix types. If the number of matrix types= 0, no testing of randomly generated examples is done β but any precomputed examples are still tested.remaining linesβ each matrix is stored on1+2*Nlines, whereNis its dimension:- The first line contains the dimension (a single integer).
- The next
Nlines contain the matrix, one row per line. - The last
Nlines correspond to each eigenvalue, and each one contains 4 real values: the real part of the eigenvalue, the imaginary part of the eigenvalue, the reciprocal condition number of the eigenvalues, and the reciprocal condition number of the eigenvector.
The end of data is indicated by dimensionN=0. Even if no data is to be tested, there must be at least one line containingN=0. no skipping, it's mandatory π«Ά
The SSX data is like SVX. First part identical to SEV, second part = test matrices with precomputed solutions. twins fr fr π
line 1β'SSX'in columns 1-3.line 2: NSIZES, INTEGERβ ifNSIZES = 0, no testing of randomly generated examples is done, but any precomputed examples are tested.line 3: NN, INTEGER array, dimension(NSIZES)line 4: NB, NBMIN, NX, NS, NBCOL, INTEGERsline 5: THRESH, REALline 6: TSTERR, LOGICALline 7: NEWSD, INTEGER
If line 7 was 2:
line 8: INTEGER array, dimension (4)lines 9 and followingβ the first line contains'SSX'in columns 1-3, followed by the number of matrix types, possibly with a second line to specify certain matrix types. If the number of matrix types= 0, no testing of randomly generated examples is done, but any precomputed examples are tested.remaining linesβ each matrix is stored on3+Nlines, whereNis its dimension:- The first line contains the dimension
Nand the dimensionMof an invariant subspace. - The second line contains
Mintegers, identifying the eigenvalues in the invariant subspace (by their position in a list of eigenvalues ordered by increasing real part). - The next
Nlines contain the matrix. - The last line contains the reciprocal condition number for the average of the selected eigenvalues, and the reciprocal condition number for the corresponding right invariant subspace.
- The first line contains the dimension
The end of data is indicated by a line containingN=0andM=0. Even if no data is to be tested, there must be at least one line containingN=0andM=0. mandatory again, sorry bestie π
SGG input file: π
line 2: NN, INTEGERβ number of values ofN.line 3: NVAL, INTEGER array, dimension (NN)β the values for the matrix dimensionN.line 4: NPARMS, INTEGERβ number of values of the parametersNB,NBMIN,NS,MAXB, andNBCOL.line 5: NBVAL, INTEGER array, dimension (NPARMS)β the values for the blocksizeNB.line 6: NBMIN, INTEGER array, dimension (NPARMS)β the values forNBMIN, the minimum row dimension for blocks.line 7: NSVAL, INTEGER array, dimension (NPARMS)β the values for the number of shifts.line 8: MXBVAL, INTEGER array, dimension (NPARMS)β the values forMAXB, used in determining minimum blocksize.line 9: IACC22, INTEGER array, dimension (NPARMS)β select structured matrix multiply:1or2)line 10: NBCOL, INTEGER array, dimension (NPARMS)β the values forNBCOL, the minimum column dimension for blocks.line 11: THRESHβ threshold value for the test ratios. Info gets printed about each test whose test ratio is greater than or equal to the threshold.line 12: TSTCHK, LOGICALβ flag for whether or not to test the LAPACK routines.line 13: TSTDRV, LOGICALβ flag for whether or not to test the driver routines.line 14: TSTERR, LOGICALβ flag for whether or not to test the error exits for the LAPACK routines and driver routines.line 15: NEWSD, INTEGERβ the seed code, one more time:= 0: set the seed to a default value before each run= 1: initialize the seed to a default value only before the first run= 2: like1, but use the seed values on the next line
If line 15 was 2:
line 16: INTEGER array, dimension (4)β four integer values for the random number seed.lines 17-EOFβ lines specifying matrix types, as for NEP. The 3-character path name is'SGG'for the generalized eigenvalue problem routines and driver routines. π
SGS and SGV input files:
line 1β'SGS'or'SGV'in columns 1 to 3.line 2: NN, INTEGERβ number of values ofN.line 3: NVAL, INTEGER array, dimension(NN)β dimensions of matrices to be tested.line 4: NB, NBMIN, NX, NS, NBCOL, INTEGERsβ these integer parameters determine how blocking is done (seeILAENVfor details)NB: block sizeNBMIN: minimum block sizeNX: minimum dimension for blockingNS: number of shifts inxHGEQRNBCOL: minimum column dimension for blocking
line 5: THRESH, REALβ the test threshold that computed residuals get compared against. Should generally be in the range from10.to20.If it's0., all test case data will be printed.line 6: TSTERR, LOGICALβ flag indicating whether or not to test the error exits.line 7: NEWSD, INTEGERβ and the seed code, as always:= 0: set the seed to a default value before each run= 1: initialize the seed to a default value only before the first run= 2: like1, but use the seed values on the next line
If line 17 was 2:
line 7: INTEGER array, dimension (4)β four integer values for the random number seed.
lines 7-EOF: lines specifying matrix types, exactly like for NEP. The 3-character path name is 'SGS' for the generalized eigenvalue problem routines and driver routines. π
β¨ SXV input files β¨
- line 1:
'SXV'chillin' in columns 1 to 3. - line 2:
N, INTEGER β the value of N, bestie. - line 3:
NB,NBMIN,NX,NS,NBCOL, all INTEGERs β these lil ints decide how the blocking gets done (go peepILAENVfor the full tea):NBβ block sizeNBMINβ minimum block sizeNXβ minimum dimension for blockingNSβ number of shifts inxHGEQRNBCOLβ minimum column dimension for blocking
- line 4:
THRESH, REAL β the test threshold your computed residuals get compared against. Lowkey should be in the range from 10. to 20. You'll get info printed about every test whose test ratio is greater than or equal to the threshold. - line 5:
TSTERR, LOGICAL β flag that says whether or not to test the error exits for the LAPACK routines and driver routines. - line 6:
NEWSD, INTEGER β a code telling us how to set the random number seed:= 0: set the seed to a default value before each run= 1: initialize the seed to a default value only before the first run= 2: like 1, but use the seed values on the next line
If line 6 was 2:
- line 7: INTEGER array, dimension
(4)β four integer values for the random number seed.
If line 2 was 0:
- line 7-EOF: precomputed examples are tested.
- remaining lines: each example is stored on
3+2*Nlines, where N is its dimension. The first line contains the dimension (a single integer). The next N lines contain the matrix A, one row per line. The next N lines contain the matrix B. The next line contains the reciprocals of the eigenvalue condition numbers. The last line contains the reciprocals of the eigenvector condition numbers. The end of data is indicated by dimension N=0. Even if no data is to be tested, there must be at least one line containing N=0. π«Ά
πΈ SGX input files πΈ
- line 1:
'SGX'in columns 1 to 3. - line 2:
N, INTEGER β the value of N. - line 3:
NB,NBMIN,NX,NS,NBCOL, INTEGERs β these integer parameters determine how blocking is done (seeILAENVfor details):NBβ block sizeNBMINβ minimum block sizeNXβ minimum dimension for blockingNSβ number of shifts inxHGEQRNBCOLβ minimum column dimension for blocking
- line 4:
THRESH, REAL β the test threshold against which computed residuals are compared. Should generally be in the range from 10. to 20. Information will be printed about each test for which the test ratio is greater than or equal to the threshold. - line 5:
TSTERR, LOGICAL β flag indicating whether or not to test the error exits for the LAPACK routines and driver routines. - line 6:
NEWSD, INTEGER β a code indicating how to set the random number seed:= 0: set the seed to a default value before each run= 1: initialize the seed to a default value only before the first run= 2: like 1, but use the seed values on the next line
If line 6 was 2:
- line 7: INTEGER array, dimension
(4)β four integer values for the random number seed.
If line 2 was 0:
- line 7-EOF: precomputed examples are tested.
- remaining lines: each example is stored on
3+2*Nlines, where N is its dimension. The first line contains the dimension (a single integer). The next line contains an integer k such that only the last k eigenvalues will be selected and appear in the leading diagonal blocks of $A$ and $B$. The next N lines contain the matrix A, one row per line. The next N lines contain the matrix B. The last line contains the reciprocal of the eigenvalue cluster condition number and the reciprocal of the deflating subspace (associated with the selected eigencluster) condition number. The end of data is indicated by dimension N=0. Even if no data is to be tested, there must be at least one line containing N=0. π
π SSB input file π
- line 2:
NN, INTEGER β number of values of N. - line 3:
NVAL, INTEGER array, dimension(NN)β the values for the matrix dimension N. - line 4:
NK, INTEGER β number of values of K. - line 5:
KVAL, INTEGER array, dimension(NK)β the values for the matrix dimension K. - line 6:
THRESHβ threshold value for the test ratios. Information will be printed about each test for which the test ratio is greater than or equal to the threshold. - line 7:
NEWSD, INTEGER β a code indicating how to set the random number seed:= 0: set the seed to a default value before each run= 1: initialize the seed to a default value only before the first run= 2: like 1, but use the seed values on the next line
If line 7 was 2:
- line 8: INTEGER array, dimension
(4)β four integer values for the random number seed. - lines 8-EOF: lines specifying matrix types, as for NEP. The 3-character path name is
'SSB'. β
π¦ SBB input file π¦
- line 2:
NN, INTEGER β number of values of M and N. - line 3:
MVAL, INTEGER array, dimension(NN)β the values for the matrix row dimension M. - line 4:
NVAL, INTEGER array, dimension(NN)β the values for the matrix column dimension N. - line 4:
NK, INTEGER β number of values of K. - line 5:
KVAL, INTEGER array, dimension(NK)β the values for the matrix bandwidth K. - line 6:
NPARMS, INTEGER β number of values of the parameterNRHS. - line 7:
NSVAL, INTEGER array, dimension(NPARMS)β the values for the number of right hand sidesNRHS. - line 8:
THRESHβ threshold value for the test ratios. Information will be printed about each test for which the test ratio is greater than or equal to the threshold. - line 9:
NEWSD, INTEGER β a code indicating how to set the random number seed:= 0: set the seed to a default value before each run= 1: initialize the seed to a default value only before the first run= 2: like 1, but use the seed values on the next line
If line 9 was 2:
- line 10: INTEGER array, dimension
(4)β four integer values for the random number seed. - lines 10-EOF: lines specifying matrix types, as for SVD. The 3-character path name is
'SBB'. π₯
π SEC input file π
- line 2:
THRESH, REAL β threshold value for the test ratios. Information will be printed about each test for which the test ratio is greater than or equal to the threshold. - lines 3-EOF: input for testing the eigencondition routines consists of a set of specially constructed test cases and their solutions. The data format is not intended to be modified by the user. π
π§ SBL and SBK input files π§
- line 1:
'SBL'in columns 1-3 to testSGEBAL, or'SBK'in columns 1-3 to testSGEBAK. - The remaining lines consist of specially constructed test cases.
π SGL and SGK input files π
- line 1:
'SGL'in columns 1-3 to testSGGBAL, or'SGK'in columns 1-3 to testSGGBAK. - The remaining lines consist of specially constructed test cases.
π GLM data file π
- line 1:
'GLM'in columns 1 to 3. - line 2:
NN, INTEGER β number of values of M, P, and N. - line 3:
MVAL, INTEGER array, dimension(NN)β values of M (row dimension). - line 4:
PVAL, INTEGER array, dimension(NN)β values of P (row dimension). - line 5:
NVAL, INTEGER array, dimension(NN)β values of N (column dimension), note M <= N <= M+P. - line 6:
THRESH, REAL β threshold value for the test ratios. Information will be printed about each test for which the test ratio is greater than or equal to the threshold. - line 7:
TSTERR, LOGICAL β flag indicating whether or not to test the error exits for the LAPACK routines and driver routines. - line 8:
NEWSD, INTEGER β a code indicating how to set the random number seed:= 0: set the seed to a default value before each run= 1: initialize the seed to a default value only before the first run= 2: like 1, but use the seed values on the next line
If line 8 was 2:
- line 9: INTEGER array, dimension
(4)β four integer values for the random number seed. - lines 9-EOF: lines specifying matrix types, as for NEP. The 3-character path name is
'GLM'for the generalized linear regression model routines. π
π GQR data file π
- line 1:
'GQR'in columns 1 to 3. - line 2:
NN, INTEGER β number of values of M, P, and N. line 3: MVAL, INTEGER array, dimension(NN)β these are the M values, bestie. πline 4: PVAL, INTEGER array, dimension(NN)β the P values, no cap. β¨line 5: NVAL, INTEGER array, dimension(NN)β the N values, periodt. πline 6: THRESH, REALβ threshold for the test ratios. Information gets printed for every test where the ratio is greater than or equal to this threshold. π¨line 7: TSTERR, LOGICALβ flag that says whether or not to test the error exits for the LAPACK routines and driver routines. πline 8: NEWSD, INTEGERβ a code that tells you how to set the random number seed:= 0: Set the seed to a default value before each run= 1: Initialize the seed to a default value only before the first run= 2: Like 1, but use the seed values on the next line
Ifline 8was2: then you get this next bit, bestie.
line 9: INTEGER array, dimension (4)β four integer values for the random number seed. π²lines 9-EOF: Lines specifying matrix types, as for NEP.β and the 3-character path name is'GQR'for the generalized QR and RQ routines. π₯
π GSV data file π
line 1: 'GSV' in columns 1 to 3.β that's the vibe check at the top, literally. πΈline 2: NN, INTEGERβ number of values of M, P, and N. β¨line 3: MVAL, INTEGER array, dimension(NN)β values of M (row dimension). πline 4: PVAL, INTEGER array, dimension(NN)β values of P (row dimension). πline 5: NVAL, INTEGER array, dimension(NN)β values of N (column dimension). π¦line 6: THRESH, REALβ threshold value for the test ratios. Info gets printed about each test for which the test ratio is greater than or equal to the threshold. π¨line 7: TSTERR, LOGICALβ flag indicating whether or not to test the error exits for the LAPACK routines and driver routines. πline 8: NEWSD, INTEGERβ a code indicating how to set the random number seed.= 0: Set the seed to a default value before each run= 1: Initialize the seed to a default value only before the first run= 2: Like 1, but use the seed values on the next line
Ifline 8was2:
line 9: INTEGER array, dimension (4)β four integer values for the random number seed. π²lines 9-EOF: Lines specifying matrix types, as for NEP.β the 3-character path name is'GSV'for the generalized SVD routines. π
π« CSD data file π«
line 1: 'CSD' in columns 1 to 3.β iconic. πline 2: NM, INTEGERβ number of values of M, P, and N. β¨line 3: MVAL, INTEGER array, dimension(NM)β values of M (row and column dimension of orthogonal matrix). πline 4: PVAL, INTEGER array, dimension(NM)β values of P (row dimension of top-left block). πline 5: NVAL, INTEGER array, dimension(NM)β values of N (column dimension of top-left block). π¦line 6: THRESH, REALβ threshold value for the test ratios. Info gets printed about each test for which the test ratio is greater than or equal to the threshold. π¨line 7: TSTERR, LOGICALβ flag indicating whether or not to test the error exits for the LAPACK routines and driver routines. πline 8: NEWSD, INTEGERβ a code indicating how to set the random number seed.= 0: Set the seed to a default value before each run= 1: Initialize the seed to a default value only before the first run= 2: Like 1, but use the seed values on the next line
Ifline 8was2:
line 9: INTEGER array, dimension (4)β four integer values for the random number seed. π²lines 9-EOF: Lines specifying matrix types, as for NEP.β the 3-character path name is'CSD'for the CSD routine. π
π§ LSE data file π§
line 1: 'LSE' in columns 1 to 3.β main character energy. πline 2: NN, INTEGERβ number of values of M, P, and N. β¨line 3: MVAL, INTEGER array, dimension(NN)β values of M. πline 4: PVAL, INTEGER array, dimension(NN)β values of P. πline 5: NVAL, INTEGER array, dimension(NN)β values of N, noteP <= N <= P+M. π¦line 6: THRESH, REALβ threshold value for the test ratios. Info gets printed about each test for which the test ratio is greater than or equal to the threshold. π¨line 7: TSTERR, LOGICALβ flag indicating whether or not to test the error exits for the LAPACK routines and driver routines. πline 8: NEWSD, INTEGERβ a code indicating how to set the random number seed.= 0: Set the seed to a default value before each run= 1: Initialize the seed to a default value only before the first run= 2: Like 1, but use the seed values on the next line
Ifline 8was2:
line 9: INTEGER array, dimension (4)β four integer values for the random number seed. π²lines 9-EOF: Lines specifying matrix types, as for NEP.β the 3-character path name is'GSV'for the generalized SVD routines. π
NMAX is currently set to 132 and must be at least 12 for some of the precomputed examples, and LWORK = NMAX*(5*NMAX+5)+1 in the parameter statements below. For SVD, we assume NRHS may be as big as N. The parameter NEED is set to 14 to allow for 14 N-by-N matrices for SGG. π₯
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 1039 of file schkee.F. π
π Author π
Generated automatically by Doxygen for LAPACK from the source code. π€
And that's the whole vibe, bestie! π This was LAPACK version 3.12.0, dated Tue Nov 28 2023 12:08:42, from the file TESTING/EIG/schkee.F(3). Mwah! πβ¨
π¦ where this page comes from
- centos-stream 10: lapack 3.12.0-7.el10 Β· BSD-3-Clause-Open-MPI
- fedora 43: lapack 3.12.0-10.fc43 Β· BSD-3-Clause-Open-MPI
- fedora 43: lapack 3.12.0-10.fc43 Β· BSD-3-Clause-Open-MPI
- fedora 44: lapack 3.12.0-11.fc44 Β· BSD-3-Clause-Open-MPI
- fedora 44: lapack 3.12.0-11.fc44 Β· BSD-3-Clause-Open-MPI
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash on : a language model's rewrite of the original, under its terms