From a86b1d918a4c31f74024aa6b42fa0021919585de Mon Sep 17 00:00:00 2001
From: Bertrand <bertrand.horel@gmail.com>
Date: Mon, 8 Feb 2016 14:40:20 +0000
Subject: débuggage et seeding
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 stratified_sampling.cpp | 38 --------------------------------------
 1 file changed, 38 deletions(-)

(limited to 'stratified_sampling.cpp')

diff --git a/stratified_sampling.cpp b/stratified_sampling.cpp
index b959776..c1c8dc1 100644
--- a/stratified_sampling.cpp
+++ b/stratified_sampling.cpp
@@ -1,6 +1,4 @@
 #include "stratified_sampling.hpp"
-#include <algorithm>
-#include <iostream>
 
 std::pair<double, double> mean_var( std::vector<double> r){
     std::pair<double, double> p;
@@ -13,39 +11,3 @@ std::pair<double, double> mean_var( std::vector<double> r){
     p.second -= p.first * p.first;
     return p;
 }
-
-//actualisation du nombre de tirages à faire par strates
-template <typename Gen>
-void stratified_sampling<Gen>::update(int Nk) {
-    int I = p.size();
-    //reinitialistation du vecteur M du nombre de tirages par strates
-    if (M.empty()) {
-        M.resize(I,1);
-    }
-    else {
-        for(int i=0; i<I; i++){
-            M[i]=1;
-        }
-    }
-
-    std::vector<double> m(I, 0); //le vecteur des m_i idéals
-
-    if (sigma.empty()) {
-        for (int i=0; i<I; i++) {
-            m[i] = (Nk-I)*p[i];
-        }
-    }
-    else {
-        double scal = std::inner_product(p.begin(), p.end(), sigma.begin(), (double) 0);
-        for (int i=0; i < I; i++) {
-            m[i] = (Nk-I)*p[i]*sigma[i]/scal;
-            std::cout<<m[i]<<std::endl;
-        }
-    }
-    M[0]+=floor(m[0]);
-    double current = m[0];
-    for (int i=1; i<I; i++){
-        M[i] += floor(current+m[i]) - floor(current);
-        current += m[i];
-    }
-}
-- 
cgit v1.2.3-70-g09d2