void vecadd(int n, float* restrict a, float* restrict b, float* restrict c)
{
#pragma acc kernels deviceptr(a,b,c)
  {
    for(int i=0; i<n; i++)
    {
      c[i] = a[i] + b[i];
    }
  }
}
void set_value(int n, float * restrict aa, float * restrict bb)
{
#pragma acc kernels deviceptr(aa,bb)
  {
    for(int i=0; i<n; i++)
    {
      aa[i] = (float)i ;
      bb[i] = (float)i*2.0f;
    }
  }
}