DPPTRI(3)

πŸ“¦ version 3.12.0, in opensuse-leap 16.0

πŸŽ€ NAME πŸŽ€

/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/SRC/dpptri.f

okay so this is literally that file bestie, and it lives in the Library Functions Manual (that's the manual name, folder it under that πŸ’…). It's the home of DPPTRI and we're about to get into it ✨


✨ SYNOPSIS ✨

Functions/Subroutines

subroutine DPPTRI (uplo, n, ap, info)

DPPTRI β€” that's her, that's the icon πŸ‘‘


πŸ’– Function/Subroutine Documentation πŸ’–

subroutine DPPTRI (character uplo, integer n, double precision, dimension( * ) ap, integer info)

DPPTRI

Purpose:

!>
!> DPPTRI computes the inverse of a real symmetric positive definite
!> matrix A using the Cholesky factorization A = U**T*U or A = L*L**T
!> computed by DPPTRF.
!>

no cap, she takes your already-Cholesky'd matrix (the one DPPTRF handed you) and flips it into an inverse. main character behavior fr πŸ”₯

Parameters

  • UPLO β€” it's a CHARACTER*1, so it's just one letter and you gotta pick a side:
!>          UPLO is CHARACTER*1
!>          = 'U':  Upper triangular factor is stored in AP;
!>          = 'L':  Lower triangular factor is stored in AP.
!>
  • N β€” an INTEGER, the order of the matrix A. and bestie, N >= 0, no negative nonsense please πŸ™ˆ
!>          N is INTEGER
!>          The order of the matrix A.  N >= 0.
!>
  • AP β€” DOUBLE PRECISION array, dimension (N*(N+1)/2). she's packed, she's columnwise, she's a whole linear array doing the most πŸ’…
!>          AP is DOUBLE PRECISION array, dimension (N*(N+1)/2)
!>          On entry, the triangular factor U or L from the Cholesky
!>          factorization A = U**T*U or A = L*L**T, packed columnwise as
!>          a linear array.  The j-th column of U or L is stored in the
!>          array AP as follows:
!>          if UPLO = 'U', AP(i + (j-1)*j/2) = U(i,j) for 1<=i<=j;
!>          if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) = L(i,j) for j<=i<=n.
!>
!>          On exit, the upper or lower triangle of the (symmetric)
!>          inverse of A, overwriting the input factor U or L.
!>
⚠️ heads up: on exit AP gets overwritten with the upper or lower triangle of the (symmetric) inverse of A. your original factor U or L? gone. she said out with the old 🫢
  • INFO β€” an INTEGER, the vibe check on how it went:
!>          INFO is INTEGER
!>          = 0:  successful exit
!>          < 0:  if INFO = -i, the i-th argument had an illegal value
!>          > 0:  if INFO = i, the (i,i) element of the factor U or L is
!>                zero, and the inverse could not be computed.
!>

like: 0 = slay, < 0 = argument i was illegal (not me...), > 0 = the (i,i) element of the factor U or L is zero so the inverse literally could not be computed 😭 that's the "cannot divide by zero" of it all.


πŸ‘‘ Author πŸ‘‘

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

🌸 Definition at line 92 of file dpptri.f. 🌸


πŸ¦‹ Author πŸ¦‹

Generated automatically by Doxygen for LAPACK from the source code. (yes bestie she was written by a robot, and she ate πŸ’€βœ¨)


okay so that's the whole page β€” LAPACK, Version 3.12.0, dated /home/abui...RC/dpptri.f(3) in the header and <rainbow>Version 3.12.0 /home/abui...RC/dpptri.f(3)</rainbow> down at the bottom in the footer πŸŽ€ go forth and invert, ily πŸ«ΆπŸ’–

πŸ“¦ where this page comes from