ExactPosition.hpp

#pragma once
 
#include <wmtk/Types.hpp>
 
#include <memory>
 
namespace wmtk {
 
template <class VecR, class VecD>
class ExactPosition
{
public:
    ExactPosition() = default;
    ExactPosition(const ExactPosition& o)
        : m_p(o.m_p ? std::make_unique<VecR>(*o.m_p) : nullptr)
    {}
    ExactPosition& operator=(const ExactPosition& o)
    {
        if (this == &o) return *this;
        if (!o.m_p) {
            m_p.reset();
        } else if (m_p) {
            *m_p = *o.m_p;
        } else {
            m_p = std::make_unique<VecR>(*o.m_p);
        }
        return *this;
    }
 
    VecR value(const VecD& posf) const { return m_p ? *m_p : to_rational(posf); }
 
    bool coord_equals(int k, const VecD& posf, const Rational& r) const
    {
        return m_p ? (*m_p)[k] == r : Rational(posf[k]) == r;
    }
 
    void store(const VecR& p)
    {
        if (m_p) {
            *m_p = p;
        } else {
            m_p = std::make_unique<VecR>(p);
        }
    }
 
    void clear() { m_p.reset(); }
 
    bool stored() const { return static_cast<bool>(m_p); }
 
private:
    std::unique_ptr<VecR> m_p;
};
 
} // namespace wmtk