--- rpl/lapack/lapack/zla_syrcond_x.f 2010/08/07 13:21:10 1.1 +++ rpl/lapack/lapack/zla_syrcond_x.f 2023/08/07 08:39:28 1.18 @@ -1,17 +1,140 @@ +*> \brief \b ZLA_SYRCOND_X computes the infinity norm condition number of op(A)*diag(x) for symmetric indefinite matrices. +* +* =========== DOCUMENTATION =========== +* +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ +* +*> \htmlonly +*> Download ZLA_SYRCOND_X + dependencies +*> +*> [TGZ] +*> +*> [ZIP] +*> +*> [TXT] +*> \endhtmlonly +* +* Definition: +* =========== +* +* DOUBLE PRECISION FUNCTION ZLA_SYRCOND_X( UPLO, N, A, LDA, AF, +* LDAF, IPIV, X, INFO, +* WORK, RWORK ) +* +* .. Scalar Arguments .. +* CHARACTER UPLO +* INTEGER N, LDA, LDAF, INFO +* .. +* .. Array Arguments .. +* INTEGER IPIV( * ) +* COMPLEX*16 A( LDA, * ), AF( LDAF, * ), WORK( * ), X( * ) +* DOUBLE PRECISION RWORK( * ) +* .. +* +* +*> \par Purpose: +* ============= +*> +*> \verbatim +*> +*> ZLA_SYRCOND_X Computes the infinity norm condition number of +*> op(A) * diag(X) where X is a COMPLEX*16 vector. +*> \endverbatim +* +* Arguments: +* ========== +* +*> \param[in] UPLO +*> \verbatim +*> UPLO is CHARACTER*1 +*> = 'U': Upper triangle of A is stored; +*> = 'L': Lower triangle of A is stored. +*> \endverbatim +*> +*> \param[in] N +*> \verbatim +*> N is INTEGER +*> The number of linear equations, i.e., the order of the +*> matrix A. N >= 0. +*> \endverbatim +*> +*> \param[in] A +*> \verbatim +*> A is COMPLEX*16 array, dimension (LDA,N) +*> On entry, the N-by-N matrix A. +*> \endverbatim +*> +*> \param[in] LDA +*> \verbatim +*> LDA is INTEGER +*> The leading dimension of the array A. LDA >= max(1,N). +*> \endverbatim +*> +*> \param[in] AF +*> \verbatim +*> AF is COMPLEX*16 array, dimension (LDAF,N) +*> The block diagonal matrix D and the multipliers used to +*> obtain the factor U or L as computed by ZSYTRF. +*> \endverbatim +*> +*> \param[in] LDAF +*> \verbatim +*> LDAF is INTEGER +*> The leading dimension of the array AF. LDAF >= max(1,N). +*> \endverbatim +*> +*> \param[in] IPIV +*> \verbatim +*> IPIV is INTEGER array, dimension (N) +*> Details of the interchanges and the block structure of D +*> as determined by ZSYTRF. +*> \endverbatim +*> +*> \param[in] X +*> \verbatim +*> X is COMPLEX*16 array, dimension (N) +*> The vector X in the formula op(A) * diag(X). +*> \endverbatim +*> +*> \param[out] INFO +*> \verbatim +*> INFO is INTEGER +*> = 0: Successful exit. +*> i > 0: The ith argument is invalid. +*> \endverbatim +*> +*> \param[out] WORK +*> \verbatim +*> WORK is COMPLEX*16 array, dimension (2*N). +*> Workspace. +*> \endverbatim +*> +*> \param[out] RWORK +*> \verbatim +*> RWORK is DOUBLE PRECISION array, dimension (N). +*> Workspace. +*> \endverbatim +* +* Authors: +* ======== +* +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. +* +*> \ingroup complex16SYcomputational +* +* ===================================================================== DOUBLE PRECISION FUNCTION ZLA_SYRCOND_X( UPLO, N, A, LDA, AF, $ LDAF, IPIV, X, INFO, $ WORK, RWORK ) * -* -- LAPACK routine (version 3.2.1) -- -* -- Contributed by James Demmel, Deaglan Halligan, Yozo Hida and -- -* -- Jason Riedy of Univ. of California Berkeley. -- -* -- April 2009 -- +* -- LAPACK computational routine -- +* -- LAPACK is a software package provided by Univ. of Tennessee, -- +* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- * -* -- LAPACK is a software package provided by Univ. of Tennessee, -- -* -- Univ. of California Berkeley and NAG Ltd. -- -* - IMPLICIT NONE -* .. * .. Scalar Arguments .. CHARACTER UPLO INTEGER N, LDA, LDAF, INFO @@ -22,60 +145,13 @@ DOUBLE PRECISION RWORK( * ) * .. * -* Purpose -* ======= -* -* ZLA_SYRCOND_X Computes the infinity norm condition number of -* op(A) * diag(X) where X is a COMPLEX*16 vector. -* -* Arguments -* ========= -* -* UPLO (input) CHARACTER*1 -* = 'U': Upper triangle of A is stored; -* = 'L': Lower triangle of A is stored. -* -* N (input) INTEGER -* The number of linear equations, i.e., the order of the -* matrix A. N >= 0. -* -* A (input) COMPLEX*16 array, dimension (LDA,N) -* On entry, the N-by-N matrix A. -* -* LDA (input) INTEGER -* The leading dimension of the array A. LDA >= max(1,N). -* -* AF (input) COMPLEX*16 array, dimension (LDAF,N) -* The block diagonal matrix D and the multipliers used to -* obtain the factor U or L as computed by ZSYTRF. -* -* LDAF (input) INTEGER -* The leading dimension of the array AF. LDAF >= max(1,N). -* -* IPIV (input) INTEGER array, dimension (N) -* Details of the interchanges and the block structure of D -* as determined by ZSYTRF. -* -* X (input) COMPLEX*16 array, dimension (N) -* The vector X in the formula op(A) * diag(X). -* -* INFO (output) INTEGER -* = 0: Successful exit. -* i > 0: The ith argument is invalid. -* -* WORK (input) COMPLEX*16 array, dimension (2*N). -* Workspace. -* -* RWORK (input) DOUBLE PRECISION array, dimension (N). -* Workspace. -* * ===================================================================== * * .. Local Scalars .. INTEGER KASE DOUBLE PRECISION AINVNM, ANORM, TMP INTEGER I, J - LOGICAL UP + LOGICAL UP, UPPER COMPLEX*16 ZDUM * .. * .. Local Arrays .. @@ -102,8 +178,15 @@ ZLA_SYRCOND_X = 0.0D+0 * INFO = 0 - IF( N.LT.0 ) THEN + UPPER = LSAME( UPLO, 'U' ) + IF( .NOT.UPPER .AND. .NOT.LSAME( UPLO, 'L' ) ) THEN + INFO = -1 + ELSE IF( N.LT.0 ) THEN INFO = -2 + ELSE IF( LDA.LT.MAX( 1, N ) ) THEN + INFO = -4 + ELSE IF( LDAF.LT.MAX( 1, N ) ) THEN + INFO = -6 END IF IF( INFO.NE.0 ) THEN CALL XERBLA( 'ZLA_SYRCOND_X', -INFO ) @@ -181,7 +264,7 @@ END DO ELSE * -* Multiply by inv(X'). +* Multiply by inv(X**T). * DO I = 1, N WORK( I ) = WORK( I ) / X( I ) @@ -211,4 +294,6 @@ * RETURN * +* End of ZLA_SYRCOND_X +* END