#ifndef SCALE_H
#define SCALE_H

#include <algorithm>

template<typename T>
class Scale
{
public:
    Scale(const T &domainMin, const T &domainMax, const T &rangeMin, const T &rangeMax)
        : m_domainMin(domainMin)
        , m_domainMax(domainMax)
        , m_rangeMin(rangeMin)
        , m_rangeMax(rangeMax)
    {
        valuesUpdated();
    }

    virtual ~Scale() {}

    T rangeMin() const  { return m_rangeMin; }
    T rangeMax() const  { return m_rangeMax; }
    T domainMin() const { return m_domainMin; }
    T domainMax() const { return m_domainMax; }

    void setRangeMin(T rangeMin)   { m_rangeMin  = rangeMin; valuesUpdated(); }
    void setRangeMax(T rangeMax)   { m_rangeMax  = rangeMax; valuesUpdated(); }
    void setDomainMin(T domainMin) { m_domainMin = domainMin; valuesUpdated(); }
    void setDomainMax(T domainMax) { m_domainMax = domainMax; valuesUpdated(); }

    void setRange(const T &rangeMin, const T &rangeMax)
    {
        m_rangeMin = rangeMin;
        m_rangeMax = rangeMax;
        valuesUpdated();
    }

    void setRange(const Scale<T> &other)
    {
        setRange(other.m_rangeMin, other.m_rangeMax);
    }

    void setDomain(const T &domainMin, const T &domainMax)
    {
        m_domainMin = domainMin;
        m_domainMax = domainMax;
        valuesUpdated();
    }

    void setDomain(const Scale<T> &other)
    {
        setDomain(other.m_domainMin, other.m_domainMax);
    }

    void inverse()
    {
        std::swap(m_rangeMin, m_domainMin);
        std::swap(m_rangeMax, m_domainMax);
        valuesUpdated();
    }

    virtual T operator()(const T &value) const = 0;

protected:
    // Called when internal values change
    virtual void valuesUpdated() {}

    T m_domainMin, m_domainMax;
    T m_rangeMin, m_rangeMax;
};

template<typename T>
class LinearScale
    : public Scale<T>
{
public:
    LinearScale(const T &domainMin, const T &domainMax, const T &rangeMin, const T &rangeMax)
        : Scale<T>(domainMin, domainMax, rangeMin, rangeMax)
    {
        valuesUpdated();
    }

    virtual T operator()(const T &value) const
    {
        return (value - Scale<T>::m_domainMin) * m_transformSlope + Scale<T>::m_rangeMin;
    }

    T slope() const { return m_transformSlope; }

    T offset() const
    {
        return Scale<T>::m_rangeMin - Scale<T>::m_domainMin * m_transformSlope;
    }

protected:
    virtual void valuesUpdated()
    {
        m_transformSlope = (Scale<T>::m_rangeMax - Scale<T>::m_rangeMin)
                         / (Scale<T>::m_domainMax - Scale<T>::m_domainMin);
    }

private:
    T m_transformSlope;
};

void updateTransform4x4(const LinearScale<float> &sx,
                        const LinearScale<float> &sy,
                        float transform[4][4]);

void updateTransform4x4(LinearScale<float> &sx,
                        LinearScale<float> &sy,
                        float offsetX,
                        float offsetY,
                        float transform[4][4]);

#endif // SCALE_H