// std::mt19937 random number generator
// uses std::function and std::bind
// g++ -std=c++0x std_mersenne.cpp utility.cpp
#include <iostream>
#include <random>
#include <ctime> 
#include <fstream>
#include <functional>
#include "utility.hpp"

using namespace std;

int main()
{
  // define generator
  mt19937 gener;
  // initialize generator
  auto seed = static_cast<uint_fast32_t> (std::time(0));
  cout<<" seed = "<<seed<<endl;
  gener.seed(seed);  

  // define distribution
  uniform_real_distribution<double> unif;
  // bind distribution and generator
  function<double()> unirand = bind(unif, gener);
  // test it
  ofstream out1("unif");
  format_output(out1,30);
  for (int i = 0; i<10000; ++i) out1<< unirand()<<endl;
  out1.close();

  normal_distribution<double> gaussian;
  function<double()> gaussrand = bind(gaussian, gener);
  ofstream out2("gauss");
  format_output(out2,30);
  for (int i = 0; i<10000; ++i) out2<< gaussrand()<<endl;
  out2.close();

  exponential_distribution<double> expon;
  function<double()> exprand = bind(expon, gener);
  ofstream out3("expon");
  format_output(out3,30);
  for (int i = 0; i<10000; ++i) out3<< exprand()<<endl;
  out3.close();


  return 0;
}