SRC/ssytri2.f(3)

hey bestie πŸ’– you just opened SRC/ssytri2.f(3) from the Library Functions Manual. buckle up, because this page is lowkey iconic and also very much doing linear algebra.

✨ NAME ✨

SRC/ssytri2.f

πŸ’… SYNOPSIS πŸ’…

🌸 Functions/Subroutines 🌸

subroutine ssytri2 (uplo, n, a, lda, ipiv, work, lwork, info)
    SSYTRI2

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

subroutine ssytri2 (character uplo, integer n, real, dimension( lda, * ) a,
    integer lda, integer, dimension( * ) ipiv, real, dimension( * ) work,
    integer lwork, integer info)
    SSYTRI2

πŸŽ€ Purpose πŸŽ€

okay so SSYTRI2 computes the inverse of a REAL symmetric indefinite matrix A using the factorization A = U*D*U**T or A = L*D*L**T computed by SSYTRF. SSYTRI2 sets the LEADING DIMENSION of the workspace before calling SSYTRI2X, which actually computes the inverse. no cap, that's the whole plot. ✨

🌈 Parameters 🌈

  • UPLO β€” UPLO is CHARACTER*1. Specifies whether the details of the factorization are stored as an upper or lower triangular matrix.
    • = 'U': Upper triangular, form is A = U*D*U**T;
    • = 'L': Lower triangular, form is A = L*D*L**T.
  • N β€” N is INTEGER. The order of the matrix A. N >= 0.
  • A β€” A is REAL array, dimension (LDA,N). On entry, the block diagonal matrix D and the multipliers used to obtain the factor U or L as computed by SSYTRF. On exit, if INFO = 0, the (symmetric) inverse of the original matrix. If UPLO = 'U', the upper triangular part of the inverse is formed and the part of A below the diagonal is not referenced; if UPLO = 'L' the lower triangular part of the inverse is formed and the part of A above the diagonal is not referenced.
  • LDA β€” LDA is INTEGER. The leading dimension of the array A. LDA >= max(1,N).
  • IPIV β€” IPIV is INTEGER array, dimension (N). Details of the interchanges and the block structure of D as determined by SSYTRF.
  • WORK β€” WORK is REAL array, dimension (MAX(1,LWORK)).
  • LWORK β€” LWORK is INTEGER. The dimension of the array WORK. If N = 0, LWORK >= 1, else LWORK >= (N+NB+1)*(NB+3). If LWORK = -1, then a workspace query is assumed; the routine calculates:
    • the optimal size of the WORK array, returns this value as the first entry of the WORK array,
    • and no error message related to LWORK is issued by XERBLA.
  • INFO β€” INFO is INTEGER.
    • = 0: successful exit
    • < 0: if INFO = -i, the i-th argument had an illegal value
    • > 0: if INFO = i, D(i,i) = 0; the matrix is singular and its inverse could not be computed.

🫢 Author 🫢

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

Definition at line 126 of file ssytri2.f.

πŸ‘‘ Author πŸ‘‘

Generated automatically by Doxygen for LAPACK from the source code.


signing off with footer energy, bestie: LAPACK, Version 3.12.0, page SRC/ssytri2.f(3). πŸ’œβœ¨
πŸ“¦ where this page comes from