#include <stdio.h>
#include <stdlib.h>
#include "pi.h"


int main(int argc, char *argv[])
{
  long N;
  int ProcessRank, P;
  double pi;
  
  if(argc != 2)
  {
    fprintf(stderr,"USAGE: %s NumberOfProcesses\n", argv[0]);
    exit(-1);
  }
  
  P = atoi(*++argv);
  
  for(N=1000000; N <= 16000000; N *= 2)
  {
    pi = 0.0;
    for(ProcessRank=0; ProcessRank<P; ProcessRank++)
    {
      pi += EstimatePi(N,ProcessRank);
    }
    
    pi /= (double) P;
    printf( "Estimate of pi from %ld samples, using %d processes, = %f\n", N, P, pi);
  }
  
 

  return 0;
}