SRC/ssycon_rook.f(3)
πΈ SRC/ssycon_rook.f(3) β from the Library Functions Manual π
okay so bestie, welcome to the glow-up of a genuinely iconic LAPACK page π this is SRC/ssycon_rook.f and she is that girl β a Fortran file that estimates condition numbers and does it flawlessly. let's get into it β¨
β¨ NAME β¨
SRC/ssycon_rook.f β literally the file, that's her name, no notes π
π¦ SYNOPSIS π¦
Functions/Subroutines
subroutine ssycon_rook (uplo, n, a, lda, ipiv, anorm, rcond, work, iwork, info)β this is the main character right here πSSYCON_ROOKβ the label, the name, the brand β¨
π Function/Subroutine Documentation π
subroutine ssycon_rook (character uplo, integer n, real, dimension( lda, * )
a, integer lda, integer, dimension( * ) ipiv, real anorm, real rcond,
real, dimension( * ) work, integer, dimension( * ) iwork, integer info)
SSYCON_ROOK β say it with your chest π
π― Purpose
!> !> SSYCON_ROOK estimates the reciprocal of the condition number (in the !> 1-norm) of a real symmetric matrix A using the factorization !> A = U*D*U**T or A = L*D*L**T computed by SSYTRF_ROOK. !> !> An estimate is obtained for norm(inv(A)), and the reciprocal of the !> condition number is computed as RCOND = 1 / (ANORM * norm(inv(A))). !>
so like. the whole vibe is: you already ran SSYTRF_ROOK, you got your factorization A = U*D*U**T or A = L*D*L**T, and now this routine slides in to give you RCOND β the reciprocal of the condition number in the 1-norm π§ it estimates norm(inv(A)) and then hits you with RCOND = 1 / (ANORM * norm(inv(A))). math but make it cute ππ
π Parameters
UPLOβ isCHARACTER*1. tells us whether the factorization details are chilling as an upper or lower triangular matrix, bestie:= 'U': upper triangular, form isA = U*D*U**T;= 'L': lower triangular, form isA = L*D*L**T.
Nβ isINTEGER. the order of the matrixA.N >= 0. periodt β¨Aβ isREAL array, dimension (LDA,N). the block diagonal matrixDplus the multipliers used to obtain the factorUorL, as computed bySSYTRF_ROOK. iconic π₯LDAβ isINTEGER. the leading dimension of the arrayA.LDA >= max(1,N). don't skip this one πIPIVβ isINTEGER array, dimension (N). all the tea about the interchanges and the block structure ofD, as determined bySSYTRF_ROOKπANORMβ isREAL. the 1-norm of the original matrixAπRCONDβ isREAL. 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. she's the output, she's the moment βWORKβ isREAL array, dimension (2*N). workspace, doing her thing πIWORKβ isINTEGER array, dimension (N). more workspace, integer flavored πINFOβ isINTEGER:= 0: successful exit, we love a win πΈ< 0: ifINFO = -i, the i-th argument had an illegal value β bestie someone passed a bad argument, call it out π
π Author
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
the whole squad showed up π«Ά
π Contributors:
!> !> December 2016, Igor Kozachenko, !> Computer Science Division, !> University of California, Berkeley !> !> September 2007, Sven Hammarling, Nicholas J. Higham, Craig Lucas, !> School of Mathematics, !> University of Manchester !> !>
shoutout to everyone who built her ππ
π Definition at line 142 of file ssycon_rook.f.
go look, she's literally right there π
π Author π
Generated automatically by Doxygen for LAPACK from the source code. the robots did this one π€β¨
and that's the whole page bestie, no cap, every flag and every line accounted for π
this cutie is from LAPACK, Version 3.12.0 and her file is SRC/ssycon_rook.f(3) β go estimate some condition numbers, you're gonna slay π¦ππ
π¦ where this page comes from
- arch rolling: lapack-doc 3.12.1-2 Β· custom
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms