/* 
 *
 */
#include <stdio.h>

#include <gsl/gsl_rng.h>

#include "mc_pi.h"

#define POINTS 500000000L /* number of random points to generate */

int main(int argc, char **argv) {
    double pi;         /* calculated value of pi */
    gsl_rng *r;

    /* allocate random number generator */
    r = gsl_rng_alloc(gsl_rng_taus2);

    /* Set seed for random number generator */
    gsl_rng_set(r, (unsigned long int) (1));

    /* calculate pi using Monte Carlo algorithm */
    pi = mc_pi(POINTS, r);

    printf("   After %3ld throws, the value of pi = %10.8f\n", POINTS, pi);

    return(0);
}