// Random numbers
#include <iomanip>
#include <fstream>
#include <algorithm>
#include <algorithm>
#include "boost_random_numbers.hpp"

using namespace std;

// utility programs:
void format_output(std::ostream& out, unsigned decimals){
  out << std::fixed<<std::setprecision(decimals);
} 

void vector_out(vector<double> v)
{
  copy(v.begin(),v.end(), // from
       ostream_iterator<double>(cout," ")); // to 
  cout<<endl;
}

int main(){  
  ofstream outunif("unif"),outgauss("gauss");
  format_output(outunif,30);
  format_output(outgauss,30);
  cout<<"U[0,1) to file 'unif'\n";
  cout<<"Normal (Gaussian) distribution to file 'gauss'\n";
  
  for(int i = 0; i < 100000; i++) {
    outunif  << unirand() << "\n";
    outgauss << gaussrand() << "\n";
  }
  outunif.close();
  outgauss.close();
  cout<<"done\n";
  cout<<endl;
  vector<double> vran(100); 
  generate(vran.begin(),vran.end(),unirand);
  cout<<"vran = \n";
  vector_out(vran);
  return 0;
}