HaltonSequence.h

Go to the documentation of this file.

 
 
#pragma once
 
#include <cmath>
#include <utility>
 
namespace INDI
{
namespace AlignmentSubsystem
{
 
class HaltonSequence
{
    public:
        HaltonSequence(int baseX = 2, int baseY = 3) : m_baseX(baseX), m_baseY(baseY) {}
 
        static double term(int i, int base)
        {
            double result = 0.0;
            double f      = 1.0;
            while (i > 0)
            {
                f /= static_cast<double>(base);
                result += (i % base) * f;
                i /= base;
            }
            return result;
        }
 
        std::pair<double, double> getRaw(int n) const
        {
            return {term(n, m_baseX), term(n, m_baseY)};
        }
 
        double getX(int n) const
        {
            return term(n, m_baseX);
        }
 
        double getY(int n) const
        {
            return term(n, m_baseY);
        }
 
    private:
        int m_baseX;
        int m_baseY;
};
 
} // namespace AlignmentSubsystem
} // namespace INDI