single_eig(3)
β¨ single_eig(3) β¨
okay bestie, sit down, hydrate, we're doing this π«Ά this is single_eig(3), a page from the LAPACK manual, and it's the full guide to the real LAPACK TESTING EIG crew. it is long. it is detailed. it is giving "phone book but make it linear algebra". but we're getting through it together, no cap. π
πΈ NAME πΈ
single_eig β real
that's the whole name. that's it. iconic. π
πΈ SYNOPSIS πΈ
Functions
program schkeeSCHKEE
πΈ Detailed Description πΈ
okay so. this is the group of real LAPACK TESTING EIG routines. that's the whole vibe check. β¨
πΈ Function Documentation πΈ
π program schkee π
SCHKEE
Purpose:
SCHKEE tests the REAL LAPACK subroutines for the matrix eigenvalue problem. The test paths in this version are
- 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:
SHSorNEPβ 21 types βSCHKHSSSTorSEPβ 21 types βSCHKST(routines)SSTorSEPβ 18 types βSDRVST(drivers)SBDorSVDβ 16 types βSCHKBD(routines)SBDorSVDβ 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.
- Number of values of
- line 3:
NVAL, INTEGER array, dimension (NN)- The values for the matrix dimension
N.
- The values for the matrix dimension
- line 4:
NPARMS, INTEGER- Number of values of the parameters
NB,NBMIN,NX,NS, andMAXB.
- Number of values of the parameters
- line 5:
NBVAL, INTEGER array, dimension (NPARMS)- The values for the blocksize
NB.
- The values for the blocksize
- line 6:
NBMIN, INTEGER array, dimension (NPARMS)- The values for the minimum blocksize
NBMIN.
- The values for the minimum blocksize
- line 7:
NXVAL, INTEGER array, dimension (NPARMS)- The values for the crossover point
NX.
- The values for the crossover point
- line 8:
INMIN, INTEGER array, dimension (NPARMS)- LAHQR vs TTQRE crossover point,
>= 11
- LAHQR vs TTQRE crossover point,
- 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, use
THRESH = 0.0.
- 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, use
- 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.
ngl the valid 3-character path names here are'NEP'or'SHS'for the nonsymmetric eigenvalue routines. π
π SEP or SSG input file π
- line 2:
NN, INTEGER- Number of values of
N.
- Number of values of
- line 3:
NVAL, INTEGER array, dimension (NN)- The values for the matrix dimension
N.
- The values for the matrix dimension
- line 4:
NPARMS, INTEGER- Number of values of the parameters
NB,NBMIN, andNX.
- Number of values of the parameters
- line 5:
NBVAL, INTEGER array, dimension (NPARMS)- The values for the blocksize
NB.
- The values for the blocksize
- line 6:
NBMIN, INTEGER array, dimension (NPARMS)- The values for the minimum blocksize
NBMIN.
- The values for the minimum blocksize
- line 7:
NXVAL, INTEGER array, dimension (NPARMS)- The values for the crossover point
NX.
- The values for the crossover point
- 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
MandN.
- Number of values of
- line 3:
MVAL, INTEGER array, dimension (NN)- The values for the matrix row dimension
M.
- The values for the matrix row dimension
- line 4:
NVAL, INTEGER array, dimension (NN)- The values for the matrix column dimension
N.
- The values for the matrix column dimension
- line 5:
NPARMS, INTEGER- Number of values of the parameter
NB,NBMIN,NX, andNRHS.
- Number of values of the parameter
- line 6:
NBVAL, INTEGER array, dimension (NPARMS)- The values for the blocksize
NB.
- The values for the blocksize
- line 7:
NBMIN, INTEGER array, dimension (NPARMS)- The values for the minimum blocksize
NBMIN.
- The values for the minimum blocksize
- line 8:
NXVAL, INTEGER array, dimension (NPARMS)- The values for the crossover point
NX.
- The values for the crossover point
- line 9:
NSVAL, INTEGER array, dimension (NPARMS)- The values for the number of right hand sides
NRHS.
- The values for the number of right hand sides
- 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β the flag that decides whether or not we're testing the LAPACK routines at all. π«Άline 12: TSTDRV, LOGICALβ same energy, but for the driver routines. πline 13: TSTERR, LOGICALβ and this one says whether or not we're testing the error exits for the LAPACK routines and the driver routines. Not me being thorough or anything. πline 14: NEWSD, INTEGERβ a code telling you how to set the random number seed. It's giving choose-your-own-adventure:= 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, you get a lil bonus line:
line 15: INTEGER array, dimension (4)β four integer values for the random number seed. Four. That's the whole vibe. π²lines 15-EOF:β lines specifying matrix types, as for NEP. The 3-character path names are'SVD'or'SBD', for both the SVD routines and the SVD driver routines. β¨
πΈ SEV and SES data files πΈ
line 1:β'SEV'or'SES'in columns 1 to 3. That's how we know who we are bestie.line 2: NSIZES, INTEGERβ number of sizes of matrices to use. Should be at least0and at most20. And heads up: ifNSIZES = 0, no testing is done β but the remaining 3 lines are still read. The dedication. πline 3: NN, INTEGER array, dimension(NSIZES)β 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 inxHSEQRNBCOL: 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.And if it's0., all test case data will be printed. All of it. Buckle up. π₯line 6: TSTERR, LOGICALβ flag indicating whether or not to test the error exits.line 7: NEWSD, INTEGERβ 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: 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 π¦
The SVX data has two parts. The first part is identical to SEV, and the second part consists of test matrices with precomputed solutions. Iconic two-act structure, no notes.
line 1:β'SVX'in columns 1-3.line 2: NSIZES, INTEGERβ ifNSIZES = 0, no testing of randomly generated examples is done, but any precomputed examples are tested. (The precomputed ones always show up. Main characters. π )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'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 tested.remaining lines:β each matrix is stored on1+2*Nlines, whereNis its dimension. The first line contains the dimension (a single integer). The nextNlines contain the matrix, one row per line. The lastNlines correspond to each eigenvalue. Each of these lastNlines 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. Non-negotiable, bestie. π
π The SSX data π
The SSX data is like SVX. The first part is identical to SEV, and the second part consists of test matrices with precomputed solutions.
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 dimensionNand the dimensionMof an invariant subspace. The second line containsMintegers, identifying the eigenvalues in the invariant subspace (by their position in a list of eigenvalues ordered by increasing real part). The nextNlines 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 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. Both zeros. Don't skip one. π€
π 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. Information will be printed about each test for which the test ratio is greater than or equal to the threshold. πline 12: TSTCHK, LOGICALβ flag indicating whether or not to test the LAPACK routines.line 13: TSTDRV, LOGICALβ flag indicating whether or not to test the driver routines.line 14: TSTERR, LOGICALβ flag indicating whether or not to test the error exits for the LAPACK routines and driver routines.line 15: 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: 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 against which computed residuals are compared. Should generally be in the range from10.to20.If it is0., 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β 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: 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, as for NEP. The 3-character path name is 'SGS' for the generalized eigenvalue problem routines and driver routines. π
periodt.
SXV input files π¦
line 1:'SXV'in columns 1 to 3. okay so first line is literally just the tag. β¨line 2:N, INTEGERβ bestie, this is the value ofN. πline 3:NB, NBMIN, NX, NS, NBCOL, INTEGERsβ these integer parameters determine how blocking is done (seeILAENVfor details). the way that blocking gets configured, iconic. π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 from10.to20.Information will be printed about each test for which the test ratio is greater than or equal to the threshold. no cap, keep it in that range. π₯line 5:TSTERR, LOGICALβ flag indicating whether or not to test the error exits for the LAPACK routines and driver routines. vibe check on error handling. πline 6:NEWSD, INTEGERβ a code indicating how to set the random number seed. choose your seed destiny π= 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, whereNis its dimension. The first line contains the dimension (a single integer). The nextNlines contain the matrixA, one row per line. The nextNlines contain the matrixB. 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 dimensionN=0. Even if no data is to be tested, there must be at least one line containingN=0. π the way that there's always aN=0line, even if empty, is so real.
SGX input files π
line 1:'SGX'in columns 1 to 3. β¨line 2:N, INTEGERβ value ofN. πline 3:NB, NBMIN, NX, NS, NBCOL, INTEGERsβ these integer parameters determine how blocking is done (seeILAENVfor details). slay π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 from10.to20.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, whereNis its dimension. The first line contains the dimension (a single integer). The next line contains an integerksuch that only the lastkeigenvalues will be selected and appear in the leading diagonal blocks of$A$and$B$. The nextNlines contain the matrixA, one row per line. The nextNlines contain the matrixB. 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 dimensionN=0. Even if no data is to be tested, there must be at least one line containingN=0. π it's giving extra credit for the eigencluster details.
SSB 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:NK, INTEGERβ number of values ofK. πline 5:KVAL, INTEGER array, dimension (NK)β the values for the matrix dimensionK. π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. main character energy π= 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 ofMandN. πline 3:MVAL, INTEGER array, dimension (NN)β the values for the matrix row dimensionM. β¨line 4:NVAL, INTEGER array, dimension (NN)β the values for the matrix column dimensionN. πline 4:NK, INTEGERβ number of values ofK. π not the line number glitching, that's literally what it saysline 5:KVAL, INTEGER array, dimension (NK)β the values for the matrix bandwidthK. π₯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. periodt, don't touch it. π§
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 ofM,P, andN. πline 3:MVAL, INTEGER array, dimension(NN)β values ofM(row dimension). πline 4:PVAL, INTEGER array, dimension(NN)β values ofP(row dimension). πline 5:NVAL, INTEGER array, dimension(NN)β values ofN(column dimension), noteM <= 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 ofM,P, andN. π
okay bestie, we're still in the data-file zone and you already know the vibe. no re-intro, just picking up at line 3 of the GQR file like nothing happened. π
line 3:MVAL,INTEGER array,dimension(NN)β the values ofM. πline 4:PVAL,INTEGER array,dimension(NN)β the values ofP. β¨line 5:NVAL,INTEGER array,dimension(NN)β the values ofN. πΈline 6:THRESH,REALβ the 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β a 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: like1, 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'GQR'for the generalized QR and RQ routines. π₯
GSV data file π¦
line 1:'GSV'in columns 1 to 3. πline 2:NN,INTEGERβ number of values ofM,P, andN. π’line 3:MVAL,INTEGER array,dimension(NN)β values ofM(row dimension). πline 4:PVAL,INTEGER array,dimension(NN)β values ofP(row dimension). πline 5:NVAL,INTEGER array,dimension(NN)β values ofN(column dimension). π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: like1, 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'GSV'for the generalized SVD routines. π₯
CSD data file π
line 1:'CSD'in columns 1 to 3. πline 2:NM,INTEGERβ number of values ofM,P, andN. π’line 3:MVAL,INTEGER array,dimension(NM)β values ofM(row and column dimension of orthogonal matrix). π§line 4:PVAL,INTEGER array,dimension(NM)β values ofP(row dimension of top-left block). π§±line 5:NVAL,INTEGER array,dimension(NM)β values ofN(column dimension of top-left block). π§±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: like1, 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'CSD'for the CSD routine. π₯
LSE data file π
line 1:'LSE'in columns 1 to 3. πline 2:NN,INTEGERβ number of values ofM,P, andN. π’line 3:MVAL,INTEGER array,dimension(NN)β values ofM. πline 4:PVAL,INTEGER array,dimension(NN)β values ofP. πline 5:NVAL,INTEGER array,dimension(NN)β values ofN, noteP <= N <= P+M. β οΈ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: like1, 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'GSV'for the generalized SVD routines. π₯
π big brain parameter tea:NMAXis currently set to132and must be at least12for some of the precomputed examples, andLWORK = NMAX*(5*NMAX+5)+1in the parameter statements below. for SVD, we assumeNRHSmay be as big asN. the parameterNEEDis set to14to allow for14N-by-Nmatrices for SGG. π
Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Author π€
Generated automatically by Doxygen for LAPACK from the source code.
and that's the whole page, bestie. Version 3.11.0, Sun Nov 27 2022, single_eig(3) β periodt. π
β¨
π¦ where this page comes from
- debian bookworm: liblapack-doc 3.11.0-2 Β· BSD-3-Clause
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash on : a language model's rewrite of the original, under its terms