SimpleRules.h

Go to the documentation of this file.

#pragma once

#include "AQ_Util.h"

namespace UA_CoMP {
namespace Num {

class SimpsonRule
{
public:
        template <class Ftor>
        static typename Ftor::result_type sum(typename Ftor::argument_type l, typename Ftor::argument_type d, Ftor f)
        {
                using Arg = typename StripConstRef<typename Ftor::argument_type>::type;
                using Res = typename StripConstRef<typename Ftor::result_type>::type;

                const Arg h{static_cast<Arg>(0.5) * d};
                const Arg r{l + d};
                return static_cast<Res>(0.333333333333333) * h * (f(l) + static_cast<Res>(4.0) * f(l + h) + f(r));
        }
};

} // namespace Num
} // namespace UA_CoMP