From bcafe5e41698c44b6aa77b8ef150f4224613c38f Mon Sep 17 00:00:00 2001
From: Bertrand <bertrand.horel@gmail.com>
Date: Tue, 12 Apr 2016 17:05:23 +0000
Subject: rectification, renommage quasi_mean en monte_carlo, renommage
 quasi_option en quasi_mean qui est déplacée dans rqmc.hpp
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 src/rqmc.hpp | 25 ++++++++++++++++++++++---
 1 file changed, 22 insertions(+), 3 deletions(-)

(limited to 'src/rqmc.hpp')

diff --git a/src/rqmc.hpp b/src/rqmc.hpp
index 03f7e81..840e069 100644
--- a/src/rqmc.hpp
+++ b/src/rqmc.hpp
@@ -56,7 +56,7 @@ mean(int n, VA X){
 //Les classes de monte-Carlo
 
 template <typename L>
-std::vector<double> quasi_mean(int n, L X, double p=0.05)
+std::vector<double> monte_carlo(int n, L X, double p=0.05)
 {
 	std::vector<double> result(3,0);
 	double x;
@@ -72,9 +72,9 @@ std::vector<double> quasi_mean(int n, L X, double p=0.05)
 };
 
 template <typename L, typename Fct>
-std::vector<double> quasi_mean(int n, Fct f, L X, double p=0.05)
+std::vector<double> monte_carlo(int n, Fct f, L X, double p=0.05)
 {
-	return quasi_mean(n, compose(f, X), p); 
+	return monte_carlo(n, compose(f, X), p); 
 };
 
 //Les classes de quasi-mean
@@ -120,6 +120,25 @@ struct gaussian_d : public var_alea<std::vector<double> >
     gaussian G;
 };
 
+template <typename Fct, typename LDS>
+struct quasi_mean : public generator<typename Fct::result_type>
+{
+    quasi_mean(int n, int d, Fct f) : n(n), d(d), f(f), G(d) {};
+        
+    typename Fct::result_type operator()() {
+        double sum =0;
+        for(int i=0; i<n; i++){
+            sum += f(G());
+            
+        }
+        return sum/n;
+    };
+    
+    private:
+        int n, d;
+        Fct f;
+        quasi_gaussian<LDS> G;
+    };
 
 
 
-- 
cgit v1.2.3-70-g09d2