MHO_UniformGridPointsCalculator.hh

Go to the documentation of this file.

#ifndef MHO_UniformGridPointsCalculator_HH__
#define MHO_UniformGridPointsCalculator_HH__
 
#include <cmath>
#include <limits>
#include <map>
#include <vector>
 
namespace hops
{
 
class MHO_UniformGridPointsCalculator
{
    public:
        MHO_UniformGridPointsCalculator();
        ~MHO_UniformGridPointsCalculator();
 
        void SetEpsilon(double eps) { fEpsilon = std::fabs(eps); };
 
        void SetDefaultGridPoints(std::size_t n) { fDefaultGridPoints = n; };
 
 
 
        void GetUniquePoints(const std::vector<double>& in_pts, double eps, 
                        std::vector<double>& out_pts, 
                        std::map< std::size_t, std::size_t >& index_map) const;
 
 
        void SetPoints(const std::vector< double >& pts);
        
        void SetPoints(const double* pts, std::size_t npts);
        
        void Calculate();
 
        void GetGridPoints(std::vector< double >* grid_pts); //fill vector with grid points
 
        double GetGridStart() const { return fStart; };
 
        double GetGridSpacing() const { return fSpacing; };
 
        double GetGridAverage() const { return fAverageLocation; }
 
        double GetSpread() const { return fSpread; }
 
        std::size_t GetNGridPoints() const { return fNGridPoints; };
 
        bool GetSpacingErrorStatus() const { return fSpacingError; };
 
 
        std::map< std::size_t, std::size_t > GetGridIndexMap() { return fIndexMap; }
 
    protected:
        void Calculate_v1(int max_pts = 8192); 
        
        void Calculate_v2();
 
        void FindStartAndMinMaxSpacing();
 
        double fEpsilon;
        double fStart;
        double fSpacing;
        double fMinSpacing;
        double fMaxSpacing;
        double fAverageLocation;
        double fNGridPoints;
        double fSpread;
        std::vector< double > fPoints;
        std::map< std::size_t, std::size_t > fIndexMap;
 
        int fDefaultGridPoints;
        bool fSpacingError;
        double fAbsEps; //use to check that value is not zero
};
 
} // namespace hops
 
#endif