integrator.hpp

Go to the documentation of this file.

// -------------------------------------------------------//
//
// SHAMROCK code for hydrodynamics
// Copyright (c) 2021-2026 Timothée David--Cléris <tim.shamrock@proton.me>
// SPDX-License-Identifier: CeCILL Free Software License Agreement v2.1
// Shamrock is licensed under the CeCILL 2.1 License, see LICENSE for more information
//
// -------------------------------------------------------//
 
#pragma once
 
#include <functional>
 
namespace shammath {
 
    template<class T, class Lambda>
    inline constexpr T integ_riemann_sum(T start, T end, T step, Lambda &&fct) {
        T acc = {};
 
        for (T x = start; x < end; x += step) {
            acc += fct(x) * step;
        }
        return acc;
    }
 
} // namespace shammath