From 71c709b0ba6d6fe9ba043348a074b612f2e1f639 Mon Sep 17 00:00:00 2001
From: Laura Messio <laura.messio@cea.fr>
Date: Mon, 22 Jul 2013 17:52:05 +0200
Subject: [PATCH] gf with products of meshes : bug corrected and commented

---
 triqs/gf/re_im_time.hpp     |  2 +-
 triqs/gf/refreq_imtime.hpp  |  2 +-
 triqs/gf/two_real_times.hpp | 12 +++++-------
 3 files changed, 7 insertions(+), 9 deletions(-)

diff --git a/triqs/gf/re_im_time.hpp b/triqs/gf/re_im_time.hpp
index 270c93a2..05c3c7b6 100644
--- a/triqs/gf/re_im_time.hpp
+++ b/triqs/gf/re_im_time.hpp
@@ -83,7 +83,7 @@ namespace triqs { namespace gf {
    static gf_t make_gf(MeshType && m) {
     typename gf_t::data_non_view_t A(m.size());
     A() =0;
-    return gf_t (m, std::move(A), nothing(), nothing(), nothing() ) ;
+    return gf_t (m, std::move(A), nothing(), nothing() ) ;
    }
 
    static gf_t make_gf(double tmin, double tmax, size_t nt, double beta, statistic_enum S, size_t ntau, mesh_kind mk=full_bins) { 
diff --git a/triqs/gf/refreq_imtime.hpp b/triqs/gf/refreq_imtime.hpp
index 9f292986..cf4c1dbf 100644
--- a/triqs/gf/refreq_imtime.hpp
+++ b/triqs/gf/refreq_imtime.hpp
@@ -82,7 +82,7 @@ namespace triqs { namespace gf {
      static gf_t make_gf(MeshType && m) {
       typename gf_t::data_non_view_t A(m.size());
       A() =0;
-      return gf_t (m, std::move(A), nothing(), nothing(), nothing() ) ;
+      return gf_t (m, std::move(A), nothing(), nothing() ) ;
      }
      
      static gf_t make_gf(double wmin, double wmax, size_t n_freq, double beta, statistic_enum S, size_t nt, mesh_kind mk=full_bins) { 
diff --git a/triqs/gf/two_real_times.hpp b/triqs/gf/two_real_times.hpp
index 349947e1..ce59563e 100644
--- a/triqs/gf/two_real_times.hpp
+++ b/triqs/gf/two_real_times.hpp
@@ -56,13 +56,11 @@ namespace triqs { namespace gf {
    struct get_closest_point <two_real_times,matrix_valued,Opt> {
    typedef typename mesh<two_real_times, Opt>::type mesh_t;
     
-   // NOT FINISHED, NOT TESTED 
-    template<typename G, typename T>
-     static typename mesh_t::index_t invoke(G const * g, closest_pt_wrap<T,T> const & p) {
-      double x = (g->mesh().kind()==half_bins ? double(p.value) :  double(p.value)+ 0.5*g->mesh().delta());
-      size_t n = std::floor(x/g->mesh().delta());
-      return n;
-     }
+//    // NOT FINISHED, NOT TESTED 
+//     template<typename G, typename T>
+//      static typename mesh_t::index_t invoke(G const * g, closest_pt_wrap<T,T> const & p) {
+//       return std::floor( double(p.value) / g->mesh().delta() + 0.5);
+//      }
 
    };