SSE. Подение программы

[Bummmm]

Привет ребята!
Такое дело, есть задание:

Обработать массивы из 8 элементов по следующему выражению:

F[i]=A[i]-B[i]*C[i]-D[i] , i=1...8;

Используются следующие массивы:
A, B и С – 8 разрядные целые знаковые числа (_int8);
D – 16 разрядные целые знаковые числа (_int16).
Полученный результат отобразить на форме с использованием соответствующих элементов. При распаковке знаковых чисел совместно с командами распаковки использовать команды сравнения (сравнивать с нулём перед распаковкой).

Всё бы ничего, но как только я пытаюсь сделать вычисления, программа падает (в 1 из 5 случаем, обрабатывает 3 элемента и зависает):
https://www.dropbox.com/s/0spd9ctqqe57d2g/SSEDeb.png

Я не могу понять, в чем дело (компилирую в Visula Studio 2010):
Код (с участком показанным на скриншоте):

#include "stdafx.h"
#include "SSE.h"
#include <iostream>
#include <malloc.h>
using namespace std;

double* SSESolve(float* a, float* b, float* c, double* d)
{
	double* resultPointer = (double*)malloc(8 * sizeof(double));
	float arr_a[8];
	float arr_b[8];
	float arr_c[8];
	double arr_d[8];
	float preresult[8];	
	double result[8];
	memcpy((void*)(&arr_a), (void*)a, 8*sizeof(float));	
	memcpy((void*)(&arr_b), (void*)b, 8*sizeof(float));
	memcpy((void*)(&arr_c), (void*)c, 8*sizeof(float));
	memcpy((void*)(&arr_d), (void*)d, 8*sizeof(double));
	memset((void*)resultPointer, 0, 8*sizeof(double));
	__asm
	{
		/*  F[i]=A[i] - B[i] * C[i] + D[i] , i=1...8; */

		/*
		Команды для чисел с плавающей точкой
Команды пересылки
Скалярные типы – MOVSS
Упакованные типы – MOVAPS, MOVUPS, MOVLPS, MOVHPS, MOVLHPS, MOVHLPS
Арифметические команды
Скалярные типы – ADDSS, SUBSS, MULSS, DIVSS, RCPSS, SQRTSS, MAXSS, MINSS, RSQRTSS
Упакованные типы – ADDPS, SUBPS, MULPS, DIVPS, RCPPS, SQRTPS, MAXPS, MINPS, RSQRTPS
Команды сравнения
Скалярные типы – CMPSS, COMISS, UCOMISS
Упакованные типы – CMPPS
Перемешивание и распаковка
Упакованные типы – SHUFPS, UNPCKHPS, UNPCKLPS
Команды для преобразования типов
Скалярные типы – CVTSI2SS, CVTSS2SI, CVTTSS2SI
Упакованные типы – CVTPI2PS, CVTPS2PI, CVTTPS2PI
Битовые логические операции
Упакованные типы – ANDPS, ORPS, XORPS, ANDNPS
Команды для целых чисел
Арифметические команды
PMULHUW, PSADBW, PAVGB, PAVGW, PMAXUB, PMINUB, PMAXSW, PMINSW
Команды пересылки - PEXTRW, PINSRW
Другие - PMOVMSKB, PSHUFW
Другие команды
Работа с регистром - MXCSR, LDMXCSR, STMXCSR
Управление кэшем и памятью
MOVNTQ, MOVNTPS, MASKMOVQ, PREFETCH0, PREFETCH1, PREFETCH2, PREFETCHNTA, SFENCE
*/
		movups xmm0, arr_b;
		movups xmm1, arr_c;
		mulps xmm0, xmm4;
		movups preresult, xmm0;
		movups xmm0, arr_b+16;
		movups xmm1, arr_c+16;
		mulps xmm0, xmm1;
		movups preresult+16, xmm0;
			
		movups xmm0, arr_a;
		movups xmm1, preresult;
		subps xmm0, xmm1;
		movups preresult, xmm0;
		movups xmm0, arr_a+16;
		movups xmm1, preresult+16;
		subps xmm0, xmm1;
		movups preresult+16, xmm0;
		
		movups xmm0, preresult;
		movaps xmm1, xmm0;
		shufpd xmm1, xmm1, 1; /*<-------------------------Вот здесь */
		cvtps2pd xmm0, xmm0;
		cvtps2pd xmm1, xmm1;
		movups xmm2, preresult+16
		movaps xmm3, xmm2;
		shufpd xmm3, xmm3, 1;
		cvtps2pd xmm2, xmm2
		cvtps2pd xmm3, xmm3;
		
		movups xmm4, arr_d;		
		movups xmm5, arr_d+16;	
		movups xmm6, arr_d+32;	
		movups xmm7, arr_d+48;	
			
		subpd xmm0, xmm4;
		subpd xmm1, xmm5;
		subpd xmm2, xmm6;
		subpd xmm3, xmm7;

		movups result, xmm0;
		movups result+16, xmm1;
		movups result+32, xmm2;
		movups result+48, xmm3;
			
	  
	}
	memcpy((void*)resultPointer, (void*)(&result), 8 * sizeof(double));
	return resultPointer;
}

SSE.h

double* SSESolve(float* a, float* b, float* c, double* d);

Весь проект полностью (55.9 MB):
https://www.dropbox.com/sh/0b1f14yjyiebw1m/Am7KMEglj-

Ребята, помогите разобраться O_o

[Bummmm]

Проблему решил.
Заменил

 movaps xmm1, xmm0;
 shufpd xmm1, xmm1, 1;

на

PSHUFD xmm1, xmm0, 00001110b

[f.hump]

странно, конечно. что это за проц, который не дружит с SSE2?

Какой на твоей машине процессор стоит?

[f.hump]

Заменил
Код:

movaps xmm1, xmm0;
shufpd xmm1, xmm1, 1;

на
Код:

PSHUFD xmm1, xmm0, 00001110b

это, что одно и тоже???

если уж делать такую странную замену, то делать так

movaps xmm1, xmm0;
PSHUFD xmm1, xmm0, 04Eh;

[Bummmm]

Сообщение от f.hump

странно, конечно. что это за проц, который не дружит с SSE2?

Какой на твоей машине процессор стоит?

Машина старая, AMD Athlon XP 2200+

[asmprog]

не лишним будет проверку воткнуть перед использованием на наличие поддержки SSE/SSE2:

mov eax,01
cpuid
test edx,03000000h
jz _SSE_AND_SSE2_NOT_SUPPORTED

[Bummmm]

Сообщение от asmprog

не лишним будет проверку воткнуть перед использованием на наличие поддержки SSE/SSE2:

mov eax,01
cpuid
test edx,03000000h
jz _SSE_AND_SSE2_NOT_SUPPORTED

О, вот за это огромное спасибо!

[f.hump]

а, старый AMD, ясна.

кстати, посмотрев на код еще раз возникает вопрос зачем использовать SHUFPD вместо MOVHLPS?

вместо этого

movaps xmm1, xmm0;
shufpd xmm1, xmm1, 1;

поставь вот это

MOVHLPS xmm1, xmm0

[f.hump]

прочитал условие задачи

void SSESolve(int * result, const char * a, const char * b, const char * c, const short * d) {
	__asm {
		
		MOV edx, [b]
		MOV ecx, [c]
				
		PMOVSXBW xmm0, QWORD PTR [edx]
		PMOVSXBW xmm1, QWORD PTR [ecx]

		PMULLLW xmm0, xmm1
		
		MOV eax, [a]
		PMOVSXBD xmm2, DWORD PTR [eax]
		PMOVSXBD xmm3, DWORD PTR [eax+4]

		MOV eax, [d]

		PMOVSXWD xmm4, QWORD PTR [eax]
		PMOVSXWD xmm5, QWORD PTR [eax+8]

		PSUBD xmm2, xmm4
		PSUBD xmm3, xmm5


		PMOVSXWD xmm4, xmm0
		PSRLDQ xmm0, 8
		PMOVSXWD xmm5, xmm0
		
		PSUBD xmm2, xmm4
		PSUBD xmm3, xmm5

		MOV edx, [result]

		MOVDQU [edx], xmm2
		MOVDQU [edx+16], xmm3
	}
}

нужен SSE4.1

[Bummmm]

Сообщение от f.hump

а, старый AMD, ясна.

кстати, посмотрев на код еще раз возникает вопрос зачем использовать SHUFPD вместо MOVHLPS?

вместо этого

movaps xmm1, xmm0;
shufpd xmm1, xmm1, 1;

поставь вот это

MOVHLPS xmm1, xmm0

Понял, просто я ещё со всеми командами не разобрался