TESTING/EIG/dget54.f(3) | Library Functions Manual | TESTING/EIG/dget54.f(3) |
NAME
TESTING/EIG/dget54.f
SYNOPSIS
Functions/Subroutines
subroutine dget54 (n, a, lda, b, ldb, s, lds, t, ldt, u,
ldu, v, ldv, work, result)
DGET54
Function/Subroutine Documentation
subroutine dget54 (integer n, double precision, dimension( lda, * ) a, integer lda, double precision, dimension( ldb, * ) b, integer ldb, double precision, dimension( lds, * ) s, integer lds, double precision, dimension( ldt, * ) t, integer ldt, double precision, dimension( ldu, * ) u, integer ldu, double precision, dimension( ldv, * ) v, integer ldv, double precision, dimension( * ) work, double precision result)
DGET54
Purpose:
DGET54 checks a generalized decomposition of the form A = U*S*V' and B = U*T* V' where ' means transpose and U and V are orthogonal. Specifically, RESULT = ||( A - U*S*V', B - U*T*V' )|| / (||( A, B )||*n*ulp )
Parameters
N
N is INTEGER The size of the matrix. If it is zero, DGET54 does nothing. It must be at least zero.
A
A is DOUBLE PRECISION array, dimension (LDA, N) The original (unfactored) matrix A.
LDA
LDA is INTEGER The leading dimension of A. It must be at least 1 and at least N.
B
B is DOUBLE PRECISION array, dimension (LDB, N) The original (unfactored) matrix B.
LDB
LDB is INTEGER The leading dimension of B. It must be at least 1 and at least N.
S
S is DOUBLE PRECISION array, dimension (LDS, N) The factored matrix S.
LDS
LDS is INTEGER The leading dimension of S. It must be at least 1 and at least N.
T
T is DOUBLE PRECISION array, dimension (LDT, N) The factored matrix T.
LDT
LDT is INTEGER The leading dimension of T. It must be at least 1 and at least N.
U
U is DOUBLE PRECISION array, dimension (LDU, N) The orthogonal matrix on the left-hand side in the decomposition.
LDU
LDU is INTEGER The leading dimension of U. LDU must be at least N and at least 1.
V
V is DOUBLE PRECISION array, dimension (LDV, N) The orthogonal matrix on the left-hand side in the decomposition.
LDV
LDV is INTEGER The leading dimension of V. LDV must be at least N and at least 1.
WORK
WORK is DOUBLE PRECISION array, dimension (3*N**2)
RESULT
RESULT is DOUBLE PRECISION The value RESULT, It is currently limited to 1/ulp, to avoid overflow. Errors are flagged by RESULT=10/ulp.
Author
Univ. of Tennessee
Univ. of California Berkeley
Univ. of Colorado Denver
NAG Ltd.
Definition at line 154 of file dget54.f.
Author
Generated automatically by Doxygen for LAPACK from the source code.
Version 3.12.0 | LAPACK |