SSE инструкции + ликбез

[Kukurudza]

нашел такую штуку в википедии:

float a[4] = { 300.0, 4.0, 4.0, 12.0 };
float b[4] = {   1.5, 2.5, 3.5,  4.5 };

__asm {
movups xmm0, a   ; // поместить 4 переменные с плавающей точкой из a в регистр xmm0
movups xmm1, b   ; // поместить 4 переменные с плавающей точкой из b в регистр xmm1

mulps xmm1, xmm0 ; // перемножить пакеты плавающих точек: xmm1=xmm1*xmm0
                 ; // xmm10 = xmm10*xmm00
                 ; // xmm11 = xmm11*xmm01
                 ; // xmm12 = xmm12*xmm02
                 ; // xmm13 = xmm13*xmm03

movups a, xmm1   ; // выгрузить результаты из регистра xmm1 по адресам a
};

так работает. но хотелось бы, чтобы массивы были размером скажем по 10000 чисел. нужно это как то делать через указатели и смещение. покажите как?
ну а идеале хотелось бы такую штуку проворачивать с STL-векторами.

как это можно сделать на простых массивах через указатели и смещение и как это сделать для векторов?

[f.hump]

__declspec(align(16)) float a[10000] = { 300.0, 4.0, 4.0, 12.0, ..... };
__declspec(align(16)) float b[10000] = {   1.5, 2.5, 3.5,  4.5, ..... };

__asm {
LEA esi, [a]
LEA edi, [b]

MOV ecx, 625

align 16
mloop:

MOVAPS xmm0, [esi]
MOVAPS xmm1, [esi+16]
MOVAPS xmm2, [esi+32]
MOVAPS xmm3, [esi+48]

MULPS xmm0, [edi]
MULPS xmm1, [edi+16]
MULPS xmm2, [edi+32]
MULPS xmm3, [edi+48]

MOVAPS [esi], xmm0
MOVAPS [esi+16], xmm1
MOVAPS [esi+32], xmm2
MOVAPS [esi+48], xmm3


ADD esi, 64
ADD edi, 64

LOOP mloop

};

enjoy


ps. думаю, что утомлять себя асм мурзилкой не нужно, потому что копилятор оптимизирует очень хорошо (особенно в случае использования STL). принимая во внимание, что MS компилятор не дружит с ассемблерными вставками в 64-х битном режиме, то смысл от подобных упражнений только общеобразовательный.

[Kukurudza]

FUN
на самом деле, интересно почему он не любит АСМ в 64.
ну я прогнал эту штуку как написано выше, и как делают по обыкновенному:

		for ( size_t j = 0; j < N; ++j )
		{
			a[j] = a[j] * b[j];
		};

прогнал в цикле например 100.000 раз. прирост скорости ошеломил, аж в 500!!! раз быстрее. для меня, напрмер, где я усиленно юзаю векторное перемножение, среди прочих операций, это не плохой плюс.

[f.hump]

прирост скорости ошеломил, аж в 500!!!

для какого N?

[Kukurudza]

static const size_t N = 1000;
static const size_t M = 100000;
	for (register size_t i = 0; i < M; ++i)
	{
		for (register size_t j = 0; j < N; ++j )
		{
			a[j] = a[j] * b[j];
		};
	};

ну и то же самое для SSE