element.h

#ifndef ELEMENT_H
#define ELEMENT_H
 
#include <assert.h>
#include <cmath>
#include <limits>
 
#include "types.h"
 
namespace svec {
 
struct Element {
    Label l;
 
    Value v;
 
    constexpr operator Value() const
    {
        return v;
    }
 
    Element& operator+=(const Value& rhs);
    Element& operator+=(const Element& rhs);
 
    Element& operator-=(const Value& rhs);
    Element& operator-=(const Element& rhs);
 
    Element& operator*=(const Value& rhs);
    Element& operator/=(const Value& rhs);
 
    bool isEnd() const;
 
    friend constexpr Element operator+(const Element& a, const Element& b);
    friend constexpr Element operator-(const Element& a, const Element& b);
    friend constexpr Element operator*(const Element& a, const Value& C);
    friend constexpr Element fma(const Element& a, const Value& C, const Element& b);
};
 
constexpr Element operator+(const Element& a, const Element& b)
{
    assert(a.l == b.l);
    return {a.l, a.v + b.v};
}
 
constexpr Element operator-(const Element& a, const Element& b)
{
    assert(a.l == b.l);
    return {a.l, a.v - b.v};
}
 
constexpr Element operator*(const Element& a, const Value& C)
{
    return {a.l, a.v * C};
}
 
constexpr Element fma(const Element& a, const Value& C, const Element& b)
{
    assert(a.l == b.l);
    return {a.l, std::fma(a.v, C, b.v)};
}
 
constexpr Element END_ELEMENT = {std::numeric_limits<Label>::max(), 0.0};
 
inline Element& Element::operator+=(const Value& rhs)
{
    this->v += rhs;
    return *this;
};
inline Element& Element::operator+=(const Element& rhs)
{
    assert(this->l == rhs.l);
    return operator+=(rhs.v);
};
 
inline Element& Element::operator-=(const Value& rhs)
{
    this->v -= rhs;
    return *this;
};
inline Element& Element::operator-=(const Element& rhs)
{
    assert(this->l == rhs.l);
    return operator-=(rhs.v);
};
 
inline Element& Element::operator*=(const Value& rhs)
{
    this->v *= rhs;
    return *this;
};
inline Element& Element::operator/=(const Value& rhs)
{
    this->v /= rhs;
    return *this;
}
 
inline bool Element::isEnd() const
{
    // checking the label is sufficient
    return this->l == END_ELEMENT.l;
}
 
} // namespace svec
 
#endif