Browse Files
      1
      2
      3
      4
      5
      6
      7
      8
      9
     10
     11
     12
     13
     14
     15
     16
     17
     18
     19
     20
     21
     22
     23
     24
     25
     26
     27
     28
     29
     30
     31
     32
     33
     34
     35
     36
     37
     38
     39
     40
     41
     42
     43
     44
     45
     46
     47
     48
     49
     50
     51
     52
     53
     54
     55
     56
     57
     58
     59
     60
     61
     62
     63
     64
     65
     66
     67
     68
     69
     70
     71
     72
     73
      SUBROUTINE ZLAIPD( N, A, INDA, VINDA )
*
*  -- LAPACK test routine (version 3.1) --
*     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
*     November 2006
*
*     .. Scalar Arguments ..
      INTEGER            INDA, N, VINDA
*     ..
*     .. Array Arguments ..
      COMPLEX*16         A( * )
*     ..
*
*  Purpose
*  =======
*
*  ZLAIPD 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.
*
*  Arguments
*  =========
*
*  N      (input) INTEGER
*         The number of diagonal elements of A.
*
*  A      (input/output) COMPLEX*16 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.
*
*  INDA   (input) 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)
*
*  VINDA  (input) 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
*
*  =====================================================================
*
*     .. Local Scalars ..
      INTEGER            I, IA, IXA
      DOUBLE PRECISION   BIGNUM
*     ..
*     .. External Functions ..
      DOUBLE PRECISION   DLAMCH
      EXTERNAL           DLAMCH
*     ..
*     .. Intrinsic Functions ..
      INTRINSIC          DBLE, DCMPLX
*     ..
*     .. Executable Statements ..
*
      BIGNUM = DLAMCH( 'Epsilon' ) / DLAMCH( 'Safe minimum' )
      IA = 1
      IXA = INDA
      DO 10 I = 1, N
         A( IA ) = DCMPLX( DBLE( A( IA ) ), BIGNUM )
         IA = IA + IXA
         IXA = IXA + VINDA
   10 CONTINUE
      RETURN
      END