claipd()

subroutine claipd	(	integer	N,
		complex, dimension( * )	A,
		integer	INDA,
		integer	VINDA
	)

CLAIPD

Purpose:

 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.

Parameters

[in]	N	N is INTEGER The number of diagonal elements of A.
[in,out]	A	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.
[in]	INDA	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)
[in]	VINDA	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

Author

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.

Date

December 2016

Definition at line 85 of file claipd.f.

*
*  -- LAPACK test routine (version 3.7.0) --
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
*     December 2016
*
*     .. 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

Here is the caller graph for this function: