Julie
e18d437924
GitHub is now enabling us to track accurately version and date. No need for this anymore.
120 lines
3.2 KiB
FortranFixed
120 lines
3.2 KiB
FortranFixed
*> \brief \b CLAIPD
|
|
*
|
|
* =========== DOCUMENTATION ===========
|
|
*
|
|
* Online html documentation available at
|
|
* http://www.netlib.org/lapack/explore-html/
|
|
*
|
|
* Definition:
|
|
* ===========
|
|
*
|
|
* SUBROUTINE CLAIPD( N, A, INDA, VINDA )
|
|
*
|
|
* .. Scalar Arguments ..
|
|
* INTEGER INDA, N, VINDA
|
|
* ..
|
|
* .. Array Arguments ..
|
|
* COMPLEX A( * )
|
|
* ..
|
|
*
|
|
*
|
|
*> \par Purpose:
|
|
* =============
|
|
*>
|
|
*> \verbatim
|
|
*>
|
|
*> CLAIPD sets the imaginary part of the diagonal elements of a complex
|
|
*> matrix A to a large value. This is used to test LAPACK routines for
|
|
*> complex Hermitian matrices, which are not supposed to access or use
|
|
*> the imaginary parts of the diagonals.
|
|
*> \endverbatim
|
|
*
|
|
* Arguments:
|
|
* ==========
|
|
*
|
|
*> \param[in] N
|
|
*> \verbatim
|
|
*> N is INTEGER
|
|
*> The number of diagonal elements of A.
|
|
*> \endverbatim
|
|
*>
|
|
*> \param[in,out] A
|
|
*> \verbatim
|
|
*> A is COMPLEX array, dimension
|
|
*> (1+(N-1)*INDA+(N-2)*VINDA)
|
|
*> On entry, the complex (Hermitian) matrix A.
|
|
*> On exit, the imaginary parts of the diagonal elements are set
|
|
*> to BIGNUM = EPS / SAFMIN, where EPS is the machine epsilon and
|
|
*> SAFMIN is the safe minimum.
|
|
*> \endverbatim
|
|
*>
|
|
*> \param[in] INDA
|
|
*> \verbatim
|
|
*> INDA is INTEGER
|
|
*> The increment between A(1) and the next diagonal element of A.
|
|
*> Typical values are
|
|
*> = LDA+1: square matrices with leading dimension LDA
|
|
*> = 2: packed upper triangular matrix, starting at A(1,1)
|
|
*> = N: packed lower triangular matrix, starting at A(1,1)
|
|
*> \endverbatim
|
|
*>
|
|
*> \param[in] VINDA
|
|
*> \verbatim
|
|
*> VINDA is INTEGER
|
|
*> The change in the diagonal increment between columns of A.
|
|
*> Typical values are
|
|
*> = 0: no change, the row and column increments in A are fixed
|
|
*> = 1: packed upper triangular matrix
|
|
*> = -1: packed lower triangular matrix
|
|
*> \endverbatim
|
|
*
|
|
* Authors:
|
|
* ========
|
|
*
|
|
*> \author Univ. of Tennessee
|
|
*> \author Univ. of California Berkeley
|
|
*> \author Univ. of Colorado Denver
|
|
*> \author NAG Ltd.
|
|
*
|
|
*> \ingroup complex_lin
|
|
*
|
|
* =====================================================================
|
|
SUBROUTINE CLAIPD( N, A, INDA, VINDA )
|
|
*
|
|
* -- LAPACK test routine --
|
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
|
*
|
|
* .. Scalar Arguments ..
|
|
INTEGER INDA, N, VINDA
|
|
* ..
|
|
* .. Array Arguments ..
|
|
COMPLEX A( * )
|
|
* ..
|
|
*
|
|
* =====================================================================
|
|
*
|
|
* .. Local Scalars ..
|
|
INTEGER I, IA, IXA
|
|
REAL BIGNUM
|
|
* ..
|
|
* .. External Functions ..
|
|
REAL SLAMCH
|
|
EXTERNAL SLAMCH
|
|
* ..
|
|
* .. Intrinsic Functions ..
|
|
INTRINSIC CMPLX, REAL
|
|
* ..
|
|
* .. Executable Statements ..
|
|
*
|
|
BIGNUM = SLAMCH( 'Epsilon' ) / SLAMCH( 'Safe minimum' )
|
|
IA = 1
|
|
IXA = INDA
|
|
DO 10 I = 1, N
|
|
A( IA ) = CMPLX( REAL( A( IA ) ), BIGNUM )
|
|
IA = IA + IXA
|
|
IXA = IXA + VINDA
|
|
10 CONTINUE
|
|
RETURN
|
|
END
|