complex_eig(3)
π¦ NAME π¦
complex_eig - complex
okay so hi bestie π you just opened complex_eig(3), straight out of the LAPACK manual β yes, the whole complex_eig situation. And ngl this page is long, so it comes in 4 parts and you're reading part 1 right now. Buckle up, it's giving documentation but make it cute. π
β¨ SYNOPSIS β¨
π Functions π
program cchkeeβCCHKEE
πΈ Detailed Description πΈ
This is the group of complex LAPACK TESTING EIG routines. Periodt. That's literally the whole vibe β it's the test squad for the complex eigenvalue problems, no more no less. π«Ά
π Function Documentation π
π§ program cchkee π§
π CCHKEE π
π Purpose π
CCHKEE tests the COMPLEX LAPACK subroutines for the matrix eigenvalue problem. The test paths in this version are
- NEP (Nonsymmetric Eigenvalue Problem): tests
CGEHRD,CUNGHR,CHSEQR,CTREVC,CHSEIN, andCUNMHR - SEP (Hermitian Eigenvalue Problem): tests
CHETRD,CUNGTR,CSTEQR,CSTERF,CSTEIN,CSTEDC, and the driversCHEEV(X),CHBEV(X),CHPEV(X),CHEEVD,CHBEVD,CHPEVD - SVD (Singular Value Decomposition): tests
CGEBRD,CUNGBR, andCBDSQR, and the driversCGESVD,CGESDD - CEV (Nonsymmetric Eigenvalue/eigenvector Driver): tests
CGEEV - CES (Nonsymmetric Schur form Driver): tests
CGEES - CVX (Nonsymmetric Eigenvalue/eigenvector Expert Driver): tests
CGEEVX - CSX (Nonsymmetric Schur form Expert Driver): tests
CGEESX - CGG (Generalized Nonsymmetric Eigenvalue Problem): tests
CGGHD3,CGGBAL,CGGBAK,CHGEQZ, andCTGEVC - CGS (Generalized Nonsymmetric Schur form Driver): tests
CGGES - CGV (Generalized Nonsymmetric Eigenvalue/eigenvector Driver): tests
CGGEV - CGX (Generalized Nonsymmetric Schur form Expert Driver): tests
CGGESX - CXV (Generalized Nonsymmetric Eigenvalue/eigenvector Expert Driver): tests
CGGEVX - CSG (Hermitian Generalized Eigenvalue Problem): tests
CHEGST,CHEGV,CHEGVD,CHEGVX,CHPGST,CHPGV,CHPGVD,CHPGVX,CHBGST,CHBGV,CHBGVD, andCHBGVX - CHB (Hermitian Band Eigenvalue Problem): tests
CHBTRD - CBB (Band Singular Value Decomposition): tests
CGBBRD - CEC (Eigencondition estimation): tests
CTRSYL,CTREXC,CTRSNA, andCTRSEN - CBL (Balancing a general matrix): tests
CGEBAL - CBK (Back transformation on a balanced matrix): tests
CGEBAK - CGL (Balancing a matrix pair): tests
CGGBAL - CGK (Back transformation on a matrix pair): tests
CGGBAK - GLM (Generalized Linear Regression Model): tests
CGGGLM - GQR (Generalized QR and RQ factorizations): tests
CGGQRFandCGGRQF - GSV (Generalized Singular Value Decomposition): tests
CGGSVD,CGGSVP,CTGSJA,CLAGS2,CLAPLL, andCLAPMT - CSD (CS decomposition): tests
CUNCSD - LSE (Constrained Linear Least Squares): tests
CGGLSE
Every single test path takes a different set of inputs, bestie β but lowkey there's a shortcut: the data sets for the driver routines xEV, xES, xVX, and xSX can all be concatenated into one single input file. π The first line of input should contain one of the 3-character path names in columns 1-3. How many lines come after that depends entirely on what's sitting 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:
Path name(s) Types Test routine
CHS or NEP 21 CCHKHS
CST or SEP 21 CCHKST (routines)
18 CDRVST (drivers)
CBD or SVD 16 CCHKBD (routines)
5 CDRVBD (drivers)
CEV 21 CDRVEV
CES 21 CDRVES
CVX 21 CDRVVX
CSX 21 CDRVSX
CGG 26 CCHKGG (routines)
CGS 26 CDRGES
CGX 5 CDRGSX
CGV 26 CDRGEV
CXV 2 CDRGVX
CSG 21 CDRVSG
CHB 15 CCHKHB
CBB 15 CCHKBB
CEC - CCHKEC
CBL - CCHKBL
CBK - CCHKBK
CGL - CCHKGL
CGK - CCHKGK
GLM 8 CCKGLM
GQR 8 CCKGQR
GSV 8 CCKGSV
CSD 3 CCKCSD
LSE 8 CCKLSE
π₯ NEP input file π₯
- line 2:
NN, INTEGERβ how many values ofNwe're testing. - line 3:
NVAL, INTEGER array, dimension (NN)β the actual values for the matrix dimensionN. - 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,1or2 - line 13:
THRESHβ threshold value for the test ratios. Info gets printed about each test where the test ratio is greater than or equal to the threshold. If you want every single test ratio printed, useTHRESH = 0.0. - line 14:
NEWSD, INTEGERβ a code telling it 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 come in sets of 1 or 2, and they let you, the user, specify the matrix types. Each line has a 3-character path name in columns 1-3, and the number of matrix types has to be the first nonblank item in columns 4-80. If the number of matrix types is at least 1 but less than the maximum number of possible types, a second line gets read to grab 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 'CHS' for the nonsymmetric eigenvalue routines.
π SEP or CSG 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, 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. Info gets printed about each test where the test ratio is greater than or equal to the threshold. - line 9:
TSTCHK, LOGICALβ flag saying whether or not to test the LAPACK routines. - line 10:
TSTDRV, LOGICALβ flag saying whether or not to test the driver routines. - line 11:
TSTERR, LOGICALβ flag saying whether or not to test the error exits for the LAPACK routines and driver routines. - line 12:
NEWSD, INTEGERβ a code telling it 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, same as for NEP. The valid 3-character path names are 'SEP' or 'CST' for the Hermitian eigenvalue routines and driver routines, and 'CSG' for the routines for the Hermitian generalized eigenvalue problem.
π¦ SVD 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 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. Info gets printed about each test where the test ratio is greater than or equal to the threshold. - line 11:
TSTCHK, LOGICALβ flag saying whether or not to test the LAPACK routines. - line 12:
TSTDRV, LOGICAL β a flag saying whether or not to test the driver routines. like, do we put the drivers through their paces or not π - line 13:
TSTERR, LOGICAL β a flag saying whether or not to test the error exits for the LAPACK routines and the driver routines. gotta make sure the error paths hit too, fr fr. - line 14:
NEWSD, INTEGER β a code that tells the test 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. the seed lore itself β¨
- lines 15-EOF: lines specifying matrix types, as for NEP. The 3-character path names are
'SVD'or'CBD'for both the SVD routines and the SVD driver routines.
πΈ CEV and CES data files πΈ
- line 1:
'CEV'or'CES'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 (although the remaining 3 lines are still read). yes bestie, we still read them, the file doesn't get to skip π - 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 (see ILAENV for details)NB: block sizeNBMIN: minimum block sizeNX: minimum dimension for blockingNS: number of shifts in xHSEQRNBCOL: minimum column dimension for blocking
- line 5:
THRESH, REAL β the test threshold against which computed residuals are compared. Should generally be in the range from 10. to 20. If it is 0., all test case data will be printed. 0. means we're dumping EVERYTHING, no filter, main character behaviour π«Ά - 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.
- lines 8 and following: lines specifying matrix types, as for NEP. The 3-character path name is
'CEV'to test CGEEV, or'CES'to test CGEES.
π The CVX data π
The CVX data has two parts. The first part is identical to CEV, and the second part consists of test matrices with precomputed solutions.
- line 1:
'CVX'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, INTEGERs - line 5:
THRESH, REAL - line 6:
NEWSD, INTEGER
If line 6 was 2:
- line 7: INTEGER array, dimension (4)
- lines 8 and following: the first line contains
'CVX'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 on
1+N+N**2lines, whereNis its dimension. The first line contains the dimensionNandISRT(two integers).ISRTindicates whether the lastNlines are sorted by increasing real part of the eigenvalue (ISRT=0) or by increasing imaginary part (ISRT=1). The nextN**2lines contain the matrix rowwise, one entry per line. The lastNlines correspond to each eigenvalue. Each of these lastNlines contains 4 real values: the real part of the eigenvalues, the imaginary part of the eigenvalue, the reciprocal condition number of the eigenvalues, and the reciprocal condition number of the vector 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 cap, at least oneN=0or the whole thing is invalid, periodt π
π¦ The CSX data π¦
The CSX data is like CVX. The first part is identical to CEV, and the second part consists of test matrices with precomputed solutions.
- line 1:
'CSX'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, INTEGERs - line 5:
THRESH, REAL - line 6:
NEWSD, INTEGER
If line 6 was 2:
- line 7: INTEGER array, dimension (4)
- lines 8 and following: the first line contains
'CSX'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 on
3+N**2lines, whereNis its dimension. The first line contains the dimensionN, the dimensionMof an invariant subspace, andISRT. The second line containsMintegers, identifying the eigenvalues in the invariant subspace (by their position in a list of eigenvalues ordered by increasing real part (ifISRT=0) or by increasing imaginary part (ifISRT=1)). The nextN**2lines contain the matrix rowwise. 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=0,M=0, andISRT = 0. Even if no data is to be tested, there must be at least one line containingN=0,M=0andISRT=0. again: that final line is MANDATORY, do not skip it, it's giving non-negotiable π
π CGG 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,NBCOL,NS, andMAXB. - 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: like 1, 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
'CGG'for the generalized eigenvalue problem routines and driver routines. generalized eigenvalue problem gang rise up π
π CGS and CGV input files π
- line 1:
'CGS'or'CGV'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 (see ILAENV for details)NB: block sizeNBMIN: minimum block sizeNX: minimum dimension for blockingNS: number of shifts in xHGEQRNBCOL: minimum column dimension for blocking
- line 5:
THRESH, REAL β the test threshold against which computed residuals are compared. Should generally be in the range from 10. to 20. If it is 0., 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: like 1, 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 'CGS' for the generalized eigenvalue problem routines and driver routines. π
π CGX input file: π
- line 1:
'CGX'in columns 1 to 3. β¨ - line 2:
N,INTEGERβ literally the value ofN, bestie. π - line 3:
NB,NBMIN,NX,NS,NBCOL,INTEGERsβ these integer parameters determine how blocking is done (peepILAENVfor the deets, iykyk):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:
INTEGERarray, dimension(4)β four integer values for the random number seed. π²
- line 7:
- If line 2 was 0:
- line 7-EOF: Precomputed examples are tested.
- remaining lines : Each example is stored on
3+2*N*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 nextN*Nlines contain the matrixA, one element per line. The nextN*Nlines 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. π
π¦ CXV input files: π¦
- line 1:
'CXV'in columns 1 to 3. β¨ - line 2:
N,INTEGERβ the value ofN. π - 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 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:
INTEGERarray, dimension(4)β four integer values for the random number seed. π²
- line 7:
- If line 2 was 0:
- line 7-EOF: Precomputed examples are tested.
- remaining lines : Each example is stored on
3+2*N*Nlines, whereNis its dimension. The first line contains the dimension (a single integer). The nextN*Nlines contain the matrixA, one element per line. The nextN*Nlines 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. π
πΈ CHB input file: πΈ
- line 2:
NN,INTEGERβ number of values ofN. π’ - line 3:
NVAL,INTEGERarray, dimension(NN)β the values for the matrix dimensionN. π - line 4:
NK,INTEGERβ number of values ofK. π’ - line 5:
KVAL,INTEGERarray, 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:= 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:
INTEGERarray, dimension(4)β four integer values for the random number seed. π²
- line 8:
- lines 8-EOF: Lines specifying matrix types, as for NEP. The 3-character path name is
'CHB'. β¨
π§ CBB input file: π§
- line 2:
NN,INTEGERβ number of values ofMandN. π’ - line 3:
MVAL,INTEGERarray, dimension(NN)β the values for the matrix row dimensionM. π - line 4:
NVAL,INTEGERarray, dimension(NN)β the values for the matrix column dimensionN. π - line 4:
NK,INTEGERβ number of values ofK. π’ - line 5:
KVAL,INTEGERarray, dimension(NK)β the values for the matrix bandwidthK. π - line 6:
NPARMS,INTEGERβ number of values of the parameterNRHS. π’ - line 7:
NSVAL,INTEGERarray, 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:
INTEGERarray, dimension(4)β four integer values for the random number seed. π²
- line 10:
- lines 10-EOF: Lines specifying matrix types, as for SVD. The 3-character path name is
'CBB'. β¨
π CEC 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. π«
π CBL and CBK input files: π
- line 1:
'CBL'in columns 1-3 to testCGEBAL, or'CBK'in columns 1-3 to testCGEBAK. β¨ - The remaining lines consist of specially constructed test cases. π§ͺ
π¦ CGL and CGK input files: π¦
- line 1:
'CGL'in columns 1-3 to testCGGBAL, or'CGK'in columns 1-3 to testCGGBAK. β¨ - 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,INTEGERarray, dimension(NN)β values ofM(row dimension). π - line 4:
PVAL,INTEGERarray, dimension(NN)β values ofP(row dimension). π - line 5:
NVAL,INTEGERarray, 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:
INTEGERarray, dimension(4)β four integer values for the random number seed. π²
- line 9:
- 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. π’ line 3:MVAL,INTEGERarray, dimension(NN) β the values of M πline 4:PVAL,INTEGERarray, dimension(NN) β the values of P β¨line 5:NVAL,INTEGERarray, dimension(NN) β the values of N πΈline 6:THRESH,REALβ the threshold value for the test ratios. bestie, if a test ratio is greater than or equal to this threshold, you're gonna get information printed about that test, no cap πline 7:TSTERR,LOGICALβ a flag that says whether or not to test the error exits for the LAPACK routines and the driver routines π§line 8:NEWSD,INTEGERβ a code telling 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
If line 8 was 2:
line 9:INTEGERarray, 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 of M, P, and Nline 3:MVAL,INTEGERarray, dimension(NN) β the values of M (row dimension)line 4:PVAL,INTEGERarray, dimension(NN) β the values of P (row dimension)line 5:NVAL,INTEGERarray, dimension(NN) β the values of N (column dimension)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β flag indicating whether or not to test the error exits for the LAPACK routines and driver routinesline 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:INTEGERarray, 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 of M, P, and Nline 3:MVAL,INTEGERarray, dimension(NM) β the values of M (row and column dimension of orthogonal matrix)line 4:PVAL,INTEGERarray, dimension(NM) β the values of P (row dimension of top-left block)line 5:NVAL,INTEGERarray, dimension(NM) β the values of N (column dimension of top-left block)line 6:THRESH,REALβ the threshold value for the test ratios. Info gets printed about each test whose test ratio is greater than or equal to the threshold πline 7:TSTERR,LOGICALβ flag saying whether or not to test the error exits for the LAPACK routines and driver routinesline 8:NEWSD,INTEGERβ a code for 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:INTEGERarray, 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 of M, P, and Nline 3:MVAL,INTEGERarray, dimension(NN) β the values of Mline 4:PVAL,INTEGERarray, dimension(NN) β the values of Pline 5:NVAL,INTEGERarray, dimension(NN) β the values of N, and heads up:P <= N <= P+Mπline 6:THRESH,REALβ the threshold value for the test ratios. Information gets printed for each test whose 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 routinesline 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:INTEGERarray, 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 π
Okay so real talk about the sizing, because this part matters: NMAX is currently set to 132 and must be at least 12 for some of the precomputed examples, and LWORK = NMAX*(5*NMAX+20) 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 CGG π₯
π 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 β we ate and left absolutely no crumbs π
This was Version 3.11.0, dated Sun Nov 27 2022, and it lives at complex_eig(3) in the LAPACK source. Go forth and test your matrix types, icon, 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