Умножение с AVX

Пожалуйста, я впервые использую AVX, и я пытаюсь выполнить простое умножение на числа с двойной точностью, но я не получаю правильные результаты.

Я получаю только первые 4 результата, а остальные — жаргон.

#include <immintrin.h>
#include <iostream>
#include <math.h>
#include <time.h>
using namespace std;

int main() {

double *a, *b;                      // data pointers
double *pA,*pB;                     // work pointer
__m256d rA_AVX, rB_AVX;     // variables for AVX

const int vector_size = 8;
a = (double*) _mm_malloc (vector_size*sizeof(double),64);
b = (double*) _mm_malloc (vector_size*sizeof(double),64);

for(int i = 0; i < vector_size; i++) {
a[i] = (rand() % 48);
b[i] = 0.0f;
cout << a[i] << endl;
}

for (int i = 0; i < vector_size; i += 8)
{
pA = a;
pB = b;
rA_AVX = _mm256_load_pd(pA);
rB_AVX = _mm256_mul_pd(rA_AVX,rA_AVX);
_mm256_store_pd(pB,rB_AVX);
pA += 8;
pB += 8;
}

for (int i=0; i<vector_size; i++){
cout << endl << b[i] << endl;
}
_mm_free(a);
_mm_free(b);

system("PAUSE");
return 0;
}