stratified_sampling.cpp


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28

#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){
	vector<double> quantile_norm(n);
	for (int i=0, i<n, i++) {
	quantile_norm[i]=gsl_cdf_gaussian_Pinv (i/n, sigma);
	}
	return quantile_norm;
}

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

  
  return 0;
}