complexGBauxiliary(3)
π complexGBauxiliary β bestie, it's the GB auxiliary girlies π
okay so real talk, this is the LAPACK page for complexGBauxiliary, aka the little support squad that hangs out behind the scenes of all your general band matrix (GB) operations. it's giving "backstage crew who are actually the main characters" energy, fr fr. β¨
πΈ NAME πΈ
complexGBauxiliary β bestie, it's the complex auxiliary function family, and it's here to help. π
π¦ SYNOPSIS π¦
β¨ Functions β¨
real function clangb (NORM, N, KL, KU, AB, LDAB, WORK)βCLANGBreturns the value of the 1-norm, Frobenius norm, infinity-norm, or the largest absolute value of any element of a general band matrix. she's a norm girlie and she knows it πsubroutine claqgb (M, N, KL, KU, AB, LDAB, R, C, ROWCND, COLCND, AMAX, EQUED)βCLAQGBscales a general band matrix, using row and column scaling factors computed bysgbequ. so basically it does the glow-up pass πΈ
π Detailed Description π
This is the group of complex auxiliary functions for GB matrices. that's it, that's the vibe. they exist to support the band matrix operations and we love them for it. π«Ά
π₯ Function Documentation π₯
β real function clangb (character NORM, integer N, integer KL, integer KU, complex, dimension( ldab, * ) AB, integer LDAB, real, dimension( * ) WORK)
CLANGB returns the value of the 1-norm, Frobenius norm, infinity-norm, or the largest absolute value of any element of general band matrix.
π Purpose:
CLANGB returns the value of the one norm, or the Frobenius norm, or
the infinity norm, or the element of largest absolute value of an
n by n band matrix A, with kl sub-diagonals and ku super-diagonals.
π Returns
CLANGB
CLANGB = ( max(abs(A(i,j))), NORM = 'M' or 'm'
(
( norm1(A), NORM = '1', 'O' or 'o'
(
( normI(A), NORM = 'I' or 'i'
(
( normF(A), NORM = 'F', 'f', 'E' or 'e'
where norm1 denotes the one norm of a matrix (maximum column sum), normI denotes the infinity norm of a matrix (maximum row sum) and normF denotes the Frobenius norm of a matrix (square root of sum of squares). Note that max(abs(A(i,j))) is not a consistent matrix norm. β οΈ so like, keep that in mind bestie.
π§ Parameters
NORMβNORMisCHARACTER*1. Specifies the value to be returned inCLANGBas described above. πNβNisINTEGER. The order of the matrix A.N >= 0. WhenN = 0,CLANGBis set to zero. (she knows when to sit quietly and return nothing, iconic) πKLβKLisINTEGER. The number of sub-diagonals of the matrix A.KL >= 0.KUβKUisINTEGER. The number of super-diagonals of the matrix A.KU >= 0.ABβABisCOMPLEX array, dimension (LDAB,N). The band matrix A, stored in rows1toKL+KU+1. The j-th column of A is stored in the j-th column of the arrayABas follows:
AB(ku+1+i-j,j) = A(i,j) for max(1,j-ku)<=i<=min(n,j+kl).
LDABβLDABisINTEGER. The leading dimension of the arrayAB.LDAB >= KL+KU+1.WORKβWORKisREAL array, dimension (MAX(1,LWORK)), whereLWORK >= NwhenNORM = 'I'; otherwise,WORKis not referenced. lowkey she only shows up when you need her π
π Author
Univ. of Tennessee
Univ. of California Berkeley
Univ. of Colorado Denver
NAG Ltd.
π subroutine claqgb (integer M, integer N, integer KL, integer KU, complex, dimension( ldab, * ) AB, integer LDAB, real, dimension( * ) R, real, dimension( * ) C, real ROWCND, real COLCND, real AMAX, character EQUED)
CLAQGB scales a general band matrix, using row and column scaling factors computed by sgbequ. it's the equilibration makeover episode β¨
π¦ Purpose:
CLAQGB equilibrates a general M by N band matrix A with KL
subdiagonals and KU superdiagonals using the row and scaling factors
in the vectors R and C.
πΈ Parameters
MβMisINTEGER. The number of rows of the matrix A.M >= 0.NβNisINTEGER. The number of columns of the matrix A.N >= 0.KLβKLisINTEGER. The number of subdiagonals within the band of A.KL >= 0.KUβKUisINTEGER. The number of superdiagonals within the band of A.KU >= 0.ABβABisCOMPLEX array, dimension (LDAB,N).
On entry, the matrix A in band storage, in rows 1 to KL+KU+1. The j-th column of A is stored in the j-th column of the array AB as follows:
AB(ku+1+i-j,j) = A(i,j) for max(1,j-ku)<=i<=min(m,j+kl)
On exit, the equilibrated matrix, in the same storage format as A. See EQUED for the form of the equilibrated matrix. π«
LDABβLDABisINTEGER. The leading dimension of the arrayAB.LDA >= KL+KU+1. (yes bestie, it saysLDAhere, we don't make the rules, we just report them faithfully π)RβRisREAL array, dimension (M). The row scale factors for A.CβCisREAL array, dimension (N). The column scale factors for A.ROWCNDβROWCNDisREAL. Ratio of the smallestR(i)to the largestR(i).COLCNDβCOLCNDisREAL. Ratio of the smallestC(i)to the largestC(i).AMAXβAMAXisREAL. Absolute value of largest matrix entry.EQUEDβEQUEDisCHARACTER*1. Specifies the form of equilibration that was done.= 'N': No equilibration= 'R': Row equilibration, i.e., A has been premultiplied bydiag(R).= 'C': Column equilibration, i.e., A has been postmultiplied bydiag(C).= 'B': Both row and column equilibration, i.e., A has been replaced bydiag(R) * A * diag(C).
π§ Internal Parameters:
THRESH is a threshold value used to decide if row or column scaling
should be done based on the ratio of the row or column scaling
factors. If ROWCND < THRESH, row scaling is done, and if
COLCND < THRESH, column scaling is done.
LARGE and SMALL are threshold values used to decide if row scaling
should be done based on the absolute size of the largest matrix
element. If AMAX > LARGE or AMAX < SMALL, row scaling is done.
π 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. (not a human, just a very diligent little robot doing her best π«Ά)
and that's the whole page bestie!! mwah π this was Version 3.11.0, dated Sun Nov 27 2022, and honestly? she ate. 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