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
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
|
*> \brief \b ZLAIPD
*
* =========== DOCUMENTATION ===========
*
* Online html documentation available at
* http://www.netlib.org/lapack/explore-html/
*
* Definition:
* ===========
*
* SUBROUTINE ZLAIPD( N, A, INDA, VINDA )
*
* .. Scalar Arguments ..
* INTEGER INDA, N, VINDA
* ..
* .. Array Arguments ..
* COMPLEX*16 A( * )
* ..
*
*
*> \par Purpose:
* =============
*>
*> \verbatim
*>
*> 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.
*> \endverbatim
*
* Arguments:
* ==========
*
*> \param[in] N
*> \verbatim
*> N is INTEGER
*> The number of diagonal elements of A.
*> \endverbatim
*>
*> \param[in,out] A
*> \verbatim
*> A is 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.
*> \endverbatim
*>
*> \param[in] INDA
*> \verbatim
*> 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)
*> \endverbatim
*>
*> \param[in] VINDA
*> \verbatim
*> 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
*> \endverbatim
*
* Authors:
* ========
*
*> \author Univ. of Tennessee
*> \author Univ. of California Berkeley
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
*> \date December 2016
*
*> \ingroup complex16_lin
*
* =====================================================================
SUBROUTINE ZLAIPD( N, A, INDA, VINDA )
*
* -- 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*16 A( * )
* ..
*
* =====================================================================
*
* .. 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
|
