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
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
|
SUBROUTINE ZTIMMM( VNAME, LAB2, NN, NVAL, NLDA, LDAVAL, TIMMIN, A,
$ B, C, RESLTS, LDR1, LDR2, NOUT )
*
* -- LAPACK timing routine (version 3.0) --
* Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd.,
* Courant Institute, Argonne National Lab, and Rice University
* September 30, 1994
*
* .. Scalar Arguments ..
CHARACTER*( * ) LAB2, VNAME
INTEGER LDR1, LDR2, NLDA, NN, NOUT
DOUBLE PRECISION TIMMIN
* ..
* .. Array Arguments ..
INTEGER LDAVAL( * ), NVAL( * )
DOUBLE PRECISION RESLTS( LDR1, LDR2, * )
COMPLEX*16 A( * ), B( * ), C( * )
* ..
*
* Purpose
* =======
*
* ZTIMMM times ZGEMM.
*
* Arguments
* =========
*
* VNAME (input) CHARACTER*(*)
* The name of the Level 3 BLAS routine to be timed.
*
* LAB2 (input) CHARACTER*(*)
* The name of the variable given in NVAL.
*
* NN (input) INTEGER
* The number of values of N contained in the vector NVAL.
*
* NVAL (input) INTEGER array, dimension (NN)
* The values of the matrix dimension N.
*
* NLDA (input) INTEGER
* The number of values of LDA contained in the vector LDAVAL.
*
* LDAVAL (input) INTEGER array, dimension (NLDA)
* The values of the leading dimension of the array A.
*
* TIMMIN (input) DOUBLE PRECISION
* The minimum time a subroutine will be timed.
*
* A (workspace) COMPLEX*16 array, dimension (LDAMAX*NMAX)
* where LDAMAX and NMAX are the maximum values permitted
* for LDA and N.
*
* B (workspace) COMPLEX*16 array, dimension (LDAMAX*NMAX)
*
* C (workspace) COMPLEX*16 array, dimension (LDAMAX*NMAX)
*
* RESLTS (output) DOUBLE PRECISION array, dimension (LDR1,LDR2,NLDA)
* The timing results for each subroutine over the relevant
* values of N and LDA.
*
* LDR1 (input) INTEGER
* The first dimension of RESLTS. LDR1 >= 1.
*
* LDR2 (input) INTEGER
* The second dimension of RESLTS. LDR2 >= max(1,NN).
*
* NOUT (input) INTEGER
* The unit number for output.
*
* =====================================================================
*
* .. Parameters ..
INTEGER NSUBS
COMPLEX*16 ONE
PARAMETER ( NSUBS = 1, ONE = ( 1.0D+0, 0.0D+0 ) )
* ..
* .. Local Scalars ..
CHARACTER*6 CNAME
INTEGER I, IC, ICL, ILDA, IN, INFO, ISUB, LDA, N
DOUBLE PRECISION OPS, S1, S2, TIME, UNTIME
* ..
* .. Local Arrays ..
LOGICAL TIMSUB( NSUBS )
CHARACTER*6 SUBNAM( NSUBS )
INTEGER IDUMMY( 1 )
* ..
* .. External Functions ..
LOGICAL LSAMEN
DOUBLE PRECISION DMFLOP, DOPBL3, DSECND
EXTERNAL LSAMEN, DMFLOP, DOPBL3, DSECND
* ..
* .. External Subroutines ..
EXTERNAL ATIMCK, DPRTBL, ZGEMM, ZTIMMG
* ..
* .. Intrinsic Functions ..
INTRINSIC DBLE
* ..
* .. Data statements ..
DATA SUBNAM / 'ZGEMM ' /
* ..
* .. Executable Statements ..
*
CNAME = VNAME
DO 10 ISUB = 1, NSUBS
TIMSUB( ISUB ) = LSAMEN( 6, CNAME, SUBNAM( ISUB ) )
IF( TIMSUB( ISUB ) )
$ GO TO 20
10 CONTINUE
WRITE( NOUT, FMT = 9999 )CNAME
GO TO 80
20 CONTINUE
*
* Check that N <= LDA for the input values.
*
CALL ATIMCK( 2, CNAME, NN, NVAL, NLDA, LDAVAL, NOUT, INFO )
IF( INFO.GT.0 ) THEN
WRITE( NOUT, FMT = 9998 )CNAME
GO TO 80
END IF
*
DO 60 ILDA = 1, NLDA
LDA = LDAVAL( ILDA )
DO 50 IN = 1, NN
N = NVAL( IN )
*
* Time ZGEMM
*
CALL ZTIMMG( 1, N, N, A, LDA, 0, 0 )
CALL ZTIMMG( 0, N, N, B, LDA, 0, 0 )
CALL ZTIMMG( 1, N, N, C, LDA, 0, 0 )
IC = 0
S1 = DSECND( )
30 CONTINUE
CALL ZGEMM( 'No transpose', 'No transpose', N, N, N, ONE, A,
$ LDA, B, LDA, ONE, C, LDA )
S2 = DSECND( )
TIME = S2 - S1
IC = IC + 1
IF( TIME.LT.TIMMIN ) THEN
CALL ZTIMMG( 1, N, N, C, LDA, 0, 0 )
GO TO 30
END IF
*
* Subtract the time used in ZTIMMG.
*
ICL = 1
S1 = DSECND( )
40 CONTINUE
S2 = DSECND( )
UNTIME = S2 - S1
ICL = ICL + 1
IF( ICL.LE.IC ) THEN
CALL ZTIMMG( 1, N, N, C, LDA, 0, 0 )
GO TO 40
END IF
*
TIME = ( TIME-UNTIME ) / DBLE( IC )
OPS = DOPBL3( 'ZGEMM ', N, N, N )
RESLTS( 1, IN, ILDA ) = DMFLOP( OPS, TIME, 0 )
50 CONTINUE
60 CONTINUE
*
* Print the table of results on unit NOUT.
*
WRITE( NOUT, FMT = 9997 )VNAME
IF( NLDA.EQ.1 ) THEN
WRITE( NOUT, FMT = 9996 )LDAVAL( 1 )
ELSE
DO 70 I = 1, NLDA
WRITE( NOUT, FMT = 9995 )I, LDAVAL( I )
70 CONTINUE
END IF
WRITE( NOUT, FMT = * )
CALL DPRTBL( ' ', LAB2, 1, IDUMMY, NN, NVAL, NLDA, RESLTS, LDR1,
$ LDR2, NOUT )
*
80 CONTINUE
RETURN
9999 FORMAT( 1X, A6, ': Unrecognized path or subroutine name', / )
9998 FORMAT( 1X, A6, ' timing run not attempted', / )
9997 FORMAT( / ' *** Speed of ', A6, ' in megaflops ***' )
9996 FORMAT( 5X, 'with LDA = ', I5 )
9995 FORMAT( 5X, 'line ', I2, ' with LDA = ', I5 )
*
* End of ZTIMMM
*
END
|
