# include <cstdlib>
# include <iostream>
# include <cmath>
# include <ctime>
//
//  For your OpenMP version, you need to include the OpenMP file <omp.h>
//
using namespace std;

int main ( int argc, char *argv[] );
double cpu_time ( void );

//****************************************************************************80

int main ( int argc, char *argv[] )

//****************************************************************************80
//
//  Purpose:
//
//    MAIN is the main program for QUAD.
//
//  Discussion:
//
//    QUAD is a sketch of a C++ program to estimate the integral of
//    a function F(X) between A and B.
//
//    This program is only a sketch.  Your job is to finish it,
//    debug it, and get it to run.
//
//    Experiment to get a value of N that gives you a 'reasonable' accuracy.
//
//    NOW modify the program so that it will run under OpenMP.  Replace
//    CPU time by wall clock time.
//
//    Show that your program STILL gets the right answer.
//
//    Then determine the speedup when you run the OpenMP version on
//    1, 2, 4, 8 processors.
//
{
  double a = 0.0;
  double b = 10.0;

  double error;
  double exact = 0.49936338107645674464;
  int flops;
  int i;
  double mflops;
  int n = 1000000;
  double pi = 3.141592653589793;
  double total;
  double wtime;
  double wtime1;
  double wtime2;
  double *x;

  cout << "\n";
  cout << "QUAD:\n";
  cout << "  C++ version\n";
  cout << "\n";
  cout << "  Estimate the integral of f(x) from A to B.\n";
  cout << "  f(x) = 50 / ( pi * ( 2500 * x * x + 1 ) ).\n";
  cout << "  A = " << a << "\n";
  cout << "  B = " << b << "\n";
  cout << "  Exact integral from 0 to 10 is 0.49936338107645674464...\n";

  x = new double[n];
//
//  Step 1:
//  Load the array X with evenly spaced values between A and B.
//
  for ( i = 0; i < n; i++ )
  {
    x[i] =  ( ( double ) ( n - i - 1 ) * a   
            + ( double ) (     i     ) * b ) 
            / ( double ) ( n     - 1 );
  }
//
//  Step 2:
//
//  Use a FOR loop to add to TOTAL the value of the function at each X.
// 
//  Get your estimate by multiplying TOTAL by ( B - A ) and dividing by N.
//
//  Set FLOPS, the number of floating point operations by replacing "?"
//  by the number of floating point operations in your calculation of TOTAL.
//
//  add the OpenMP include file,
//
//  insert an OpenMP directive before this loop of the form
//
//      # pragma omp parallel for private ( ... ) shared ( ... ) reduction ( + : ... )
//
  wtime1 = omp_get_wtime ( );
  total = 0.0;

  MISSING OPENMP DIRECTIVE
  MISSING FOR STATEMENT
  MISSING CALCULATION OF TOTAL

  flops = ? * n;
  wtime2 = omp_get_wtime ( );
  total = ( b - a ) * total / ( double ) n;
//
//  Step 3:
//  Print quadrature estimate, Error, CPU time, MFLOPS rate
//
  error = fabs ( total - exact );
  wtime = wtime2 - wtime1;
  mflops = ( double ) ( flops ) / 1000000.0 / wtime;

  cout << "\n";
  cout << "  Estimate = " << total << "\n";
  cout << "  Error    = " << error << "\n";
  cout << "  W time   = " << wtime << "\n";
  cout << "  FLOPS    = " << flops << "\n";
  cout << "  MFLOPS   = " << mflops << "\n";

  delete [] n;

  return 0;
}