vecacc-aligned.c

Go to the documentation of this file.

#include <stdio.h>
#include <stdlib.h>
#include <unistd.h>
//#include "tools.h"
 
void vecacc(int n, float* src1, float* src2, float* result)
{
    int i;
#pragma vector aligned
#pragma ivdep
    for(i=0;i<n;i++)
        result[i]=src1[i]*src2[i];
}
 
int main(int argc, char **argv)
{
        int n;
        float *src1, *src2, *result;
        n = atoi(argv[1]);
 
#define xmalloc(p)                                                      \
        do {                                                            \
                if (posix_memalign((void **) &p, 32, n * sizeof(float))) \
                        return 3;                                       \
        } while(0);
#define vmalloc(p)\
        do {\
                p=(float*)malloc(n*sizeof(float));\
        } while(0)
 
        xmalloc(src1);
        xmalloc(src2);
        xmalloc(result);
 
        init_data_file(argv[2]);
        init_data_float(src1, n);
        init_data_float(src2, n);
        close_data_file();
 
        vecacc(n, src1, src2, result);
 
        print_array_float("result", result, n);
 
        free(src1);
        free(src2);
        free(result);
 
        return 0;
}