DSYSV_AA(3)
π NAME
okay so this whole page is about the file /home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/SRC/dsysv_aa.f β that's the source file, bestie. It lives in the Library Functions Manual as dsysv_aa.f(3), and yes, the filename is the main character.
β¨ SYNOPSIS
π§ Functions/Subroutines
subroutine DSYSV_AA (uplo, n, nrhs, a, lda, ipiv, b, ldb, work, lwork,
info)
DSYSV_AA computes the solution to system of linear equations A * X
= B for SY matrices
Like, that's the call shape. DSYSV_AA computes the solution to system of linear equations A * X = B for SY matrices. Periodt.
π Function/Subroutine Documentation
subroutine DSYSV_AA (character uplo, integer n, integer nrhs, double
precision, dimension( lda, * ) a, integer lda, integer, dimension( * )
ipiv, double precision, dimension( ldb, * ) b, integer ldb, double
precision, dimension( * ) work, integer lwork, integer info)
DSYSV_AA computes the solution to system of linear equations A * X = B for SY matrices.
πΈ Purpose
okay so DSYSV β yes, the Purpose text says DSYSV, keep that energy β computes the solution to a real system of linear equations:
A * X = B,
where A is an N-by-N symmetric matrix and X and B are N-by-NRHS matrices.
Aasen's algorithm is used to factor A as
A = U**T * T * U, if UPLO = 'U', or
A = L * T * L**T, if UPLO = 'L',
where U (or L) is a product of permutation and unit upper (lower) triangular matrices, and T is symmetric tridiagonal. The factored form of A is then used to solve the system of equations A * X = B.
tiny source-flavor note: the raw file wraps all this in Doxygen !> comment markers. Cute, but the math is the moment.
π Parameters
UPLOβUPLOisCHARACTER*1.= 'U':upper triangle ofAis stored;= 'L':lower triangle ofAis stored.NβNisINTEGER. The number of linear equations, i.e., the order of the matrixA.N >= 0.NRHSβNRHSisINTEGER. The number of right hand sides, i.e., the number of columns of the matrixB.NRHS >= 0.AβAisDOUBLE PRECISIONarray, dimension(LDA,N). On entry, the symmetric matrixA. IfUPLO = 'U', the leadingN-by-Nupper triangular part ofAcontains the upper triangular part of the matrixA, and the strictly lower triangular part ofAis not referenced. IfUPLO = 'L', the leadingN-by-Nlower triangular part ofAcontains the lower triangular part of the matrixA, and the strictly upper triangular part ofAis not referenced. On exit, ifINFO = 0, the tridiagonal matrixTand the multipliers used to obtain the factorUorLfrom the factorizationA = U**T*T*UorA = L*T*L**Tas computed byDSYTRF.LDAβLDAisINTEGER. The leading dimension of the arrayA.LDA >= max(1,N).IPIVβIPIVisINTEGERarray, dimension(N). On exit, it contains the details of the interchanges, i.e., the row and columnkofAwere interchanged with the row and columnIPIV(k).BβBisDOUBLE PRECISIONarray, dimension(LDB,NRHS). On entry, theN-by-NRHSright hand side matrixB. On exit, ifINFO = 0, theN-by-NRHSsolution matrixX.LDBβLDBisINTEGER. The leading dimension of the arrayB.LDB >= max(1,N).WORKβWORKisDOUBLE PRECISIONarray, dimension(MAX(1,LWORK)). On exit, ifINFO = 0,WORK(1)returns the optimalLWORK.LWORKβLWORKisINTEGER. The length ofWORK.LWORK >= MAX(1,2*N,3*N-2), and for the best performance,LWORK >= MAX(1,N*NB), whereNBis the optimal blocksize forDSYTRF_AA. IfLWORK = -1, then a workspace query is assumed; the routine only calculates the optimal size of theWORKarray, returns this value as the first entry of theWORKarray, and no error message related toLWORKis issued byXERBLA.INFOβINFOisINTEGER.= 0:successful exit.< 0:ifINFO = -i, the i-th argument had an illegal value.> 0:ifINFO = i,D(i,i)is exactly zero. The factorization has been completed, but the block diagonal matrixDis exactly singular, so the solution could not be computed.
π Author
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 160 of file dsysv_aa.f.
π¦ Author
Generated automatically by Doxygen for LAPACK from the source code.
okay bestie, that's the whole page. The footer stamp is LAPACK, Version 3.12.0, and it points back to /home/abui.../dsysv_aa.f(3) β no separate date printed in that footer, that's just the source/version vibe. πβ¨
π¦ where this page comes from
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms