From 355e4567e68a76356714e2e58a42dcd78533cf6c Mon Sep 17 00:00:00 2001
From: Bertrand <bertrand.horel@gmail.com>
Date: Fri, 19 Feb 2016 14:36:20 +0000
Subject: un peu de nettoyage

---
 rqmc.hpp | 14 +++++++++-----
 1 file changed, 9 insertions(+), 5 deletions(-)

(limited to 'rqmc.hpp')

diff --git a/rqmc.hpp b/rqmc.hpp
index 415ab00..1b33713 100644
--- a/rqmc.hpp
+++ b/rqmc.hpp
@@ -79,22 +79,26 @@ std::vector<double> monte_carlo(int n, Fct f, L X)
 
 //Les classes de quasi-mean
 
-
+template <typename LDS>
 struct quasi_gaussian : public var_alea<std::vector<double> >
 {
     quasi_gaussian (int d, double mean = 0, double std =1)
-        : d(d), mean(mean), std(std), sob(d) {};
+        : d(d), mean(mean), std(std), s(d), U(0,1), seed(d) {
+            for(int i=0; i<d; i++) {seed[i]=U();}
+                };
     std::vector<double> operator ()(){
-        std::vector<double> result = sob();
+        std::vector<double> result = s();
         for(int i=0;i<d; i++){
-            result[i] = mean + gsl_cdf_gaussian_Pinv(result[i], std);
+            result[i] = mean + gsl_cdf_gaussian_Pinv(frac_part(result[i]+seed[i]), std);
         }
         return value = result;
     }
     private: 
     int d;
     double mean, std;
-    sobol sob;
+    LDS s;
+    uniform U;
+    std::vector<double> seed;
 };
 
 
-- 
cgit v1.2.3-70-g09d2