#include <stdio.h>
#include <time.h>
#include <omp.h>

int main() {
  long long int i, n=50000;
  float a=2.0;
  float x[n];
  float y[n];
  double start_time, end_time;

#pragma omp parallel for private(i)
  for (i = 0; i < n; i++){
    x[i] = 1.0;
    y[i] = 2.0;
  }

  start_time = omp_get_wtime();
//#pragma omp parallel for schedule(static,100) private(i)
//#pragma omp parallel for private(i)
#pragma omp parallel for shared(n,a) private(i)
  for (i = 0; i < n; i++){
    y[i] = a*x[i] + y[i];
  }
  end_time = omp_get_wtime();

  printf ("SAXPY Time: %f\n", end_time - start_time);
  printf ("%f\n", y[0]);
}