Trapezoidal rule. More...

#include <benchmarkUtil.h>

Collaboration diagram for olb::util::TrapezRuleInt1D< T >:

Public Member Functions
	TrapezRuleInt1D (AnalyticalF1D< T, T > &f)

T	integrate (T min, T max, int nSteps)

Protected Attributes
AnalyticalF1D< T, T > &	_f

Detailed Description

template<typename T>
class olb::util::TrapezRuleInt1D< T >

Trapezoidal rule.

Definition at line 141 of file benchmarkUtil.h.

Constructor & Destructor Documentation

◆ TrapezRuleInt1D()

template<typename T >

olb::util::TrapezRuleInt1D< T >::TrapezRuleInt1D ( AnalyticalF1D< T, T > & f )

Definition at line 297 of file benchmarkUtil.hh.

                                                         : _f(f)
{
}

Member Function Documentation

◆ integrate()

template<typename T >

T olb::util::TrapezRuleInt1D< T >::integrate	(	T	min,
		T	max,
		int	nSteps )

Definition at line 302 of file benchmarkUtil.hh.

{
  T integral[1], tmp[1], min_[1], max_[1], x[1], deltax[0];
  min_[0] = min;
  max_[0] = max;
  deltax[0] = (max - min) / nSteps;
  _f(tmp, min);
  integral[0] = 0.5 * tmp[0];
  x[0] = min;
  for (int i = 1; i < nSteps; i++) {
    x[0] += deltax[0];
    _f(tmp, x);
    integral[0] += tmp[0];
  }
 
  _f(tmp, max);
  integral[0] += 0.5*tmp[0];
  integral[0] *= deltax[0];
 
  std::cout << "Itnegral=" <<  integral[0] << std::endl;
 
  return integral[0];
}