From d2b133901a65244934eb642ec8e20c797efaf650 Mon Sep 17 00:00:00 2001
From: Bertrand <bertrand.horel@gmail.com>
Date: Fri, 19 Feb 2016 15:03:51 +0000
Subject: nettoyage du dépôt
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 main.cpp | 53 -----------------------------------------------------
 1 file changed, 53 deletions(-)
 delete mode 100644 main.cpp

(limited to 'main.cpp')

diff --git a/main.cpp b/main.cpp
deleted file mode 100644
index f81f1dd..0000000
--- a/main.cpp
+++ /dev/null
@@ -1,53 +0,0 @@
-#include <iostream>
-#include <gsl/gsl_rng.h>
-#include <vector>
-#include <gsl/gsl_cdf.h>
-#include <gsl/gsl_math.h>
-#include "stratified_sampling.hpp"
-#include <cmath>
-#include <algorithm>
-
-using namespace std;
-//--génération quantiles--
-vector<double> quantile_norm(int n, double sigma){
-    vector<double> q(n);
-    for (int i=0; i<n; i++) {
-        q[i] = gsl_cdf_gaussian_Pinv ((double)(i+1)/n, sigma);
-    }
-    return q;
-}
-
-int main()
-{
-    //--- GSL random init ---
-    gsl_rng_env_setup();
-    vector<double> q = quantile_norm(10, 1);
-    vector<double> p(10, 0.1);
-    vector<gaussian_truncated> rvar;
-    rvar.push_back(gaussian_truncated(GSL_NEGINF, q[0],0,1,0));
-    for (int i=1; i<10; i++){
-        rvar.push_back(gaussian_truncated(q[i-1], q[i],0,1,i));
-    }
-    stratified_sampling<gaussian_truncated> S(p,rvar);
-    //S.update(100);
-    S.draw(100);
-    double x = 1.64*S.estimator().second;
-		cout<<"l'estimateur de la moyenne est :"<<S.estimator().first<<endl;
-		cout<<"Son intervalle de confiance à 95% est :"<<"["<<S.estimator().first-(x/10)<<" ,"<<S.estimator().first+(x/10)<<"]"<<endl;
-    //S.update(1000);
-    S.draw(1000);
-		x = 1.64*S.estimator().second;
-		cout<<"l'estimateur de la moyenne est :"<<S.estimator().first<<endl;
-		cout<<"Son intervalle de confiance à 95% est :"<<"["<<S.estimator().first-(x/sqrt(1100))<<" ,"<<S.estimator().first+(x/sqrt(1100))<<"]"<<endl;
-		//S.update(10000);
-		S.draw(10000);
-		x = 1.64*S.estimator().second;
-		cout<<"l'estimateur de la moyenne est :"<<S.estimator().first<<endl;
-		cout<<"Son intervalle de confiance à 95% est :"<<"["<<S.estimator().first-(x/sqrt(11100))<<" ,"<<S.estimator().first+(x/sqrt(11100))<<"]"<<endl;
-
-    //~ S.draw();
-    //~ for(int i=0;i<10;i++){
-        //~ cout<<S.get_mean()[i]<<endl;
-    //~ }
-    return 0;
-}
-- 
cgit v1.2.3-70-g09d2