SRC/dtbcon.f(3)

πŸ’… SRC/dtbcon.f(3) πŸ’…

welcome to the Library Functions Manual corner of the group chat, bestie 🫢 this page is all about the Fortran subroutine dtbcon, and no cap it is giving serious condition-number-detective energy.


🌸 NAME 🌸

SRC/dtbcon.f


✨ SYNOPSIS ✨

πŸŽ€ Functions/Subroutines πŸŽ€

       subroutine dtbcon (norm, uplo, diag, n, kd, ab, ldab, rcond, work,
           iwork, info)
           DTBCON

πŸ¦‹ Function/Subroutine Documentation πŸ¦‹

okay so here is the full signature, copied exactly because we do not mess with code πŸ’…

   subroutine dtbcon (character norm, character uplo, character diag, integer
       n, integer kd, double precision, dimension( ldab, * ) ab, integer ldab,
       double precision rcond, double precision, dimension( * ) work, integer,
       dimension( * ) iwork, integer info)
       DTBCON

πŸ”₯ Purpose πŸ”₯

DTBCON estimates the reciprocal of the condition number of a triangular band matrix A, in either the 1-norm or the infinity-norm. periodt πŸ’œ

The norm of A is computed and an estimate is obtained for norm(inv(A)), then the reciprocal of the condition number is computed as

RCOND = 1 / ( norm(A) * norm(inv(A)) ).

literally the whole vibe: get norm(A), estimate norm(inv(A)), flip it into RCOND ✨

πŸ“ Parameters πŸ“

  • NORM β€” CHARACTER*1. Specifies whether the 1-norm condition number or the infinity-norm condition number is required:
    • = '1' or 'O': 1-norm;
    • = 'I': Infinity-norm.
  • UPLO β€” CHARACTER*1.
    • = 'U': A is upper triangular;
    • = 'L': A is lower triangular.
  • DIAG β€” CHARACTER*1.
    • = 'N': A is non-unit triangular;
    • = 'U': A is unit triangular.
  • N β€” INTEGER. The order of the matrix A. N >= 0.
  • KD β€” INTEGER. The number of superdiagonals or subdiagonals of the triangular band matrix A. KD >= 0.
  • AB β€” DOUBLE PRECISION array, dimension (LDAB,N). The upper or lower triangular band matrix A, stored in the first kd+1 rows of the array. The j-th column of A is stored in the j-th column of the array AB as follows:
if UPLO = 'U', AB(kd+1+i-j,j) = A(i,j) for max(1,j-kd)<=i<=j;
if UPLO = 'L', AB(1+i-j,j)    = A(i,j) for j<=i<=min(n,j+kd).

If DIAG = 'U', the diagonal elements of A are not referenced and are assumed to be 1.

  • LDAB β€” INTEGER. The leading dimension of the array AB. LDAB >= KD+1.
  • RCOND β€” DOUBLE PRECISION. The reciprocal of the condition number of the matrix A, computed as RCOND = 1/(norm(A) * norm(inv(A))).
  • WORK β€” DOUBLE PRECISION array, dimension (3*N).
  • IWORK β€” INTEGER array, dimension (N).
  • INFO β€” 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 141 of file dtbcon.f.


🫢 Author 🫢

Generated automatically by Doxygen for LAPACK from the source code.


okay bye bestie πŸ’œ signing off from LAPACK, Version 3.12.0, page SRC/dtbcon.f(3) β€” slay and stay safe out there πŸŽ€

πŸ“¦ where this page comes from