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
|
/*
* PDTIME.C - PDBLib performance test
*
* Source Version: 9.0
* Software Release #92-0043
*
*/
#include "cpyright.h"
#include "pdb.h"
#if 1
#define PDB_TIME
#endif
/*--------------------------------------------------------------------------*/
/*--------------------------------------------------------------------------*/
/* MAIN - start it all here */
/* To find the peak rate for a given array length, start with a number of
* iterations such that the elapsed time is around 1 second and then vary it.
* For shorter elapsed times, timing resolution effects cause noisy results.
* For longer elapsed times, file length effects(?) cause byte rates to drop.
*/
int main(argc, argv)
int argc;
char **argv;
{int i, len, n, is, ia, size, bytes;
char filename[20];
double dt, *f;
#ifdef PDB_TIME
char s[20], **names;
PDBfile *pdbf;
#else
FILE *binf;
#endif
SIGNAL(SIGINT, SIG_DFL);
setvbuf(stdout, NULL, _IONBF, 0);
if (argc < 3)
{printf("\nUsage: pdtime <array-len> <n-iter> [<std> <align>]\n\n");
printf("<array-len> is the number of elements in the array\n");
printf("<n-iter> is the number of times the array is written\n");
printf("<std> is the pdb file data standard\n");
printf("<align> is the pdb file data alignment\n\n");
exit(0);};
len = atoi(argv[1]);
n = atoi(argv[2]);
if (argc == 5)
{is = atoi(argv[3]);
ia = atoi(argv[4]);
REQ_STANDARD = PD_std_standards[is - 1];
REQ_ALIGNMENT = PD_std_alignments[ia - 1];
printf("Data Standard = %d\nData Alignment = %d\n", is, ia);};
#ifdef PDB_TIME
names = MAKE_N(char *, n);
for (i = 0; i < n; i++)
{sprintf(s, "f%d-%d(%d)", len, i, len);
names[i] = SC_strsave(s);};
SC_zero_space(0);
strcpy(filename, "time.pdb");
if ((pdbf = PD_open(filename, "w")) == NULL)
{printf("Time couldn't create PDB file %s\r\n", filename);
exit(1);};
#else
strcpy(filename, "time.bin");
if ((binf = fopen(filename, BINARY_MODE_W)) == NULL)
{printf("Time couldn't create binary file %s\r\n", filename);
exit(1);};
#endif
size = sizeof(double);
bytes = len*n*size;
printf("\nArray Size Iterations Time Rate\n");
printf("(elements) (sec) (bytes/sec)\n");
if ((f = MAKE_N(double, len)) == NULL)
{printf("\nCAN'T ALLOCATE %d DOUBLES\n\n", len);
exit(1);};
for (i = 0; i < len; i++)
f[i] = i;
dt = SC_wall_clock_time();
#ifdef PDB_TIME
for (i = 0; i < n; i++)
PD_write(pdbf, names[i], "double", f);
#else
for (i = 0; i < n; i++)
fwrite(f, (size_t) size, (size_t) len, binf);
#endif
dt = SC_wall_clock_time() - dt;
printf("%8d %8d %11.3e %11.3e\n\n", len, n,
dt, ((double) bytes)/dt);
SFREE(f);
#ifdef PDB_TIME
PD_close(pdbf);
#else
fclose(binf);
#endif
exit(0);}
/*--------------------------------------------------------------------------*/
/*--------------------------------------------------------------------------*/
|
