BLAS/SRC/dspr.f(3)
π BLAS/SRC/dspr.f(3) β Library Functions Manual π
okay bestie, buckle up β this is the man page for BLAS/SRC/dspr.f, straight out of the Library Functions Manual, and honestly? she's iconic. let's get into it β¨
πΈ NAME πΈ
BLAS/SRC/dspr.f β that's the file, that's the vibe, that's the whole name π
β¨ SYNOPSIS β¨
π¦ Functions/Subroutines
subroutine dspr (uplo, n, alpha, x, incx, ap)
DSPR β literally the main character of this page π
π§ Function/Subroutine Documentation π§
subroutine dspr (character uplo, integer n, double precision alpha, double
precision, dimension(*) x, integer incx, double precision, dimension(*)
ap)
DSPR
π Purpose:
!>
!> DSPR performs the symmetric rank 1 operation
!>
!> A := alpha*x*x**T + A,
!>
!> where alpha is a real scalar, x is an n element vector and A is an
!> n by n symmetric matrix, supplied in packed form.
!>
so basically she does a symmetric rank 1 operation, no cap π
the formula is right there: A := alpha*x*x**T + A, where alpha is a real scalar, x is an n element vector, and A is an n by n symmetric matrix supplied in packed form. slay.
π Parameters π
UPLOβUPLOisCHARACTER*1. on entry,UPLOtells you whether the upper or lower triangular part of the matrixAis supplied in the packed arrayAP, like this:UPLO = 'U'or'u'β the upper triangular part ofAis supplied inAP.UPLO = 'L'or'l'β the lower triangular part ofAis supplied inAP.
NβNisINTEGER. on entry,Nspecifies the order of the matrixA.Nmust be at least zero. (zero is allowed, bestie, she's not gonna judge π«Ά)ALPHAβALPHAisDOUBLE PRECISION.on entry,ALPHAspecifies the scalar alpha. πXβXisDOUBLE PRECISION array, dimension at least( 1 + ( n - 1 )*abs( INCX ) ). before entry, the incremented arrayXmust contain the n element vectorx. πINCXβINCXisINTEGER. on entry,INCXspecifies the increment for the elements ofX.INCXmust not be zero. β οΈ like, genuinely not zero. don't do it.APβAPisDOUBLE PRECISION array, dimension at least( ( n*( n + 1 ) )/2 ). before entry withUPLO = 'U'or'u', the arrayAPmust contain the upper triangular part of the symmetric matrix packed sequentially, column by column, so thatAP( 1 )containsa( 1, 1 ),AP( 2 )andAP( 3 )containa( 1, 2 )anda( 2, 2 )respectively, and so on. on exit, the arrayAPis overwritten by the upper triangular part of the updated matrix. π before entry withUPLO = 'L'or'l', the arrayAPmust contain the lower triangular part of the symmetric matrix packed sequentially, column by column, so thatAP( 1 )containsa( 1, 1 ),AP( 2 )andAP( 3 )containa( 2, 1 )anda( 3, 1 )respectively, and so on. on exit, the arrayAPis overwritten by the lower triangular part of the updated matrix. π₯
π Author π
the legends themselves, fr fr:
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
β Further Details:
!>
!> Level 2 Blas routine.
!>
!> -- Written on 22-October-1986.
!> Jack Dongarra, Argonne National Lab.
!> Jeremy Du Croz, Nag Central Office.
!> Sven Hammarling, Nag Central Office.
!> Richard Hanson, Sandia National Labs.
!>
Level 2 Blas routine, written on 22-October-1986 by Jack Dongarra (Argonne National Lab), Jeremy Du Croz (Nag Central Office), Sven Hammarling (Nag Central Office) and Richard Hanson (Sandia National Labs). vintage β¨ we stan an old-school icon.
π Definition π
definition at line 126 of file dspr.f π
π«Ά Author π«Ά
generated automatically by Doxygen for LAPACK from the source code. not a human, but she still ate π
ngl bestie, if you're passing INCX = 0 that's on you β the page literally says it must not be zero π
LAPACK β Version 3.12.0 β BLAS/SRC/dspr.f(3)
and that's the whole page, bestie β mwah, go forth and pack those triangular arrays πβ¨ππ¦πΈπ π«Άππ₯ππβππππ§π
π¦ where this page comes from
- arch rolling: lapack-doc 3.12.1-2 Β· custom
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms