TESTING/LIN/cpot03.f(3)
hey bestie, welcome to TESTING/LIN/cpot03.f(3) β that's the page name β from the Library Functions Manual, aka the LAPACK corner of the internet. πβ¨
πΈ NAME πΈ
TESTING/LIN/cpot03.f
π¦ SYNOPSIS π¦
β¨ Functions/Subroutines β¨
Functions/Subroutines
subroutine cpot03 (uplo, n, a, lda, ainv, ldainv, work, ldwork, rwork,
rcond, resid)
CPOT03
π Function/Subroutine Documentation π
subroutine cpot03 (character uplo, integer n, complex, dimension( lda, * )
a, integer lda, complex, dimension( ldainv, * ) ainv, integer ldainv,
complex, dimension( ldwork, * ) work, integer ldwork, real, dimension( *
) rwork, real rcond, real resid)
CPOT03
π Purpose π
The code comment says it best, and it says:
!>
!> CPOT03 computes the residual for a Hermitian matrix times its
!> inverse:
!> norm( I - A*AINV ) / ( N * norm(A) * norm(AINV) * EPS ),
!> where EPS is the machine epsilon.
!>
So yes, CPOT03 is out here computing the residual for a Hermitian matrix times its inverse. That formula is the whole vibe check, no cap. π₯
π Parameters π
UPLOβ this isCHARACTER*1, and it tells you which triangular part of the Hermitian matrixAis stored:= 'U': Upper triangular= 'L': Lower triangular
Nβ this isINTEGER, and it's the number of rows and columns of the matrixA. Constraint:N >= 0.Aβ this isCOMPLEX array, dimension (LDA,N). It's the original Hermitian matrixA, bestie.LDAβ this isINTEGER, and it's the leading dimension of the arrayA. Constraint:LDA >= max(1,N).AINVβ this isCOMPLEX array, dimension (LDAINV,N). On entry, it's the inverse of the matrixA, stored as a Hermitian matrix in the same format asA. In this version,AINVis expanded into a full matrix and multiplied byA, so the opposing triangle ofAINVwill be changed β like, if the upper triangular part ofAINVis stored, the lower triangular part gets used as work space. π«ΆLDAINVβ this isINTEGER, and it's the leading dimension of the arrayAINV. Constraint:LDAINV >= max(1,N).WORKβ this isCOMPLEX array, dimension (LDWORK,N).LDWORKβ this isINTEGER, and it's the leading dimension of the arrayWORK. Constraint:LDWORK >= max(1,N).RWORKβ this isREAL array, dimension (N).RCONDβ this isREAL, and it's the reciprocal of the condition number ofA, computed as( 1/norm(A) ) / norm(AINV).RESIDβ this isREAL, and it'snorm(I - A*AINV) / ( N * norm(A) * norm(AINV) * EPS ).
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π Definition at line 124 of file cpot03.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
okay bye bestie β this page is from LAPACK Version 3.12.0, and the footer's page name is TESTING/LIN/cpot03.f(3). Those are the source/version deets, periodt. π
β¨
π¦ 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