#include <iostream>
#include <gsl/gsl_rng.h>

#include "rtnorm.hpp"
#include <gsl/gsl_cdf.h>


//--génération quantiles--
std::vector<double> quantile_norm(int n, double sigma){
  std::vector<double> q(n);
  for (int i=0; i<n; i++) {
    q[i] = gsl_cdf_gaussian_Pinv ((double)(i+1)/n, sigma);
  }
  return q;
}

//--tirage de normale tronquée entre les quantile de taille 1/n i et i+1--
double quantile_truncate_normal (int i, int n, double mu, 
 double sigma, gsl_rng *gen) {
	 std::vector<double> q;
	 q = quantile_norm(n, sigma);
	 std::pair<double, double> p;
	 p = rtnorm (gen, q[i], q[i+1], mu, sigma);
	 return p.first;
 }
  	 


int main()
{
  //--- GSL random init ---
  gsl_rng_env_setup();                          // Read variable environnement
  const gsl_rng_type* type = gsl_rng_default;   // Default algorithm 'twister'
  gsl_rng *gen = gsl_rng_alloc (type);          // Rand generator allocation
  std::vector<double> q;
  q = quantile_norm(10, 1);
  std::cout<<q[5]<<std::endl;
  double a;
  int i;
  for (i=0; i<20; i++){
  a= quantile_truncate_normal(4, 10, 0, 1, gen);
  std::cout<<a<<std::endl;}
  return 0;
}