Integral.h

Go to the documentation of this file.

#pragma once

#include <functional>

namespace UA_CoMP {
namespace Num {

template <typename QuadAlgorithm, typename Integrand>
class Integral : public std::binary_function<typename Integrand::argument_type, typename Integrand::argument_type,
                                             typename Integrand::result_type>
{
public:
        Integral(QuadAlgorithm const& quad, Integrand const& ftor) : fQuadAlgorithm(quad), fIntegrand(ftor) {}

        typename Integral::result_type operator()(typename Integral::first_argument_type l,
                                                  typename Integral::first_argument_type r) const
        {
                return fQuadAlgorithm(l, r, fIntegrand);
        }

private:
        QuadAlgorithm fQuadAlgorithm;
        Integrand fIntegrand;
};

template <typename QuadAlgorithm, typename Integrand>
inline Integral<QuadAlgorithm, Integrand> make_integral(QuadAlgorithm const& quad, Integrand const& ftor)
{
        return Integral<QuadAlgorithm, Integrand>(quad, ftor);
}

} // namespace Num
} // namespace UA_CoMP