TESTING/LIN/cpot03.f(3)

πŸ“¦
version 3.12.1 Β· 3 versions
, in arch rolling

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 is CHARACTER*1, and it tells you which triangular part of the Hermitian matrix A is stored:
    • = 'U': Upper triangular
    • = 'L': Lower triangular
  • N β€” this is INTEGER, and it's the number of rows and columns of the matrix A. Constraint: N >= 0.
  • A β€” this is COMPLEX array, dimension (LDA,N). It's the original Hermitian matrix A, bestie.
  • LDA β€” this is INTEGER, and it's the leading dimension of the array A. Constraint: LDA >= max(1,N).
  • AINV β€” this is COMPLEX array, dimension (LDAINV,N). On entry, it's the inverse of the matrix A, stored as a Hermitian matrix in the same format as A. In this version, AINV is expanded into a full matrix and multiplied by A, so the opposing triangle of AINV will be changed β€” like, if the upper triangular part of AINV is stored, the lower triangular part gets used as work space. 🫢
  • LDAINV β€” this is INTEGER, and it's the leading dimension of the array AINV. Constraint: LDAINV >= max(1,N).
  • WORK β€” this is COMPLEX array, dimension (LDWORK,N).
  • LDWORK β€” this is INTEGER, and it's the leading dimension of the array WORK. Constraint: LDWORK >= max(1,N).
  • RWORK β€” this is REAL array, dimension (N).
  • RCOND β€” this is REAL, and it's the reciprocal of the condition number of A, computed as ( 1/norm(A) ) / norm(AINV).
  • RESID β€” this is REAL, and it's norm(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