ZPBEQU(3)
Okay bestie, welcome to ZPBEQU(3) β from the Library Functions Manual β where we compute scalings for a Hermitian positive definite band matrix and make it behave. No cap, sheβs giving linear algebra main character energy. π
π NAME π
The source file for this routine is literally:
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/SRC/zpbequ.f
β¨ SYNOPSIS β¨
Functions/Subroutines
subroutine ZPBEQU (uplo, n, kd, ab, ldab, s, scond, amax, info)
ZPBEQU
π«Ά Function/Subroutine Documentation π«Ά
Full signature, character-for-character:
subroutine ZPBEQU (character uplo, integer n, integer kd, complex*16,
dimension( ldab, * ) ab, integer ldab, double precision, dimension( * )
s, double precision scond, double precision amax, integer info)
ZPBEQU
π― Purpose π―
ZPBEQU computes row and column scalings intended to equilibrate a Hermitian positive definite band matrix A and reduce its condition number (with respect to the two-norm). S contains the scale factors, S(i) = 1/sqrt(A(i,i)), chosen so that the scaled matrix B with elements B(i,j) = S(i)*A(i,j)*S(j) has ones on the diagonal. This choice of S puts the condition number of B within a factor N of the smallest possible condition number over all possible diagonal scalings. Obsessed with a diagonal of ones, honestly. β¨
π Parameters π
UPLOβUPLOisCHARACTER*1.= 'U':Upper triangular ofAis stored;= 'L':Lower triangular ofAis stored.
NβNisINTEGER. The order of the matrixA.N >= 0.KDβKDisINTEGER. The number of superdiagonals of the matrixAifUPLO = 'U', or the number of subdiagonals ifUPLO = 'L'.KD >= 0.ABβABisCOMPLEX*16array, dimension(LDAB,N). The upper or lower triangle of the Hermitian band matrixA, stored in the firstKD+1rows of the array. The j-th column ofAis stored in the j-th column of the arrayABas follows:- if
UPLO = 'U',AB(kd+1+i-j,j) = A(i,j)formax(1,j-kd)<=i<=j; - if
UPLO = 'L',AB(1+i-j,j) = A(i,j)forj<=i<=min(n,j+kd).
- if
LDABβLDABisINTEGER. The leading dimension of the arrayA.LDAB >= KD+1.SβSisDOUBLE PRECISIONarray, dimension(N). IfINFO = 0,Scontains the scale factors forA.SCONDβSCONDisDOUBLE PRECISION. IfINFO = 0,Scontains the ratio of the smallestS(i)to the largestS(i). IfSCOND >= 0.1andAMAXis neither too large nor too small, it is not worth scaling byS.AMAXβAMAXisDOUBLE PRECISION. Absolute value of largest matrix element. IfAMAXis very close to overflow or very close to underflow, the matrix should be scaled.INFOβINFOisINTEGER.= 0:successful exit< 0:ifINFO = -i, the i-th argument had an illegal value.> 0:ifINFO = i, the i-th diagonal element is nonpositive.
π Author π
Respect, respect:
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 129 of file zpbequ.f.
πΈ Author πΈ
Generated automatically by Doxygen for LAPACK from the source code.
Catch you later, bestie β this page is brought to you by LAPACK Version 3.12.0, with footer path /home/abui...RC/zpbequ.f(3). Slay responsibly. Periodt. π
π¦ 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