combustionprocess.h

#ifndef zh2_combustionprocess
#define zh2_combustionprocess

#include "integrable.h"

#include <vector>
#include <unordered_map>

namespace h2comb {

class Reaction;
class Compound;
class Substance;

class CombustionProcess : public Integrable
{
public:
    CombustionProcess(Compound *c, std::ostream *o);// : compound(c), output(o) {}
    
    virtual int getDim();
    
    void diffEq(const double *const x, const double t, double *const dx);
    
    double getState(double *const x); 
    void setState(const double *const x, const double t);
    virtual void peek(const double *const, const double);
    
    void listSubstances();
    
    void setAdiabatic(){adiabatic = true;}
    
protected:
    std::unordered_map<Substance *, int> indices;
    int numberOfSubstances;
    
private:
    Compound *compound;
    
    bool frozen = false;
    
    std::vector<Substance *> substances;
    int* quantities;
    
    std::ostream *output;
    double saytime;
    bool adiabatic;
};

} // namespace

#endif