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 CTIMMM( 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
REAL TIMMIN
* ..
* .. Array Arguments ..
INTEGER LDAVAL( * ), NVAL( * )
REAL RESLTS( LDR1, LDR2, * )
COMPLEX A( * ), B( * ), C( * )
* ..
*
* Purpose
* =======
*
* CTIMMM times CGEMM.
*
* 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) REAL
* The minimum time a subroutine will be timed.
*
* A (workspace) COMPLEX array, dimension (LDAMAX*NMAX)
* where LDAMAX and NMAX are the maximum values permitted
* for LDA and N.
*
* B (workspace) COMPLEX array, dimension (LDAMAX*NMAX)
*
* C (workspace) COMPLEX array, dimension (LDAMAX*NMAX)
*
* RESLTS (output) REAL 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 ONE
PARAMETER ( NSUBS = 1, ONE = ( 1.0E+0, 0.0E+0 ) )
* ..
* .. Local Scalars ..
CHARACTER*6 CNAME
INTEGER I, IC, ICL, ILDA, IN, INFO, ISUB, LDA, N
REAL OPS, S1, S2, TIME, UNTIME
* ..
* .. Local Arrays ..
LOGICAL TIMSUB( NSUBS )
CHARACTER*6 SUBNAM( NSUBS )
INTEGER IDUMMY( 1 )
* ..
* .. External Functions ..
LOGICAL LSAMEN
REAL SECOND, SMFLOP, SOPBL3
EXTERNAL LSAMEN, SECOND, SMFLOP, SOPBL3
* ..
* .. External Subroutines ..
EXTERNAL ATIMCK, CGEMM, CTIMMG, SPRTBL
* ..
* .. Intrinsic Functions ..
INTRINSIC REAL
* ..
* .. Data statements ..
DATA SUBNAM / 'CGEMM ' /
* ..
* .. 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 CGEMM
*
CALL CTIMMG( 1, N, N, A, LDA, 0, 0 )
CALL CTIMMG( 0, N, N, B, LDA, 0, 0 )
CALL CTIMMG( 1, N, N, C, LDA, 0, 0 )
IC = 0
S1 = SECOND( )
30 CONTINUE
CALL CGEMM( 'No transpose', 'No transpose', N, N, N, ONE, A,
$ LDA, B, LDA, ONE, C, LDA )
S2 = SECOND( )
TIME = S2 - S1
IC = IC + 1
IF( TIME.LT.TIMMIN ) THEN
CALL CTIMMG( 1, N, N, C, LDA, 0, 0 )
GO TO 30
END IF
*
* Subtract the time used in CTIMMG.
*
ICL = 1
S1 = SECOND( )
40 CONTINUE
S2 = SECOND( )
UNTIME = S2 - S1
ICL = ICL + 1
IF( ICL.LE.IC ) THEN
CALL CTIMMG( 1, N, N, C, LDA, 0, 0 )
GO TO 40
END IF
*
TIME = ( TIME-UNTIME ) / REAL( IC )
OPS = SOPBL3( 'CGEMM ', N, N, N )
RESLTS( 1, IN, ILDA ) = SMFLOP( 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 SPRTBL( ' ', 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 CTIMMM
*
END
|
