DPBEQU(3)
π DPBEQU(3) β Library Functions Manual π
okay so bestie, welcome to DPBEQU(3) from the Library Functions Manual β this is the LAPACK routine page for DPBEQU, and we are gonna get into the scaling tea ββ¨
πΈ NAME πΈ
the file for this routine lives here, no cap:
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/SRC/dpbequ.f
β¨ SYNOPSIS β¨
π¦ Functions/Subroutines π¦
subroutine DPBEQU (uplo, n, kd, ab, ldab, s, scond, amax, info)
DPBEQU
π Function/Subroutine Documentation π
the full signature is giving very much βI know my typesβ energy:
subroutine DPBEQU (character uplo, integer n, integer kd, double precision,
dimension( ldab, * ) ab, integer ldab, double precision, dimension( * )
s, double precision scond, double precision amax, integer info)
DPBEQU
πΈ Purpose πΈ
DPBEQU computes row and column scalings intended to equilibrate a symmetric 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. periodt π
π Parameters π
UPLOβCHARACTER*1= 'U': Upper triangular ofAis stored;= 'L': Lower triangular ofAis stored.
NβINTEGER; the order of the matrixA.N >= 0.KDβINTEGER; the number of superdiagonals of the matrixAifUPLO = 'U', or the number of subdiagonals ifUPLO = 'L'.KD >= 0.ABβDOUBLE PRECISION array, dimension (LDAB,N); the upper or lower triangle of the symmetric band matrixA, stored in the firstKD+1rows of the array. Thej-th column ofAis stored in thej-th column of the arrayABas follows:
if UPLO = 'U', AB(kd+1+i-j,j) = A(i,j) for max(1,j-kd)<=i<=j;
if UPLO = 'L', AB(1+i-j,j) = A(i,j) for j<=i<=min(n,j+kd).
LDABβINTEGER; the leading dimension of the arrayA.LDAB >= KD+1.SβDOUBLE PRECISION array, dimension (N); ifINFO = 0,Scontains the scale factors forA.SCONDβDOUBLE PRECISION; ifINFO = 0, the page literally saysScontains 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βDOUBLE PRECISION; absolute value of largest matrix element. IfAMAXis very close to overflow or very close to underflow, the matrix should be scaled.INFOβINTEGER= 0: successful exit< 0: ifINFO = -i, thei-th argument had an illegal value.> 0: ifINFO = i, thei-th diagonal element is nonpositive.
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 128 of file dpbequ.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
okay thatβs the whole vibe check π«Ά this footer moment is brought to you by LAPACK, Version 3.12.0, with the path /home/abui...RC/dpbequ.f(3) β byeee, ily, go scale responsibly π
β¨
π¦ 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