TESTING/EIG/zget22.f(3)
π TESTING/EIG/zget22.f(3) β Library Functions Manual π
okay so bestie, gather round, because this is the manual page for TESTING/EIG/zget22.f and honestly? it's giving "eigenvector quality control" and I'm lowkey obsessed π¦
β¨ NAME β¨
TESTING/EIG/zget22.f β yes that's the whole name, it's just out here living its best file-path life, no cap π
πΈ SYNOPSIS πΈ
Functions/Subroutines
subroutine zget22 (transa, transe, transw, n, a, lda, e, lde, w, work,
rwork, result)
ZGET22 β that's her, that's the icon, that's the whole subroutine π
π Function/Subroutine Documentation π
subroutine zget22 (character transa, character transe, character transw,
integer n, complex*16, dimension( lda, * ) a, integer lda, complex*16,
dimension( lde, * ) e, integer lde, complex*16, dimension( * ) w,
complex*16, dimension( * ) work, double precision, dimension( * ) rwork,
double precision, dimension( 2 ) result)
ZGET22 π
π§ Purpose:
!>!> ZGET22 does an eigenvector check.!>!> The basic test is:!>!> RESULT(1) = | A E - E W | / ( |A| |E| ulp )!>!> using the 1-norm. It also tests the normalization of E:!>!> RESULT(2) = max | m-norm(E(j)) - 1 | / ( n ulp )!> j!>!> where E(j) is the j-th eigenvector, and m-norm is the max-norm of a!> vector. The max-norm of a complex n-vector x in this case is the!> maximum of |re(x(i)| + |im(x(i)| over i = 1, ..., n.!>
so basically, bestie, it does an eigenvector check β two whole tests, no cap. Test one (RESULT(1)) checks | A E - E W | / ( |A| |E| ulp ) using the 1-norm. Test two (RESULT(2)) checks the normalization of E with max | m-norm(E(j)) - 1 | / ( n ulp ) over j. and E(j) is the j-th eigenvector, m-norm is the max-norm of a vector, and for a complex n-vector x that max-norm is literally the maximum of |re(x(i)| + |im(x(i)| over i = 1, ..., n. periodt β¨
π Parameters π
TRANSAβ it's aCHARACTER*1and it says whether A is transposed or not, bestie:= 'N': No transpose= 'T': Transpose= 'C': Conjugate transpose
TRANSEβCHARACTER*1, controls whether E is transposed:= 'N': No transpose, eigenvectors are in columns of E= 'T': Transpose, eigenvectors are in rows of E= 'C': Conjugate transpose, eigenvectors are in rows of E
TRANSWβCHARACTER*1, whether W gets transposed:= 'N': No transpose= 'T': Transpose, same as TRANSW = 'N'(yes, literally the same, don't ask π)= 'C': Conjugate transpose, use -WI(j) instead of WI(j)
NβINTEGER, the order of the matrix A.N >= 0. keep it non-negative or it's a no from me πAβCOMPLEX*16 array, dimension (LDA,N). it's the matrix whose eigenvectors are chilling in E π«ΆLDAβINTEGER, the leading dimension of the array A.LDA >= max(1,N). respect the constraint fr frEβCOMPLEX*16 array, dimension (LDE,N), the matrix of eigenvectors. ifTRANSE = 'N'the eigenvectors are stored in the columns of E; ifTRANSE = 'T'or'C', they're stored in the rows of E. iconic either way πLDEβINTEGER, the leading dimension of the array E.LDE >= max(1,N).WβCOMPLEX*16 array, dimension (N), the eigenvalues of A πWORKβCOMPLEX*16 array, dimension (N*N)RWORKβDOUBLE PRECISION array, dimension (N)RESULTβDOUBLE PRECISION array, dimension (2), and it holds:RESULT(1) = | A E - E W | / ( |A| |E| ulp )RESULT(2) = max | m-norm(E(j)) - 1 | / ( n ulp )j
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
the academic royalty of it all, honestly π
β Author (again, but different) β
Generated automatically by Doxygen for LAPACK from the source code. so she didn't even write herself, she was born from the source code, main character energy for real π¦
π Definition
Definition at line 142 of file zget22.f. β that's where she lives, line 142, remember it π
ok that's the whole page bestie, you ate that π
LAPACK β Version 3.12.0 π bye bestie, go check those eigenvectors, and remember: the max-norm is the maximum of |re(x(i)| + |im(x(i)| over i = 1, ..., n, no cap πβ¨
π¦ 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