PLASMA  2.4.5
PLASMA - Parallel Linear Algebra for Scalable Multi-core Architectures
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Groups
zlaipd.f
Go to the documentation of this file.
1  SUBROUTINE zlaipd( N, A, INDA, VINDA )
2 *
3 * -- LAPACK test routine (version 3.1) --
4 * Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
5 * November 2006
6 *
7 * .. Scalar Arguments ..
8  INTEGER inda, n, vinda
9 * ..
10 * .. Array Arguments ..
11  COMPLEX*16 a( * )
12 * ..
13 *
14 * Purpose
15 * =======
16 *
17 * ZLAIPD sets the imaginary part of the diagonal elements of a complex
18 * matrix A to a large value. This is used to test LAPACK routines for
19 * complex Hermitian matrices, which are not supposed to access or use
20 * the imaginary parts of the diagonals.
21 *
22 * Arguments
23 * =========
24 *
25 * N (input) INTEGER
26 * The number of diagonal elements of A.
27 *
28 * A (input/output) COMPLEX*16 array, dimension
29 * (1+(N-1)*INDA+(N-2)*VINDA)
30 * On entry, the complex (Hermitian) matrix A.
31 * On exit, the imaginary parts of the diagonal elements are set
32 * to BIGNUM = EPS / SAFMIN, where EPS is the machine epsilon and
33 * SAFMIN is the safe minimum.
34 *
35 * INDA (input) INTEGER
36 * The increment between A(1) and the next diagonal element of A.
37 * Typical values are
38 * = LDA+1: square matrices with leading dimension LDA
39 * = 2: packed upper triangular matrix, starting at A(1,1)
40 * = N: packed lower triangular matrix, starting at A(1,1)
41 *
42 * VINDA (input) INTEGER
43 * The change in the diagonal increment between columns of A.
44 * Typical values are
45 * = 0: no change, the row and column increments in A are fixed
46 * = 1: packed upper triangular matrix
47 * = -1: packed lower triangular matrix
48 *
49 * =====================================================================
50 *
51 * .. Local Scalars ..
52  INTEGER i, ia, ixa
53  DOUBLE PRECISION bignum
54 * ..
55 * .. External Functions ..
56  DOUBLE PRECISION dlamch
57  EXTERNAL dlamch
58 * ..
59 * .. Intrinsic Functions ..
60  INTRINSIC dble, dcmplx
61 * ..
62 * .. Executable Statements ..
63 *
64  bignum = dlamch( 'Epsilon' ) / dlamch( 'Safe minimum' )
65  ia = 1
66  ixa = inda
67  DO 10 i = 1, n
68  a( ia ) = dcmplx( dble( a( ia ) ), bignum )
69  ia = ia + ixa
70  ixa = ixa + vinda
71  10 continue
72  return
73  END