DSPCON(3)
π DSPCON(3) β like, the Library Functions Manual, but make it cute π
okay so bestie, welcome to the Library Functions Manual entry for DSPCON(3), straight out of the LAPACK-verse β¨ the file lives at /home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/SRC/dspcon.f and honestly? it's giving "quietly doing the most important math in the room" π
β¨ SYNOPSIS β¨
Functions/Subroutines β the cast π
subroutine DSPCON (uplo, n, ap, ipiv, anorm, rcond, work, iwork, info)
DSPCON β that's her, that's the main character π
Function/Subroutine Documentation π
subroutine DSPCON (character uplo, integer n, double precision, dimension( * )
ap, integer, dimension( * ) ipiv, double precision anorm, double
precision rcond, double precision, dimension( * ) work, integer,
dimension( * ) iwork, integer info)
yep, DSPCON again β she said "billing twice, I earned it" π
πΈ PURPOSE πΈ
so what does she do? DSPCON estimates the reciprocal of the condition number (in the 1-norm) of a real symmetric packed matrix A, using the factorization A = U*D*U**T or A = L*D*L**T computed by DSPTRF π«Ά
an estimate is obtained for norm(inv(A)), and the reciprocal of the condition number is computed as RCOND = 1 / (ANORM * norm(inv(A))) β that's the whole formula, no cap, it's literally right there πβ¨
π¦ PARAMETERS π¦
every single one of them, bestie, no skips π
UPLOβUPLO is CHARACTER*1. Specifies whether the details of the factorization are stored as an upper or lower triangular matrix ποΈ= 'U': Upper triangular, form isA = U*D*U**T;= 'L': Lower triangular, form isA = L*D*L**T.
NβN is INTEGER. The order of the matrixA. Constraint, and she's strict about it:N >= 0π¨APβAP is DOUBLE PRECISION array, dimension (N*(N+1)/2). The block diagonal matrixDand the multipliers used to obtain the factorUorLas computed byDSPTRF, stored as a packed triangular matrix π§© (packed!! like a little suitcase of math)IPIVβIPIV is INTEGER array, dimension (N). Details of the interchanges and the block structure ofDas determined byDSPTRFπANORMβANORM is DOUBLE PRECISION. The 1-norm of the original matrixAπRCONDβRCOND is DOUBLE PRECISION. The reciprocal of the condition number of the matrixA, computed asRCOND = 1/(ANORM * AINVNM), whereAINVNMis an estimate of the 1-norm ofinv(A)computed in this routine π«WORKβWORK is DOUBLE PRECISION array, dimension (2*N)π οΈIWORKβIWORK is INTEGER array, dimension (N)π’INFOβINFO is INTEGERπ£= 0: successful exit π< 0: ifINFO = -i, the i-th argument had an illegal value π (so likeβ¦ check your inputs, bestie)
β οΈ heads up: that INFO < 0 situation means you passed something illegal into the i-th argument. it's not a "vibe", it's an error. respect the constraint π
π AUTHOR π
iconic institutions only, we stan credentialed excellence:
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π Definition at line 123 of file dspcon.f.
π AUTHOR (the meta one) π
Generated automatically by Doxygen for LAPACK from the source code π€β¨ (yes bestie, even the docs have docs)
and that's the whole page!! π this lovely little document is LAPACK, Version 3.12.0, and it signs off from /home/abui...RC/dspcon.f(3) β mwah, drink water, and may your RCOND always be reasonable ππ
π¦ 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