wmapplication_axisymmetric_imhd_1D.h

Go to the documentation of this file.

#ifndef WMAPPLICATION_AXISYMMETRIC_IMHD_1D_H
#define WMAPPLICATION_AXISYMMETRIC_IMHD_1D_H
 
// wxm includes
#include "apps/wmapplication.h"
 
class WmApplication_AXISYMMETRIC_IMHD_1D : public WmApplication
{
public:
    WmApplication_AXISYMMETRIC_IMHD_1D();
 
    ~WmApplication_AXISYMMETRIC_IMHD_1D() override;
 
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getInputVariableIndexes(int flag) const override
    {
        return _inputVariables;
    }
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return _outputVariables;
    }
 
    // const std::vector<int>& getAuxiliaryVariableIndexes(int flag) const override
    // {
    //     return _aux_variables;
    // }
 
    // redefinition of fluxes from parent
    real numerical_flux(const real* q_l,
                        const real* q_r,
                        const real* aux_l,
                        const real* aux_r,
                        const solverVariables_t* pFV,
                        real* numericalFlux) const override;
    real internal_flux(const real* q,
                       const real* aux,
                       const solverVariables_t* pSV,
                       std::vector<std::vector<real>>& internalFlux) const override;
 
    // redefinition of geometric source terms from parent
    real source(const real* q,
                const real* aux,
                const elementGeometry_t* pEG,
                real* source) const override;
 
protected:
    real _gas_gamma;
 
    std::string _numerical_flux_type; // this is the type of numerical flux you want to
                                      // use (rusanov, roe, hll, hlld, etc.)
 
    std::vector<int> _inputVariables;
    std::vector<int> _outputVariables;
    // std::vector<int> _aux_variables;
 
    // real _c_h, _c_p, _c_r;
 
    real _min_density_floor;
    real _min_pressure_floor;
 
private:
    WmApplication_AXISYMMETRIC_IMHD_1D&
    operator=(const WmApplication_AXISYMMETRIC_IMHD_1D& var);
    WmApplication_AXISYMMETRIC_IMHD_1D(const WmApplication_AXISYMMETRIC_IMHD_1D& var);
};
 
#endif // WMAPPLICATION_AXISYMMETRIC_IMHD_1D_H