SSE2 - 16 -байтовая выровненная динамическое распределение памяти
-
18-09-2019 - |
Вопрос
РЕДАКТИРОВАТЬ:
Это продолжение Ошибка компилятора SSE2
Это настоящий ошибка, которую я испытывал раньше и воспроизводил ниже, изменив _mm_malloc
заявление как Майкл Берр предложенный:
Неправомерное исключение на 0x00415116 в So.exe: 0xc0000005: Расположение считывания нарушений доступа 0xffffffff.
В линии label: movdqa xmm0, xmmword ptr [t1+eax]
Я пытаюсь динамически распределить t1
а также t2
а также Согласно этому уроку, Я использовал _mm_malloc
:
#include <emmintrin.h>
int main(int argc, char* argv[])
{
int *t1, *t2;
const int n = 100000;
t1 = (int*)_mm_malloc(n*sizeof(int),16);
t2 = (int*)_mm_malloc(n*sizeof(int),16);
__m128i mul1, mul2;
for (int j = 0; j < n; j++)
{
t1[j] = j;
t2[j] = (j+1);
} // set temporary variables to random values
_asm
{
mov eax, 0
label: movdqa xmm0, xmmword ptr [t1+eax]
movdqa xmm1, xmmword ptr [t2+eax]
pmuludq xmm0, xmm1
movdqa mul1, xmm0
movdqa xmm0, xmmword ptr [t1+eax]
pshufd xmm0, xmm0, 05fh
pshufd xmm1, xmm1, 05fh
pmuludq xmm0, xmm1
movdqa mul2, xmm0
add eax, 16
cmp eax, 100000
jnge label
}
_mm_free(t1);
_mm_free(t2);
return 0;
}
Решение
Я думаю, что 2 -я проблема заключается в том, что вы читаете в смещении с переменной указателя (а не смещение от того, на что указывает указатель).
Изменять:
label: movdqa xmm0, xmmword ptr [t1+eax]
Что -то вроде:
mov ebx, [t1]
label: movdqa xmm0, xmmword ptr [ebx+eax]
И аналогично вашим доступу через указатель T2.
Это может быть еще лучше (хотя у меня не было возможности проверить это, так что это может даже не сработать):
_asm
{
mov eax, [t1]
mov ebx, [t1]
lea ecx, [eax + (100000*4)]
label: movdqa xmm0, xmmword ptr [eax]
movdqa xmm1, xmmword ptr [ebx]
pmuludq xmm0, xmm1
movdqa mul1, xmm0
movdqa xmm0, xmmword ptr [eax]
pshufd xmm0, xmm0, 05fh
pshufd xmm1, xmm1, 05fh
pmuludq xmm0, xmm1
movdqa mul2, xmm0
add eax, 16
add ebx, 16
cmp eax, ecx
jnge label
}
Другие советы
Вы не распределяете достаточно памяти:
t1 = (int*)_mm_malloc(n * sizeof( int),16);
t2 = (int*)_mm_malloc(n * sizeof( int),16);
Возможно:
t1 = (int*)_mm_malloc(n*sizeof(int),16);