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
| #include <stdio.h>
| #include <stdlib.h>
| #include <memory.h>
| #include <time.h>
| #include "config.h"
| #ifdef HAVE_GETOPT_H
| #include <getopt.h>
| #endif
| #include "fec.h"
|
| #if HAVE_GETOPT_LONG
| struct option Options[] = {
| {"frame-length",1,NULL,'l'},
| {"frame-count",1,NULL,'n'},
| {"verbose",0,NULL,'v'},
| {"force-altivec",0,NULL,'a'},
| {"force-port",0,NULL,'p'},
| {"force-mmx",0,NULL,'m'},
| {"force-sse",0,NULL,'s'},
| {"force-sse2",0,NULL,'t'},
| {NULL},
| };
| #endif
|
| int Verbose = 0;
|
| int main(int argc,char *argv[]){
| signed short *buf;
| int i,d,trial,trials=10000;
| int bufsize = 2048;
| long long port_sum,simd_sum;
| time_t t;
| int timetrials=0;
|
| find_cpu_mode();
| time(&t);
| srandom(t);
|
| #if HAVE_GETOPT_LONG
| while((d = getopt_long(argc,argv,"vapmstl:n:T",Options,NULL)) != EOF){
| #else
| while((d = getopt(argc,argv,"vapmstl:n:T")) != EOF){
| #endif
| switch(d){
| case 'a':
| Cpu_mode = ALTIVEC;
| break;
| case 'p':
| Cpu_mode = PORT;
| break;
| case 'm':
| Cpu_mode = MMX;
| break;
| case 's':
| Cpu_mode = SSE;
| break;
| case 't':
| Cpu_mode = SSE2;
| break;
| case 'l':
| bufsize = atoi(optarg);
| break;
| case 'n':
| trials = atoi(optarg);
| break;
| case 'v':
| Verbose++;
| break;
| case 'T':
| timetrials++;
| break;
| }
| }
|
| buf = (signed short *)calloc(bufsize,sizeof(signed short));
| if(timetrials){
| for(trial=0;trial<trials;trial++){
| (void)sumsq(buf,bufsize);
| }
| } else {
| for(trial=0;trial<trials;trial++){
| int length,offset;
|
| offset = random() & 7;
| length = (random() % bufsize) - offset;
| if(length <= 0)
| continue;
| for(i=0;i<bufsize;i++)
| buf[i] = random();
|
| port_sum = sumsq_port(buf+offset,length);
| simd_sum = sumsq(buf+offset,length);
| if(port_sum != simd_sum){
| printf("offset %d len %d port_sum = %lld simd_sum = %lld ",offset,length,port_sum,simd_sum);
|
| printf("ERROR! diff = %lld\n",simd_sum-port_sum);
| }
| }
| }
| exit(0);
| }
|
|
