85{
86 int i, packed = 0, sse = 0;
87 float a[4] = { 1.0, 2.0, 3.0, 4.0 };
88 float b[4] = { 2.0, 3.0, 4.0, 5.0 };
89 float c[4] = { 0.0, 0.0, 0.0, 0.0 };
90 double d[4] = { 1.0, 2.0, 3.0, 4.0 };
91 double e[4] = { 2.0, 3.0, 4.0, 5.0 };
92 double f[4] = { 0.0, 0.0, 0.0, 0.0 };
93
94 if ( argc != 3 ) {
95 bail:
96 printf( "Usage %s: <packed|unpacked> <sse|sse2>\n", argv[0] );
97 exit( 1 );
98 }
99 if ( strcasecmp( argv[1], "packed" ) == 0 )
100 packed = 1;
101 else if ( strcasecmp( argv[1], "unpacked" ) == 0 )
102 packed = 0;
103 else
104 goto bail;
105 if ( strcasecmp( argv[2], "sse" ) == 0 )
106 sse = 1;
107 else if ( strcasecmp( argv[2], "sse2" ) == 0 )
108 sse = 0;
109 else
110 goto bail;
111
112#if 0
113 if ( ( sse ) &&
114 ( system( "cat /proc/cpuinfo | grep sse > /dev/null" ) != 0 ) ) {
115 printf( "This processor does not have SSE.\n" );
116 exit( 1 );
117 }
118 if ( ( sse == 0 ) &&
119 ( system( "cat /proc/cpuinfo | grep sse2 > /dev/null" ) != 0 ) ) {
120 printf( "This processor does not have SSE2.\n" );
121 exit( 1 );
122 }
123#endif
124
125 printf(
"Vector 1: %f %f %f %f\n",
a[0],
a[1],
a[2],
a[3] );
126 printf(
"Vector 2: %f %f %f %f\n\n",
b[0],
b[1],
b[2],
b[3] );
127
128 if ( ( packed == 0 ) && ( sse == 1 ) ) {
131 }
132 printf(
"%d SSE Unpacked Adds: Result %f\n",
NUMBER,
c[0] );
133
136 }
137 printf(
"%d SSE Unpacked Muls: Result %f\n",
NUMBER,
c[0] );
138 }
139 if ( ( packed == 1 ) && ( sse == 1 ) ) {
142 }
143 printf(
"%d SSE Packed Adds: Result %f %f %f %f\n",
NUMBER,
c[0],
c[1],
147 }
148 printf(
"%d SSE Packed Muls: Result %f %f %f %f\n",
NUMBER,
c[0],
c[1],
150 }
151
152 if ( ( packed == 0 ) && ( sse == 0 ) ) {
155 }
156 printf(
"%d SSE2 Unpacked Adds: Result %f\n",
NUMBER,
c[0] );
157
160 }
161 printf(
"%d SSE2 Unpacked Muls: Result %f\n",
NUMBER,
c[0] );
162 }
163 if ( ( packed == 1 ) && ( sse == 0 ) ) {
166 }
167 printf(
"%d SSE2 Packed Adds: Result %f\n",
NUMBER,
c[0] );
168
171 }
172 printf(
"%d SSE2 Packed Muls: Result %f\n",
NUMBER,
c[0] );
173 }
174
175
176 exit( 0 );
177}
static double a[MATRIX_SIZE][MATRIX_SIZE]
static double b[MATRIX_SIZE][MATRIX_SIZE]
static double c[MATRIX_SIZE][MATRIX_SIZE]
void inline_packed_sse_add(float *aa, float *bb, float *cc)
void inline_unpacked_sse2_mul(double *aa, double *bb, double *cc)
void inline_packed_sse2_add(double *aa, double *bb, double *cc)
void inline_unpacked_sse2_add(double *aa, double *bb, double *cc)
void inline_unpacked_sse_mul(float *aa, float *bb, float *cc)
void inline_packed_sse_mul(float *aa, float *bb, float *cc)
void inline_packed_sse2_mul(double *aa, double *bb, double *cc)
void inline_unpacked_sse_add(float *aa, float *bb, float *cc)