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

int main ( int argc, char *argv[] )
{
  int i, n = 10000000;
  double s, wtime, *x, *y;

  wtime = omp_get_wtime ( );

  s = 123.456;

  x = ( double * ) malloc ( n * sizeof ( double ) );
  y = ( double * ) malloc ( n * sizeof ( double ) );

# pragma omp parallel
  {
# pragma omp for
    for ( i = 0; i < n; i++ )
    {
      x[i] = ( double ) ( ( i + 1 ) % 17 );
      y[i] = ( double ) ( ( i + 1 ) % 31 );
    }

# pragma omp for
    for ( i = 0; i < n; i++ )
    {
      y[i] = y[i] + s * x[i];
    }
  }

  free ( x );
  free ( y );

  wtime = omp_get_wtime ( ) - wtime;
  printf ( "%g seconds.\n", wtime );

  return 0;
}