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
|