DGET53(3)
DGET53(3) π
okay so this is DGET53(3), straight out of the Library Functions Manual β and yes, bestie, we're about to get into some eigenvalue-checking lore. π¦
π NAME π
the name on this page is literally this source path, no cap:
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/TESTING/EIG/dget53.f
β¨ SYNOPSIS β¨
Functions/Subroutines
subroutine DGET53 (a, lda, b, ldb, scale, wr, wi, result, info)
DGET53
π Function/Subroutine Documentation π
here's the full signature, bestie:
subroutine DGET53 (double precision, dimension( lda, * ) a, integer lda,
double precision, dimension( ldb, * ) b, integer ldb, double precision
scale, double precision wr, double precision wi, double precision
result, integer info)
DGET53
π― Purpose
okay so DGET53 checks the generalized eigenvalues computed by DLAG2. iconic behavior. π
The basic test for an eigenvalue is this absolute cinema:
| det( s A - w B ) |
RESULT = ---------------------------------------------------
ulp max( s norm(A), |w| norm(B) )*norm( s A - w B )
Two are performed:
ulp*max( s*norm(A), |w|*norm(B) )must be at leastsafe_minimum. This insures that the test performed is not essentiallydet(0*A + 0*B)=0.s*norm(A) + |w|*norm(B)must be less than1/safe_minimum. This insures thats*A - w*Bwill not overflow.
If these tests are not passed, then s and w are scaled and tested anyway, if this is possible. resilient queen energy. β¨
π§Έ Parameters
AβAisDOUBLE PRECISIONarray, dimension(LDA, 2). The 2x2 matrixA.LDAβLDAisINTEGER. The leading dimension ofA. It must be at least2.BβBisDOUBLE PRECISIONarray, dimension(LDB, N). The 2x2 upper-triangular matrixB.LDBβLDBisINTEGER. The leading dimension ofB. It must be at least2.SCALEβSCALEisDOUBLE PRECISION. Thesin the formulas A - w B. It is assumed to be non-negative.WRβWRisDOUBLE PRECISION. The real part of the eigenvaluewin the formulas A - w B.WIβWIisDOUBLE PRECISION. The imaginary part of the eigenvaluewin the formulas A - w B.RESULTβRESULTisDOUBLE PRECISION. IfINFOis2or less, the value computed by the test described above. IfINFO=3, this will just be1/ulp.INFOβINFOisINTEGER. Here's the tea:=0: The input data pass the . (yes, the page saysthe .β mysterious main character behavior)=1:s*norm(A) + |w|*norm(B) > 1/safe_minimum.=2:ulp*max( s*norm(A), |w|*norm(B) ) < safe_minimum=3: same asINFO=2, butsandwcould not be scaled so as to compute the test.
π Author
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 125 of file dget53.f.
πΈ AUTHOR πΈ
This page was generated automatically by Doxygen for LAPACK from the source code.
Catch you later, bestie β this was LAPACK Version 3.12.0, from /home/abui...IG/dget53.f(3). The footer didn't give us a date, so that's the whole story. Mwah πβ¨
π¦ 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