Integrator.hh

//
// This file is part of MARLEY (Model of Argon Reaction Low Energy Yields)
//
// MARLEY is free software: you can redistribute it and/or modify it under the
// terms of version 3 of the GNU General Public License as published by the
// Free Software Foundation.
//
// For the full text of the license please see COPYING or
// visit http://opensource.org/licenses/GPL-3.0
//
// Please respect the MCnet academic usage guidelines. See GUIDELINES
// or visit https://www.montecarlonet.org/GUIDELINES for details.
 
#pragma once
#include <functional>
#include <vector>
 
#include "marley/marley_utils.hh"
 
namespace marley {
 
  class Integrator {
    public:
 
      Integrator(size_t num = N_DEFAULT_);
 
      double num_integrate(const std::function<double(double)>& f,
        double a, double b) const;
 
    private:
 
      size_t N_;
 
      std::vector<double> weights_;
 
      std::vector<double> offsets_;
 
      static constexpr size_t N_DEFAULT_ = 20; 
  };
 
}