#include "blaswrap.h"
/*  -- translated by f2c (version 19990503).
   You must link the resulting object file with the libraries:
	-lf2c -lm   (in that order)
*/

#include "f2c.h"

/* Subroutine */ int claipd_(integer *n, complex *a, integer *inda, integer *
	vinda)
{
    /* System generated locals */
    integer i__1, i__2, i__3;
    real r__1;
    complex q__1;

    /* Local variables */
    static integer i__, ia;
    extern doublereal slamch_(char *);
    static real bignum;
    static integer ixa;


/*  -- LAPACK test routine (version 3.0) --   
       Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd.,   
       Courant Institute, Argonne National Lab, and Rice University   
       September 30, 1994   


    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.   

    Arguments   
    =========   

    N      (input) INTEGER   
           The number of diagonal elements of A.   

    A      (input/output) 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.   

    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   

    =====================================================================   


       Parameter adjustments */
    --a;

    /* Function Body */
    bignum = slamch_("Epsilon") / slamch_("Safe minimum");
    ia = 1;
    ixa = *inda;
    i__1 = *n;
    for (i__ = 1; i__ <= i__1; ++i__) {
	i__2 = ia;
	i__3 = ia;
	r__1 = a[i__3].r;
	q__1.r = r__1, q__1.i = bignum;
	a[i__2].r = q__1.r, a[i__2].i = q__1.i;
	ia += ixa;
	ixa += *vinda;
/* L10: */
    }
    return 0;
} /* claipd_ */

