SRC/cgtcon.f(3)

Bestie, welcome to SRC/cgtcon.f(3) β€” yes, that’s the header title β€” and it lives in the Library Functions Manual. Buckle up, because this page is lowkey a whole vibe: it’s the tea on CGTCON, and it’s giving numerical linear algebra excellence. πŸ’…βœ¨

πŸ’– NAME πŸ’–

SRC/cgtcon.f

🌈 SYNOPSIS 🌈

🧩 Functions/Subroutines 🧩

  • subroutine cgtcon (norm, n, dl, d, du, du2, ipiv, anorm, rcond, work, info) β€” CGTCON

πŸ› οΈ Function/Subroutine Documentation πŸ› οΈ

subroutine cgtcon (character norm, integer n, complex, dimension( * ) dl, complex, dimension( * ) d, complex, dimension( * ) du, complex, dimension( * ) du2, integer, dimension( * ) ipiv, real anorm, real rcond, complex, dimension( * ) work, integer info)
CGTCON

🎯 Purpose 🎯

Okay so CGTCON estimates the reciprocal of the condition number of a complex tridiagonal matrix A using the LU factorization as computed by CGTTRF. ✨

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 main-character formula, no cap. πŸ’–

πŸŽ›οΈ Parameters πŸŽ›οΈ

  • NORM β€” NORM is CHARACTER*1. It specifies whether the 1-norm condition number or the infinity-norm condition number is required:
    • = '1' or 'O': 1-norm;
    • = 'I': Infinity-norm.
  • N β€” N is INTEGER. This is the order of the matrix A. N >= 0.
  • DL β€” DL is COMPLEX array, dimension (N-1). These are the (n-1) multipliers that define the matrix L from the LU factorization of A as computed by CGTTRF.
  • D β€” D is COMPLEX array, dimension (N). These are the n diagonal elements of the upper triangular matrix U from the LU factorization of A.
  • DU β€” DU is COMPLEX array, dimension (N-1). These are the (n-1) elements of the first superdiagonal of U.
  • DU2 β€” DU2 is COMPLEX array, dimension (N-2). These are the (n-2) elements of the second superdiagonal of U.
  • IPIV β€” IPIV is INTEGER array, dimension (N). These are the pivot indices; for 1 <= i <= n, row i of the matrix was interchanged with row IPIV(i). IPIV(i) will always be either i or i+1; IPIV(i) = i indicates a row interchange was not required.
  • ANORM β€” ANORM is REAL. If NORM = '1' or 'O', this is the 1-norm of the original matrix A. If NORM = 'I', this is the infinity-norm of the original matrix A.
  • RCOND β€” RCOND is REAL. This is the reciprocal of the condition number of the matrix A, computed as RCOND = 1/(ANORM * AINVNM), where AINVNM is an estimate of the 1-norm of inv(A) computed in this routine.
  • WORK β€” WORK is COMPLEX array, dimension (2*N).
  • INFO β€” INFO is INTEGER:
    • = 0: successful exit
    • < 0: if INFO = -i, the i-th argument had an illegal value

πŸ‘‘ Author πŸ‘‘

  • Univ. of Tennessee
  • Univ. of California Berkeley
  • Univ. of Colorado Denver
  • NAG Ltd.

Definition at line 139 of file cgtcon.f.

πŸ‘‘ Author πŸ‘‘

Generated automatically by Doxygen for LAPACK from the source code.


Stay reciprocal and radiant, bestie πŸ«ΆπŸ’– β€” the footer is giving LAPACK Version 3.12.0 and SRC/cgtcon.f(3), with the date not shown on this page. Periodt. πŸ’…βœ¨

πŸ“¦ where this page comes from