在C/C++代码中使用SSE等指令集的指令(5)SSE进行加法运算简单的性能测试

下面是一个简单的测试SSE指令性能的程序，可以看到明显的性能提升。

（说明：程序中的timing.h使用的是/article/8516656.html中介绍的时间间隔获取方法）

[cpp]
view plaincopyprint?

#define WIN
#include "timing.h"
#include <intrin.h>
#include <stdlib.h>
#include <math.h>

#define N 4*100000 // 注意：必须是4的倍数，否则使用SSE指令计算，要进行一些处理，从而保证正确。
_MM_ALIGN16 float op1
;
_MM_ALIGN16 float op2
;
_MM_ALIGN16 float result1
;
_MM_ALIGN16 float result2
;

void init()
{
for(int i = 0;i < N; i++)
{
op1[i] = (float)rand()/(float)RAND_MAX;
op2[i] = (float)rand()/(float)RAND_MAX;
}
}

void checkResult(int debug)
{
bool isSame = true;
for(int i = 0;i < N; i++)
{
if (debug)
{
printf("%lf %lf\n", result1[i], result2[i]);
}
else
{
if (fabs(result1[i] - result2[i]) > 0.000001)
{
isSame = false;
break;
}
}
}
if (!debug) {
if (isSame)
printf("Result is Same\n");
else
printf("Result is not same\n");
}
}

void add1()
{
for(int i = 0; i < N;i++)
result1[i] = op1[i] + op2[i];
}

void add2()
{
__m128 a;
__m128 b;
__m128 c;

for(int i = 0; i < N;i = i + 4)
{
// Load
a = _mm_load_ps(op1 + i);
b = _mm_load_ps(op2 + i);

c = _mm_add_ps(a, b); // c = a + b

_mm_store_ps(result2 + i, c);
}
}

int main(int argc, char* argv[])
{
init();
srand((unsigned int)time(NULL));

printf("Add a vector:\n");
startTiming();
add1();
stopWithPrintTiming();

printf("\n");
printf("Add a vector with SSE instructions:\n");
startTiming();
add2();
stopWithPrintTiming();

printf("\n");
checkResult(0);

return 0;
}

#define WIN
#include "timing.h"
#include <intrin.h>
#include <stdlib.h>
#include <math.h>

#define N 4*100000		// 注意：必须是4的倍数，否则使用SSE指令计算，要进行一些处理，从而保证正确。
_MM_ALIGN16 float op1
;
_MM_ALIGN16 float op2
;
_MM_ALIGN16 float result1
;
_MM_ALIGN16 float result2
;

void init()
{
for(int i = 0;i < N; i++)
{
op1[i] = (float)rand()/(float)RAND_MAX;
op2[i] = (float)rand()/(float)RAND_MAX;
}
}

void checkResult(int debug)
{
bool isSame = true;
for(int i = 0;i < N; i++)
{
if (debug)
{
printf("%lf		%lf\n", result1[i], result2[i]);
}
else
{
if (fabs(result1[i] - result2[i]) > 0.000001)
{
isSame = false;
break;
}
}
}
if (!debug) {
if (isSame)
printf("Result is Same\n");
else
printf("Result is not same\n");
}
}

void add1()
{
for(int i = 0; i < N;i++)
result1[i] = op1[i] + op2[i];
}

void add2()
{
__m128  a;
__m128  b;
__m128  c;

for(int i = 0; i < N;i = i + 4)
{
// Load
a = _mm_load_ps(op1 + i);
b = _mm_load_ps(op2 + i);

c = _mm_add_ps(a, b);	// c = a + b

_mm_store_ps(result2 + i, c);
}
}

int main(int argc, char* argv[])
{
init();
srand((unsigned int)time(NULL));

printf("Add a vector:\n");
startTiming();
add1();
stopWithPrintTiming();

printf("\n");
printf("Add a vector with SSE instructions:\n");
startTiming();
add2();
stopWithPrintTiming();

printf("\n");
checkResult(0);

return 0;
}

转自：http://blog.csdn.net/gengshenghong/article/details/7011373