TESTING/MATGEN/clatm3.f(3)
π NAME π
TESTING/MATGEN/clatm3.f β the page title is TESTING/MATGEN/clatm3.f(3), and weβre in the Library Functions Manual. Welcome to the CLATM3 page, bestie. π«Ά
πΈ SYNOPSIS πΈ
Functions/Subroutines
complex function clatm3 (m, n, i, j, isub, jsub, kl, ku, idist, iseed,
d, igrade, dl, dr, ipvtng, iwork, sparse)
CLATM3
π§ Function/Subroutine Documentation π§
complex function clatm3 (integer m, integer n, integer i, integer j, integer
isub, integer jsub, integer kl, integer ku, integer idist, integer,
dimension( 4 ) iseed, complex, dimension( * ) d, integer igrade,
complex, dimension( * ) dl, complex, dimension( * ) dr, integer ipvtng,
integer, dimension( * ) iwork, real sparse)
CLATM3
β¨ Purpose β¨
okay so hereβs the tea: CLATM3 returns the (ISUB,JSUB) entry of a random matrix of dimension (M, N) described by the other parameters. The (ISUB,JSUB) spot is the final position of the (I,J) entry after pivoting according to IPVTNG and IWORK. CLATM3 gets called by CLATMR to build random test matrices. No error checking on parameters is done here, because this routine is called in a tight loop by CLATMR, which has already checked the parameters. no cap. π
The way CLATM3 works differs from CLATM2 in the order in which the random number generator is called to fill in random matrix entries. With CLATM2, the generator is called to fill in the pivoted matrix columnwise. With CLATM3, the generator is called to fill in the matrix columnwise, after which it is pivoted. Thus, CLATM3 can be used to construct random matrices which differ only in their order of rows and/or columns. CLATM2 is used to construct band matrices while avoiding calling the random number generator for entries outside the band (and therefore generating random numbers in different orders for different pivot orders). iconic behavior, honestly. π¦
The matrix whose (ISUB,JSUB) entry is returned is constructed as follows (this routine only computes one entry):
- If
ISUBis outside(1..M)orJSUBis outside(1..N), return zero (this is convenient for generating matrices in band format). - Generate a matrix
Awith random entries of distributionIDIST. - Set the diagonal to
D. - Grade the matrix, if desired, from the left (by
DL) and/or from the right (byDRorDL) as specified byIGRADE. - Permute, if desired, the rows and/or columns as specified by
IPVTNGandIWORK. - Band the matrix to have lower bandwidth
KLand upper bandwidthKU. - Set random entries to zero as specified by
SPARSE.
π Parameters π
MβINTEGER. Number of rows of matrix. Not modified.NβINTEGER. Number of columns of matrix. Not modified.IβINTEGER. Row of unpivoted entry to be returned. Not modified.JβINTEGER. Column of unpivoted entry to be returned. Not modified.ISUBβINTEGER. Row of pivoted entry to be returned. Changed on exit.JSUBβINTEGER. Column of pivoted entry to be returned. Changed on exit.KLβINTEGER. Lower bandwidth. Not modified.KUβINTEGER. Upper bandwidth. Not modified.IDISTβINTEGER. On entry,IDISTspecifies the type of distribution to be used to generate a random matrix. Not modified.1=> real and imaginary parts eachUNIFORM( 0, 1 )2=> real and imaginary parts eachUNIFORM( -1, 1 )3=> real and imaginary parts eachNORMAL( 0, 1 )4=> complex number uniform inDISK( 0 , 1 )
ISEEDβINTEGERarray of dimension( 4 ). Seed for random number generator. Changed on exit.DβCOMPLEXarray of dimension( MIN( I , J ) ). Diagonal entries of matrix. Not modified.IGRADEβINTEGER. Specifies grading of matrix as follows. Not modified.0=> no grading1=> matrix premultiplied bydiag( DL )2=> matrix postmultiplied bydiag( DR )3=> matrix premultiplied bydiag( DL )and postmultiplied bydiag( DR )4=> matrix premultiplied bydiag( DL )and postmultiplied byinv( diag( DL ) )5=> matrix premultiplied bydiag( DL )and postmultiplied bydiag( CONJG(DL) )6=> matrix premultiplied bydiag( DL )and postmultiplied bydiag( DL )
DLβCOMPLEXarray( I or J, as appropriate ). Left scale factors for grading matrix. Not modified.DRβCOMPLEXarray( I or J, as appropriate ). Right scale factors for grading matrix. Not modified.IPVTNGβINTEGER. On entry specifies pivoting permutations as follows. Not modified.0=> none.1=> row pivoting.2=> column pivoting.3=> full pivoting, i.e., on both sides.
IWORKβINTEGERarray( I or J, as appropriate ). This array specifies the permutation used. The row (or column) originally in positionKis in positionIWORK( K )after pivoting. This differs fromIWORKforCLATM2. Not modified.SPARSEβREALbetween0.and1.β on entry specifies the sparsity of the matrix if sparse matrix is to be generated.SPARSEshould lie between 0 and 1. A uniform( 0, 1 )random numberxis generated and compared toSPARSE; ifxis larger the matrix entry is unchanged and ifxis smaller the entry is set to zero. Thus on the average a fractionSPARSEof the entries will be set to zero. Not modified.
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 226 of file clatm3.f.
π Author π
This page was generated automatically by Doxygen for LAPACK from the source code.
Footer vibes, bestie: LAPACK Version 3.12.0, source page TESTING/MATGEN/clatm3.f(3) β no explicit date is listed in the footer, so we are absolutely not inventing one. periodt. π
β¨
π¦ where this page comes from
- arch rolling: lapack-doc 3.12.1-2 Β· custom
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms