-
Notifications
You must be signed in to change notification settings - Fork 0
/
ap2.c
152 lines (119 loc) · 3.26 KB
/
ap2.c
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
#include <pmmintrin.h>
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <unistd.h>
// Cache L1: 512 líneas
// Cache L2: 4096 líneas
// Cache L3: 131072 líneas
#define CLS 64
// Double: 8 bytes
// Tamaño de línea: 64
// En cada línea 8 doubles
#define NUMD 8
// C=4
// F = L * 2
// C=8
// F = L
// C=20
// F = (floor(L / 5) * 2) + if(L % 5 <= 3, 1, 2)
// C=40
// F = ceil(L / 5)
// #define L 65536 // 512 768 2048 6144 65536 98304 262144 524288
// #define F 65536
// #define C 8 // 4 8 20 40
void start_counter();
double get_counter();
double mhz();
/* Initialize the cycle counter */
static unsigned cyc_hi = 0;
static unsigned cyc_lo = 0;
/* Set *hi and *lo to the high and low order bits of the cycle counter.
Implementation requires assembly code to use the rdtsc instruction. */
void access_counter(unsigned *hi, unsigned *lo)
{
asm("rdtsc; movl %%edx,%0; movl %%eax,%1" /* Read cycle counter */
: "=r"(*hi), "=r"(*lo) /* and move results to */
: /* No input */ /* the two outputs */
: "%edx", "%eax");
}
/* Record the current value of the cycle counter. */
void start_counter()
{
access_counter(&cyc_hi, &cyc_lo);
}
/* Return the number of cycles since the last call to start_counter. */
double get_counter()
{
unsigned ncyc_hi, ncyc_lo;
unsigned hi, lo, borrow;
double result;
/* Get cycle counter */
access_counter(&ncyc_hi, &ncyc_lo);
/* Do double precision subtraction */
lo = ncyc_lo - cyc_lo;
borrow = lo > ncyc_lo;
hi = ncyc_hi - cyc_hi - borrow;
result = (double)hi * (1 << 30) * 4 + lo;
if (result < 0)
{
fprintf(stderr, "Error: counter returns neg value: %.0f\n", result);
}
return result;
}
double mhz(int verbose, int sleeptime)
{
double rate;
start_counter();
sleep(sleeptime);
rate = get_counter() / (1e6 * sleeptime);
if (verbose)
printf("\n Processor clock rate = %.1f MHz\n", rate);
return rate;
}
int main(int argc, char ** argv)
{
if (argc != 3) {
printf("Uso: %s [F] [C]\n", argv[0]);
exit(1);
}
// int L = atoi(argv[1]);
int F = atoi(argv[1]);
int C = atoi(argv[2]);
srand(time(0));
double ck;
// Matriz de valores con reserva dinámica
double **M;
M = (double **) _mm_malloc(sizeof(double *) * F, CLS);
for (int i = 0; i < F; i++) *(M+i) = (double *) malloc(sizeof(double) * C);
for (int i = 0; i < F; i++) {
for (int j = 0; j < C; j++) {
M[i][j] = rand();
}
}
start_counter();
/* Poñer aquí o código a medir */
double red[10];
for (int i = 0; i < 10; i++) {
for (int j = 0; j < F; j++) {
int start = 0;
while (start < C) {
red[i] += M[j][start];
start += NUMD;
}
}
}
ck = get_counter() / 10;
for (int i = 0; i < 10; i++){
printf("%lf\t", red[i]);
}
printf("\n Clocks=%1.10lf \n", ck);
/* Esta rutina imprime a frecuencia de reloxo estimada coas rutinas start_counter/get_counter */
// mhz(1, 1);
for (int i = 0; i < F; i++) free(M[i]);
_mm_free(M);
FILE *f = fopen("res.txt", "w");
fprintf(f, "%lf", ck);
fclose(f);
return 0;
}