WeightedSum.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
namespace MathLib
{
 
namespace detail {
 
template <unsigned I, typename Method>
struct SUM
{
    template <typename F>
    static
    double
    add(F const& f)
    {
        return f(Method::X[I - 1]) * Method::W[I - 1] + SUM<I - 1, Method>::add(f);
    }
};

template <typename Method>
struct SUM<0, Method>
{
    template <typename F>
    static double add(F const& /*unused*/)
    {
        return 0;
    }
};
 
}   // namespace detail

template <typename Method>
struct WeightedSum
{
    template <typename Func>
    static
    double
    add(Func const& f)
    {
        return detail::SUM<Method::Order, Method>::add(f);
    }
};
 
}   // namespace MathLib