/home/abui...N/dlatm3.f(3)
This is the dlatm3.f(3) page from the Library Functions Manual β yep, the path goes wild, but that's literally the vibe. Buckle up, bestie, it's a LAPACK test-matrix situation. π
πΈ NAME πΈ
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/TESTING/MATGEN/dlatm3.f
Like, that's the whole name section. She's long, she's iconic, she contains the entire build path. Periodt. π
β¨ SYNOPSIS β¨
π§ Functions/Subroutines π§
double precision function DLATM3 (m, n, i, j, isub, jsub, kl, ku, idist,
iseed, d, igrade, dl, dr, ipvtng, iwork, sparse)
DLATM3
π Function/Subroutine Documentation π
Okay so here's the full signature, bestie β not me copying it exactly because that's the whole point:
double precision function DLATM3 (integer m, integer n, integer i, integer
j, integer isub, integer jsub, integer kl, integer ku, integer idist,
integer, dimension( 4 ) iseed, double precision, dimension( * ) d,
integer igrade, double precision, dimension( * ) dl, double precision,
dimension( * ) dr, integer ipvtng, integer, dimension( * ) iwork, double
precision sparse)
DLATM3
Yep: DLATM3. That's her name. Main character energy fr fr. π
π¦ Purpose π¦
DLATM3returns the(ISUB,JSUB)entry of a random matrix of dimension(M, N)described by the other parameters.(ISUB,JSUB)is the final position of the(I,J)entry after pivoting according toIPVTNGandIWORK.DLATM3is called by theDLATMRroutine in order to build random test matrices. No error checking on parameters is done, because this routine is called in a tight loop byDLATMRwhich has already checked the parameters.
Lowkey important: the page says use of DLATM3 differs from SLATM2 in the order in which the random number generator is called to fill in random matrix entries. (Yes, it literally says SLATM2 there β keep that in your notes, bestie.) With DLATM2, the generator is called to fill in the pivoted matrix columnwise. With DLATM3, the generator is called to fill in the matrix columnwise, after which it is pivoted. Thus, DLATM3 can be used to construct random matrices which differ only in their order of rows and/or columns. DLATM2 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. π
The matrix whose (ISUB,JSUB) entry is returned is constructed as follows β and bestie, this routine only computes one entry, that's the whole assignment:
- 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.
That's the whole construction pipeline. Slay. π
π 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:1=>UNIFORM( 0, 1 )2=>UNIFORM( -1, 1 )3=>NORMAL( 0, 1 )Not modified. π
ISEEDβINTEGERarray of dimension( 4 ). Seed for random number generator. Changed on exit. π±DβDOUBLE PRECISIONarray of dimension( MIN( I , J ) ). Diagonal entries of matrix. Not modified. πIGRADEβINTEGER. Specifies grading of matrix as follows: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( DL )Not modified. β¨
DLβDOUBLE PRECISIONarray( I or J, as appropriate ). Left scale factors for grading matrix. Not modified. π¦DRβDOUBLE PRECISIONarray( I or J, as appropriate ). Right scale factors for grading matrix. Not modified. π¦IPVTNGβINTEGER. On entry specifies pivoting permutations as follows:0=> none.1=> row pivoting.2=> column pivoting.3=> full pivoting, i.e., on both sides. Not modified. π
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 fromIWORKforDLATM2. Not modified. πSPARSEβDOUBLE PRECISIONbetween0.and1.. On entry specifies the sparsity of the matrix if sparse matrix is to be generated.SPARSEshould lie between0and1. 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 223 of file dlatm3.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
Mwah π that's the page. Footer flavor for the record: LAPACK, Version 3.12.0, page path /home/abui...N/dlatm3.f(3) β and there's no extra date listed on this one, so don't come for me, bestie. π«Άβ¨
π¦ where this page comes from
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms