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

---
 src/rqmc.cpp | 38 ++++++++++++++++++++++++++++++++++++++
 1 file changed, 38 insertions(+)
 create mode 100644 src/rqmc.cpp

(limited to 'src/rqmc.cpp')

diff --git a/src/rqmc.cpp b/src/rqmc.cpp
new file mode 100644
index 0000000..7109d82
--- /dev/null
+++ b/src/rqmc.cpp
@@ -0,0 +1,38 @@
+#include "rqmc.hpp"
+
+
+double frac_part(double x){
+    return x - floor(x);
+}
+
+
+double mean_rqmc(int N, double X) {
+    sobol s(1);
+    double sum = 0;
+    for(int i=0; i<N; i++){
+        sum +=gsl_cdf_gaussian_Pinv (frac_part(X+s()[0]), 1);
+    }
+    return sum/N;
+}
+
+//~ int main() {
+    //~ init_alea(0);
+    //~ int I=100;
+    //~ int N= 10000;
+    //~ uniform U;
+    //~ double m = 0;
+    //~ double s = 0;
+    //~ double temp;
+    
+    //~ for(int i=0;i<I;i++){
+        //~ temp = mean_rqmc(N,U());
+        //~ m+=temp;
+        //~ s+=temp*temp;
+    //~ }
+    //~ m = m/I;
+    //~ s = s/I - m*m;
+    
+    //~ std::cout<<"espérance "<<m<<" taille de l'IC "<<sqrt(s)*1.96/10<<std::endl;
+
+    //~ return 0;
+//~ }
-- 
cgit v1.2.3-70-g09d2