discretizer.hpp

#pragma once
 
#include <armadillo>
 
class Pipeline;
 
class Discretizer
{
public:
    virtual ~Discretizer();
 
    Discretizer(
            const arma::uword nGridPoints,
            const arma::uword nEquationsAndVariables);
 
    virtual void discretize(
            const arma::uword dt,
            const Pipeline& currentState,
            const Pipeline& newState) = 0;
 
    const arma::cube& term_i() const { return m_term_i; } 
    const arma::cube& term_ipp() const { return m_term_ipp; } 
    const arma::mat& boundaryTerms() const { return m_boundaryTerm; } 
 
protected:
    // TODO: better description of how m_term_i and m_term_ipp are organized
    arma::cube m_term_i;
 
    arma::cube m_term_ipp;
 
    arma::mat m_boundaryTerm;
 
    double m_gravity = 9.81; 
};