SPTTRSV(l)
β¨ SPTTRSV(l) β¨
okay bestie sit DOWN because we're about to talk about SPTTRSV β a LAPACK routine, specifically from LAPACK version 2.0 π this is the SPTTRSV(l) manual page and she's giving linear algebra main character energy, no cap.
πΈ NAME πΈ
SPTTRSV β solve one of the triangular systems L**T* X = B, or L * X = B,
(yes the comma is staying, we don't edit her π )
π¦ SYNOPSIS π¦
SUBROUTINE SPTTRSV( TRANS, N, NRHS, D, E, B, LDB, INFO )
CHARACTER TRANS
INTEGER INFO, LDB, N, NRHS
REAL D( * )
REAL B( LDB, * ), E( * )
π PURPOSE π
so like, what does she even DO? SPTTRSV solves one of the triangular systems
L**T* X = B, or L * X = B,
where L is the Cholesky factor of a Hermitian positive definite tridiagonal matrix A such that
A = L*D*L**H
and that factorization is computed by SPTTRF, just so you know where she came from π«Ά
π ARGUMENTS π
TRANS(input)CHARACTERβ tells us which form of the system of equations we're working with, like pick your fighter:= 'N': L * X = B(No transpose)= 'T': L**T * X = B(Transpose)
N(input)INTEGERβ the order of the tridiagonal matrixA. constraint:N >= 0. periodt.NRHS(input)INTEGERβ the number of right hand sides, aka the number of columns of the matrixB. constraint:NRHS >= 0. yes she counted them herself πD(input)REALarray, dimension(N)β thendiagonal elements of the diagonal matrixDfrom the factorization computed bySPTTRF. πE(input)COMPLEXarray, dimension(N-1)β the(n-1)off-diagonal elements of the unit bidiagonal factorUorLfrom the factorization computed bySPTTRF(seeUPLO). π₯B(input/output)COMPLEXarray, dimension(LDB,NRHS)β on entry, the right hand side matrixB. On exit, the solution matrixX. she literally transforms in front of you, iconic πLDB(input)INTEGERβ the leading dimension of the arrayB. constraint:LDB >= max(1,N). do NOT skip this one, she WILL be upset πINFO(output)INTEGERβ the vibe check report:= 0:successful exit β¨< 0:ifINFO = -i, thei-th argument had an illegal value π
bestie listen β ifINFOcomes back negative, something you passed in was not it. go back and look at argument numberi. not me being the one to tell you, but she's judging a little. π
π that's the whole page, slay π
this was a modified LAPACK routine, dated 12 May 1997, and you just read the SPTTRSV(l) page top to bottom. obsessed with you honestly. go solve a triangular system, you've got this πππ§
π¦ 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