integrator.hpp

#pragma once
 
#include <nda/nda.hpp>
 
#include <utility>
 
namespace triqs::experimental::utility {


  template <typename D, int I> auto integrate(auto const &integrator1d, auto expr_to_integrate, nda::clef::pair<I, D> const &p) {
    auto f = nda::clef::make_function(std::move(expr_to_integrate), nda::clef::placeholder<I>{});
    return integrator1d(std::move(f), p.rhs);
  }

 
} // namespace triqs::experimental::utility