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
| /* Test dot-product function */
|
| #include <stdio.h>
| #include <stdlib.h>
| #include <memory.h>
| #include <math.h>
| #include "config.h"
| #ifdef HAVE_GETOPT_H
| #include <getopt.h>
| #endif
| #include "fec.h"
|
| #if HAVE_GETOPT_LONG
| struct option Options[] = {
| {"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'},
| {"trials",0,NULL,'n'},
| {NULL},
| };
| #endif
|
| int main(int argc,char *argv[]){
| short coeffs[512];
| short input[2048];
| int trials=1000,d;
| int errors = 0;
|
| #if HAVE_GETOPT_LONG
| while((d = getopt_long(argc,argv,"apmstn:",Options,NULL)) != EOF){
| #else
| while((d = getopt(argc,argv,"apmstn:")) != 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 'n':
| trials = atoi(optarg);
| break;
| }
| }
|
| while(trials--){
| long port_result;
| long simd_result;
| int ntaps;
| int i;
| int csum = 0;
| int offset;
| void *dp_simd,*dp_port;
|
| /* Generate set of coefficients
| * limit sum of absolute values to 32767 to avoid overflow
| */
| memset(coeffs,0,sizeof(coeffs));
| for(i=0;i<512;i++){
| double gv;
|
| gv = normal_rand(0.,100.);
| if(csum + fabs(gv) > 32767)
| break;
| coeffs[i] = gv;
| csum += fabs(gv);
| }
| ntaps = i;
|
| /* Compare results to portable C version for a bunch of random data buffers and offsets */
| dp_simd = initdp(coeffs,ntaps);
| dp_port = initdp_port(coeffs,ntaps);
|
| for(i=0;i<2048;i++)
| input[i] = random();
|
| offset = random() & 511;
|
| simd_result = dotprod(dp_simd,input+offset);
| port_result = dotprod_port(dp_port,input+offset);
| if(simd_result != port_result){
| errors++;
| }
| }
| printf("dtest: %d errors\n",errors);
| exit(0);
| }
|
|
