// -- numeric.hpp -- Boost Lambda Library -----------------------------------

// Copyright (C) 2002 Jaakko Järvi (jaakko.jarvi@cs.utu.fi)

// Copyright (C) 2002 Gary Powell (gwpowell@hotmail.com)

//

// Distributed under the Boost Software License, Version 1.0. (See

// accompanying file LICENSE_1_0.txt or copy at

// http://www.boost.org/LICENSE_1_0.txt)

//

// For more information, see http://www.boost.org

#ifndef BOOST_LAMBDA_NUMERIC_HPP

#define BOOST_LAMBDA_NUMERIC_HPP

#include "boost/lambda/core.hpp"

#include <numeric>

namespace boost {

  namespace lambda {

namespace ll {

// accumulate ---------------------------------

struct accumulate {

  template <class Args>

  struct sig { 

    typedef typename boost::remove_const<

        typename boost::tuples::element<3, Args>::type 

     >::type type; 

  };

  template <class A, class B, class C>

  C

  operator()(A a, B b, C c) const

  { return ::std::accumulate(a, b, c); }

  template <class A, class B, class C, class D>

  C

  operator()(A a, B b, C c, D d) const

  { return ::std::accumulate(a, b, c, d); }

};

// inner_product ---------------------------------

struct inner_product {

  template <class Args>

  struct sig { 

    typedef typename boost::remove_const<

        typename boost::tuples::element<4, Args>::type 

     >::type type; 

  };

  template <class A, class B, class C, class D>

  D

  operator()(A a, B b, C c, D d) const

  { return ::std::inner_product(a, b, c, d); }

  template <class A, class B, class C, class D, class E, class F>

  D

  operator()(A a, B b, C c, D d, E e, F f) const

  { return ::std::inner_product(a, b, c, d, e, f); }

};

// partial_sum ---------------------------------

struct partial_sum {

  template <class Args>

  struct sig { 

    typedef typename boost::remove_const<

        typename boost::tuples::element<3, Args>::type 

     >::type type; 

  };

  template <class A, class B, class C>

  C

  operator()(A a, B b, C c) const

  { return ::std::partial_sum(a, b, c); }

  template <class A, class B, class C, class D>

  C

  operator()(A a, B b, C c, D d) const

  { return ::std::partial_sum(a, b, c, d); }

};

// adjacent_difference ---------------------------------

struct adjacent_difference {

  template <class Args>

  struct sig { 

    typedef typename boost::remove_const<

        typename boost::tuples::element<3, Args>::type 

     >::type type; 

  };

  template <class A, class B, class C>

  C

  operator()(A a, B b, C c) const

  { return ::std::adjacent_difference(a, b, c); }

  template <class A, class B, class C, class D>

  C

  operator()(A a, B b, C c, D d) const

  { return ::std::adjacent_difference(a, b, c, d); }

};

} // end of ll namespace

} // end of lambda namespace

} // end of boost namespace

#endif