composition.hpp

#pragma once
 
#include <ostream>
#include <limits>
#include <armadillo>
 
class Composition
{
public:
    Composition();
 
    explicit Composition(const arma::vec& composition);
 
    inline double operator()(const arma::uword i) const { return m_composition(i); }
 
    inline double& operator()(const arma::uword i) { return m_composition(i); }
 
    // user-defined conversion (implicit conversion)
    // this allows Composition to be a drop-in replacement for the arma::vec
    // in most cases, which is both positive and negative
    inline operator arma::vec() const { return m_composition; }
 
    inline const arma::vec& vec() const { return m_composition; }
 
    inline arma::vec vec() { return m_composition; }
 
    friend std::ostream& operator << (std::ostream &out, const Composition& c);
 
    friend bool operator==(const Composition& lhs, const Composition& rhs);
 
    bool isNormalized() const;
 
    Composition& normalize();
 
    Composition normalized() const;
 
    const static arma::uword n_elem = 10;
 
    static const Composition defaultComposition; // initialized in .cpp
 
    // getters
    inline double C1()  const { return m_composition[0]; } 
    inline double C2()  const { return m_composition[1]; } 
    inline double C3()  const { return m_composition[2]; } 
    inline double iC4() const { return m_composition[3]; } 
    inline double nC4() const { return m_composition[4]; } 
    inline double iC5() const { return m_composition[5]; } 
    inline double nC5() const { return m_composition[6]; } 
    inline double C6()  const { return m_composition[7]; } 
    inline double N2()  const { return m_composition[8]; } 
    inline double CO2() const { return m_composition[9]; } 
 
    // setters
    inline double& C1()  { return m_composition[0]; } 
    inline double& C2()  { return m_composition[1]; } 
    inline double& C3()  { return m_composition[2]; } 
    inline double& iC4() { return m_composition[3]; } 
    inline double& nC4() { return m_composition[4]; } 
    inline double& iC5() { return m_composition[5]; } 
    inline double& nC5() { return m_composition[6]; } 
    inline double& C6()  { return m_composition[7]; } 
    inline double& N2()  { return m_composition[8]; } 
    inline double& CO2() { return m_composition[9]; } 
 
private:
    arma::vec::fixed<10> m_composition;
};