Chinaunix

标题: g++ 编译ｓｓｅ出错 [打印本页]

作者: yyfn风辰 时间: 2010-06-27 20:05
标题: g++ 编译ｓｓｅ出错
代码如下：
#include <iostream>
#include <math.h>
#include <stdio.h>
#include <xmmintrin.h>

using namespace std;

typedef union{
__m128 m;
float f[4];
} sseTemp;

float dotWithSSE(float *a, float *b, const int len){
int step = len/4;
__m128* one = (__m128*)a;
__m128* two = (__m128*)b;
if(NULL == one || NULL == two){
printf("%d\n", __LINE__);
exit(1);
}

__m128 result = _mm_setzero_ps();
__m128 temp = _mm_set_ps(0.1f, 0.1f, 0.1f, 0.1f);

for(int i = 0; i < step -1; i++){
temp = _mm_add_ps(*one, *two);
temp = _mm_mul_ps(temp, temp);
result = _mm_add_ps(result, temp);
}

sseTemp s;
s.m = result;

float r = 0.0f;

r += s.f[0] + s.f[1] + s.f[2] + s.f[3];

if(0 != (len & 3)){
for(int i = 4*step; i < len; i++){
r += (a[i] - b[i])*(a[i] - b[i]);
}
}

return r;
}

float dotWithoutSSE(float* __restrict__ a, float* __restrict__ b, const int len){
float result = 0.0f, temp;

for(int i = 0; i < len; i++){
temp = (a[i] - b[i]);
result += temp*temp;
}

return result;
}

int runTest(const int len){
float *a, *b;
a = (float*) malloc(sizeof(float)*len);
b = (float*) malloc(sizeof(float)*len);
if(NULL == a || NULL == b){
printf("malloc fail\n");
exit(1);
}

for(int i = 0; i < len; i++){
a[i] = 1.0f*rand()/RAND_MAX;
b[i] = 1.0f*rand()/RAND_MAX;
}

clock_t start = clock();
float sse = dotWithSSE(a, b, len);
clock_t end = clock();
clock_t timeSSE = end - start;

start = clock();
float withoutSSE = dotWithoutSSE(a, b, len);
end = clock();
clock_t timeWithoutSSE = end - start;

if(fabs(sse - withoutSSE) < 0.00001){
printf("test pass\n");
}
printf("withSSE = %f, wihtoutSSE = %f", sse, withoutSSE);

free(a);
free(b);
}

int main(){

runTest(1000);
return 0;
}

使用ｇｃｃ测试时，代码 temp = _mm_add_ps(*one, *two);越界，这另我百思不得其解，先谢谢了！

作者: 变异老鼠 时间: 2010-06-27 22:40
好像是地址没对齐，SSE 要求 16-Byte 对齐的。

作者: yyfn风辰 时间: 2010-06-27 22:57
回复 2# 变异老鼠
谢谢你！

那在ｇｃｃ中怎么将ｍａｌｌｏｃ分配的空间对齐到１６字节啊？我试过__attribute__(aligned(16))不行啊！

作者: hellioncu 时间: 2010-06-28 09:12

好像是地址没对齐，SSE 要求 16-Byte 对齐的。
变异老鼠发表于 2010-06-27 22:40

这个运行期才能判断出来的。LZ的问题是编译错误

作者: 变异老鼠 时间: 2010-06-28 10:22
回复 3# yyfn风辰

试试这个？

NAME
posix_memalign, memalign, valloc - Allocate aligned memory
SYNOPSIS
#include <stdlib.h>
int posix_memalign(void **memptr, size_t alignment, size_t size);
#include <malloc.h>
void *valloc(size_t size);
void *memalign(size_t boundary, size_t size);
Feature Test Macro Requirements for glibc (see feature_test_macros(7)):
posix_memalign(): _POSIX_C_SOURCE >= 200112L || _XOPEN_SOURCE >= 600

复制代码

作者: yyfn风辰 时间: 2010-06-28 19:51
谢谢！

作者: yyfn风辰 时间: 2010-07-05 22:23
回复 5# 变异老鼠
非常感谢你的回复，我测试了，完全可用。另外，我查到valloc分配的空间是按页对齐的，这个应当会多占空间，这个函数有什么好处呢？

欢迎光临 Chinaunix (http://bbs.chinaunix.net/)