eos.h

Go to the documentation of this file.

#ifndef __EOS_H__
#define __EOS_H__

#include <vector>


class EOS {
protected:
    int m_iEOSChoice; 
public:
    virtual ~EOS() {};

    int getEOSChoice() {return m_iEOSChoice;}
    
    virtual void getParameters(std::vector<double>& params) = 0;

    virtual double getEnergy(double pressure, double density) = 0;
    
    virtual double getSoundSpeed(double pressure, double density) = 0;
};










class PolytropicGasEOS : public EOS {
protected:
    double m_fGamma; 

public:
    PolytropicGasEOS(double gamma) : m_fGamma(gamma) {m_iEOSChoice=1;}
    
    // Destructor
    virtual ~PolytropicGasEOS() {}  

    virtual void getParameters(std::vector<double>& params) {params.push_back(m_fGamma);}

    virtual double getEnergy(double pressure, double density);

    virtual double getSoundSpeed(double pressure, double density);
};












class StiffPolytropicGasEOS : public EOS {
protected:
    double m_fGamma; 
    double m_fPinf; 
    double m_fEinf; 

public:
    StiffPolytropicGasEOS(double gamma, double pinf, double einf):
        m_fGamma(gamma), m_fPinf(pinf), m_fEinf(einf) {m_iEOSChoice=2;}
    
    virtual ~StiffPolytropicGasEOS() {}
    
    virtual void getParameters(std::vector<double>& params) { 
        params.push_back(m_fGamma);
        params.push_back(m_fPinf);
        params.push_back(m_fEinf);
    }

    virtual double getEnergy(double pressure, double density);
    
    virtual double getSoundSpeed(double pressure, double density);
};



#endif // __EOS_H__