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
|
/* -*- mode: C; tab-width: 2; indent-tabs-mode: nil; fill-column: 79; coding: iso-latin-1-unix -*- */
#include <hpcc.h>
int
HPCC_StarDGEMM(HPCC_Params *params) {
int myRank, commSize;
double localGflops, minGflops, maxGflops, avgGflops;
int n;
int rv, errCount, failure, failureAll;
FILE *outputFile;
MPI_Comm comm = MPI_COMM_WORLD;
localGflops = minGflops = maxGflops = avgGflops = 0.0;
MPI_Comm_size( comm, &commSize );
MPI_Comm_rank( comm, &myRank );
rv = HPCC_TestDGEMM( params, 0 == myRank ? 1 : 0, &localGflops, &n, &failure );
params->DGEMM_N = n;
MPI_Reduce( &rv, &errCount, 1, MPI_INT, MPI_SUM, 0, comm );
MPI_Allreduce( &failure, &failureAll, 1, MPI_INT, MPI_MAX, comm );
if (failureAll) params->Failure = 1;
MPI_Reduce( &localGflops, &minGflops, 1, MPI_DOUBLE, MPI_MIN, 0, comm );
MPI_Reduce( &localGflops, &avgGflops, 1, MPI_DOUBLE, MPI_SUM, 0, comm );
MPI_Reduce( &localGflops, &maxGflops, 1, MPI_DOUBLE, MPI_MAX, 0, comm );
avgGflops /= (double)commSize;
MPI_Bcast( &avgGflops, 1, MPI_DOUBLE, 0, comm ); params->StarDGEMMGflops = avgGflops;
BEGIN_IO( myRank, params->outFname, outputFile);
fprintf( outputFile, "Node(s) with error %d\n", errCount );
fprintf( outputFile, "Minimum Gflop/s %.6f\n", minGflops );
fprintf( outputFile, "Average Gflop/s %.6f\n", avgGflops );
fprintf( outputFile, "Maximum Gflop/s %.6f\n", maxGflops );
END_IO( myRank, outputFile );
return 0;
}
int
HPCC_SingleDGEMM(HPCC_Params *params) {
int myRank, commSize;
int rv, errCount, rank, failure = 0;
double localGflops;
int n;
double scl = 1.0 / RAND_MAX;
FILE *outputFile;
MPI_Comm comm = MPI_COMM_WORLD;
localGflops = 0.0;
MPI_Comm_size( comm, &commSize );
MPI_Comm_rank( comm, &myRank );
srand(time(NULL));
scl *= (double)commSize;
/* select a node at random, but not node 0 (unless there is just one node) */
if (1 == commSize)
rank = 0;
else
for (rank = 0; ; rank = (int)(scl * rand())) {
if (rank > 0 && rank < commSize) break;
}
MPI_Bcast( &rank, 1, MPI_INT, 0, comm ); /* broadcast the rank selected on node 0 */
if (myRank == rank) /* if this node has been selected */
rv = HPCC_TestDGEMM( params, 0 == myRank ? 1 : 0, &localGflops, &n, &failure );
MPI_Bcast( &rv, 1, MPI_INT, rank, comm ); /* broadcast error code */
MPI_Bcast( &failure, 1, MPI_INT, rank, comm ); /* broadcast failure indication */
errCount = rv;
if (failure) params->Failure = 1;
/* broadcast result */
MPI_Bcast( &localGflops, 1, MPI_DOUBLE, rank, comm );
params->SingleDGEMMGflops = localGflops;
BEGIN_IO( myRank, params->outFname, outputFile);
fprintf( outputFile, "Node(s) with error %d\n", errCount );
fprintf( outputFile, "Node selected %d\n", rank );
fprintf( outputFile, "Single DGEMM Gflop/s %.6f\n", localGflops );
END_IO( myRank, outputFile );
return 0;
}
|
