#include <upc.h>
#include <stdio.h>
#include <stdlib.h>
#include <math.h>

shared int all_hits[THREADS];

int main (int argc, char **argv) {
  int i, hits = 0, trials = 0, my_trials;
  double pi;
  
  if (argc != 2) trials = 10000000;
  else trials = atoi (argv[1]);
  
  my_trials = trials / THREADS;
  
  srand (MYTHREAD * 13);
  
  for (i=0; i<my_trials; ++i)
    all_hits[MYTHREAD] += hit ();

  upc_barrier;

  if (MYTHREAD == 0) {
    for (i=0; i<THREADS; ++i) hits += all_hits[i];
    printf ("PI estimated to %f.\n", 4.0 * hits / trials);    
  }
}
  
int hit () {
	double x = rand() / (double) RAND_MAX;
	double y = rand() / (double) RAND_MAX;
	if ((x*x + y*y) <= 1.0) return (1);
	else return(0);
}