NumCpp  2.10.1
A Templatized Header Only C++ Implementation of the Python NumPy Library
trapazoidal.hpp
Go to the documentation of this file.
1 #pragma once
33 
34 #include <functional>
35 
36 #include "NumCpp/Core/Types.hpp"
37 
38 namespace nc::integrate
39 {
40  //============================================================================
41  // Method Description:
51  inline double trapazoidal(const double low,
52  const double high,
53  const uint32 n,
54  const std::function<double(double)>& f) noexcept
55  {
56  const double width = (high - low) / static_cast<double>(n);
57 
58  double trapezoidal_integral = 0.;
59  for (uint32 step = 0; step < n; ++step)
60  {
61  const double x1 = low + static_cast<double>(step) * width;
62  const double x2 = low + static_cast<double>(step + 1) * width;
63 
64  trapezoidal_integral += 0.5 * (x2 - x1) * (f(x1) + f(x2));
65  }
66 
67  return trapezoidal_integral;
68  }
69 } // namespace nc::integrate
Definition: gauss_legendre.hpp:43
double trapazoidal(const double low, const double high, const uint32 n, const std::function< double(double)> &f) noexcept
Definition: trapazoidal.hpp:51
dtype f(GeneratorType &generator, dtype inDofN, dtype inDofD)
Definition: f.hpp:56
std::uint32_t uint32
Definition: Types.hpp:40