#include <stdlib.h>
#include <xmmintrin.h>  

#ifndef N
#error N is not defined
#endif

void compute (float *Y, float *X1, float *X2){
  int i;
  for (i=0;i<N/4;i++)
    ((__m128 *)Y)[i] = _mm_mul_ps(((__m128 *)X1)[i], ((__m128 *)X2)[i]);
}

int main(){
  int i;
  float* Y = (float*)_mm_malloc(N*sizeof(float),16);
  float* X1 = (float*)_mm_malloc(N*sizeof(float),16);
  float* X2 = (float*)_mm_malloc(N*sizeof(float),16);

  for(i=0;i<N;i++){
    Y[i]=X1[i]=X2[i]=i;
  }

  compute(Y, X1, X2);
  return 0;
}