PZDTTRSV(l)
β¨ NAME β¨
okay so bestie, meet PZDTTRSV(l) β the LAPACK routine that solves a tridiagonal triangular system of linear equations:
A(1:N, JA:JA+N-1) * X = B(IB:IB+N-1, 1:NRHS)
no cap, that's the whole name vibe. π
π¦ SYNOPSIS π¦
the lineup, exactly as it appears:
SUBROUTINE PZDTTRSV( UPLO, TRANS, N, NRHS, DL, D, DU, JA, DESCA, B, IB,
DESCB, AF, LAF, WORK, LWORK, INFO )
CHARACTER TRANS, UPLO
INTEGER IB, INFO, JA, LAF, LWORK, N, NRHS
INTEGER DESCA( * ), DESCB( * )
COMPLEX*16 AF( * ), B( * ), D( * ), DL( * ), DU( * ), WORK( * )
like, that's the whole cast. UPLO, TRANS, N, NRHS, DL, D, DU, JA, DESCA, B, IB, DESCB, AF, LAF, WORK, LWORK, INFO. iconic. π
πΈ PURPOSE πΈ
PZDTTRSV solves a tridiagonal triangular system of linear equations, and it's giving either:
A(1:N, JA:JA+N-1) * X = B(IB:IB+N-1, 1:NRHS)
or
A(1:N, JA:JA+N-1)^H * X = B(IB:IB+N-1, 1:NRHS)
where A(1:N, JA:JA+N-1) is a tridiagonal triangular matrix factor produced by the Gaussian elimination code PZ@(dom_pre)TTRF and is stored in A(1:N,JA:JA+N-1) and AF.
the matrix stored in A(1:N, JA:JA+N-1) is either upper or lower triangular according to UPLO, and the choice of solving A(1:N, JA:JA+N-1) or A(1:N, JA:JA+N-1)^H is dictated by the user by the parameter TRANS.
β οΈ routine PZDTTRF MUST be called first. like, no skipping. periodt. π
π that's the tea. this was LAPACK routine version 1.5, dated 12 May 1997, and the page is PZDTTRSV(l). bye bestie, go slay those tridiagonal systems π«Άβ¨
π¦ where this page comes from
- ubuntu jammy: scalapack-doc 1.5-11
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash on : a language model's rewrite of the original, under its terms