integrals.H

Go to the documentation of this file.

#ifndef INTEGRALS_H
#define INTEGRALS_H
#include <AMReX_REAL.H>
 
namespace amr_wind::utils {
 
template <typename Function>
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real trapz(
    const amrex::Real xa, const amrex::Real xb, const int n, const Function& f)
{
 
    const amrex::Real dx = (xb - xa) / n;
 
    // first and last terms
    amrex::Real sum = 0.5 * (f(xa) + f(xb));
 
    for (int i = 1; i < n; i++) {
        sum += f(xa + i * dx);
    }
    return dx * sum;
}
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real trapz( {…}
 
} // namespace amr_wind::utils
#endif