TESTING/LIN/zsyt03.f(3)
okay bestie, welcome to the TESTING/LIN/zsyt03.f(3) page from the Library Functions Manual β yes, the manual name is folded in right here like a little bow π
π NAME π
TESTING/LIN/zsyt03.f
β¨ SYNOPSIS β¨
πΈ Functions/Subroutines πΈ
subroutine zsyt03 (uplo, n, a, lda, ainv, ldainv, work, ldwork, rwork, rcond, resid)βZSYT03
subroutine zsyt03 (uplo, n, a, lda, ainv, ldainv, work, ldwork, rwork,
rcond, resid)
ZSYT03
π Function/Subroutine Documentation π
subroutine zsyt03 (character uplo, integer n, complex*16, dimension( lda, *
) a, integer lda, complex*16, dimension( ldainv, * ) ainv, integer
ldainv, complex*16, dimension( ldwork, * ) work, integer ldwork, double
precision, dimension( * ) rwork, double precision rcond, double
precision resid)
ZSYT03
π₯ Purpose π₯
okay so ZSYT03 computes the residual for a complex symmetric matrix times its inverse:
norm( I - A*AINV ) / ( N * norm(A) * norm(AINV) * EPS )
where EPS is the machine epsilon. that's the whole vibe check π«Ά
π Parameters π
UPLOβUPLO is CHARACTER*1. Specifies whether the upper or lower triangular part of the complex symmetric matrixAis stored:= 'U': Upper triangular= 'L': Lower triangular
NβN is INTEGER. The number of rows and columns of the matrixA.N >= 0.AβA is COMPLEX*16 array, dimension (LDA,N). The original complex symmetric matrixA.LDAβLDA is INTEGER. The leading dimension of the arrayA.LDA >= max(1,N)AINVβAINV is COMPLEX*16 array, dimension (LDAINV,N). On entry, the inverse of the matrixA, stored as a symmetric matrix in the same format asA. In this version,AINVis expanded into a full matrix and multiplied byA, so the opposing triangle ofAINVwill be changed; i.e., if the upper triangular part ofAINVis stored, the lower triangular part will be used as work space.LDAINVβLDAINV is INTEGER. The leading dimension of the arrayAINV.LDAINV >= max(1,N).WORKβWORK is COMPLEX*16 array, dimension (LDWORK,N)LDWORKβLDWORK is INTEGER. The leading dimension of the arrayWORK.LDWORK >= max(1,N).RWORKβRWORK is DOUBLE PRECISION array, dimension (N)RCONDβRCOND is DOUBLE PRECISION. The reciprocal of the condition number ofA, computed asRCOND = 1/ (norm(A) * norm(AINV)).RESIDβRESID is DOUBLE PRECISION.norm(I - A*AINV) / ( N * norm(A) * norm(AINV) * EPS )
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 124 of file zsyt03.f.
π¦ Author π¦
Generated automatically by Doxygen for LAPACK from the source code.
and that's the whole page, bestie! the footer stamp is LAPACK Version 3.12.0 for TESTING/LIN/zsyt03.f(3) β mwah πβ¨
π¦ 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