#include "axpby.h"

__attribute__ ((target(mic))) void axpby1(double *z, double a, double *x, double b, double *y, int n) {
	#pragma omp parallel
	{
		int i;
		#pragma omp for
		for (i = 0; i < n; i++)
			z[i] = a*x[i] + b*y[i];
	}
}

void axpby2(double *z, double a, double *x, double b, double *y, int n) {
	#pragma offload target(mic) in(x, y, z : length(0) REUSE)
	{
		#pragma omp parallel
		{
			int i;
			#pragma omp for
			for (i = 0; i < n; i++)
				z[i] = a*x[i] + b*y[i];
		}
	}
}