.TH "TESTING/EIG/dget53.f" 3 "Version 3.12.0" "LAPACK" \" -*- nroff -*- .ad l .nh .SH NAME TESTING/EIG/dget53.f .SH SYNOPSIS .br .PP .SS "Functions/Subroutines" .in +1c .ti -1c .RI "subroutine \fBdget53\fP (a, lda, b, ldb, scale, wr, wi, result, info)" .br .RI "\fBDGET53\fP " .in -1c .SH "Function/Subroutine Documentation" .PP .SS "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)" .PP \fBDGET53\fP .PP \fBPurpose:\fP .RS 4 .PP .nf !> !> DGET53 checks the generalized eigenvalues computed by DLAG2\&. !> !> The basic test for an eigenvalue is: !> !> | det( s A - w B ) | !> RESULT = --------------------------------------------------- !> ulp max( s norm(A), |w| norm(B) )*norm( s A - w B ) !> !> Two are performed: !> !> (1) ulp*max( s*norm(A), |w|*norm(B) ) must be at least !> safe_minimum\&. This insures that the test performed is !> not essentially det(0*A + 0*B)=0\&. !> !> (2) s*norm(A) + |w|*norm(B) must be less than 1/safe_minimum\&. !> This insures that s*A - w*B will not overflow\&. !> !> If these tests are not passed, then s and w are scaled and !> tested anyway, if this is possible\&. !> .fi .PP .RE .PP \fBParameters\fP .RS 4 \fIA\fP .PP .nf !> A is DOUBLE PRECISION array, dimension (LDA, 2) !> The 2x2 matrix A\&. !> .fi .PP .br \fILDA\fP .PP .nf !> LDA is INTEGER !> The leading dimension of A\&. It must be at least 2\&. !> .fi .PP .br \fIB\fP .PP .nf !> B is DOUBLE PRECISION array, dimension (LDB, N) !> The 2x2 upper-triangular matrix B\&. !> .fi .PP .br \fILDB\fP .PP .nf !> LDB is INTEGER !> The leading dimension of B\&. It must be at least 2\&. !> .fi .PP .br \fISCALE\fP .PP .nf !> SCALE is DOUBLE PRECISION !> The s in the formula s A - w B \&. It is !> assumed to be non-negative\&. !> .fi .PP .br \fIWR\fP .PP .nf !> WR is DOUBLE PRECISION !> The real part of the eigenvalue w in the formula !> s A - w B \&. !> .fi .PP .br \fIWI\fP .PP .nf !> WI is DOUBLE PRECISION !> The imaginary part of the eigenvalue w in the formula !> s A - w B \&. !> .fi .PP .br \fIRESULT\fP .PP .nf !> RESULT is DOUBLE PRECISION !> If INFO is 2 or less, the value computed by the test !> described above\&. !> If INFO=3, this will just be 1/ulp\&. !> .fi .PP .br \fIINFO\fP .PP .nf !> INFO is INTEGER !> =0: The input data pass the \&. !> =1: s*norm(A) + |w|*norm(B) > 1/safe_minimum\&. !> =2: ulp*max( s*norm(A), |w|*norm(B) ) < safe_minimum !> =3: same as INFO=2, but s and w could not be scaled so !> as to compute the test\&. !> .fi .PP .RE .PP \fBAuthor\fP .RS 4 Univ\&. of Tennessee .PP Univ\&. of California Berkeley .PP Univ\&. of Colorado Denver .PP NAG Ltd\&. .RE .PP .PP Definition at line \fB125\fP of file \fBdget53\&.f\fP\&. .SH "Author" .PP Generated automatically by Doxygen for LAPACK from the source code\&.