pulses_fluid5.h

Go to the documentation of this file.

#ifndef wmapplication_function_five_moment_pulses_h
#define wmapplication_function_five_moment_pulses_h
 
// WARPXM includes
#include "apps/wmapplication.h"
 
namespace wxm
{
namespace apps
{
namespace functions
{
namespace five_moment
{
 
class PulsesFluid5 : public WmApplication
{
public:
    PulsesFluid5();
 
    ~PulsesFluid5() override;
 
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return _output_fluid;
    }
 
    void evaluate_function(const real* q,
                           const real* aux,
                           const solverVariables_t* pSV,
                           real* result) const override;
 
protected:
    std::vector<int> _output_fluid;
 
    real calcRadiusPosition(const real* xc, const real* x) const;
    bool isInPulse(const real& a, const real& r) const;
    real gaussianPulse(const real& a, const real& r) const;
    real polynomialPulse(const real& a, const real& r) const;
 
    real _mass;
    real _n0;
    real _n1;
    real _temperature;
    real _gas_gamma;
 
    real _radius;
    real _num_cols;
    bool _gaussian;
    bool _polynomial;
    real _k;
 
    std::vector<real> _center_x;
    std::vector<real> _center_y;
    std::vector<real> _center_z;
 
private:
    PulsesFluid5& operator=(const PulsesFluid5& var);
    PulsesFluid5(const PulsesFluid5& var);
};
} // namespace five_moment
} // namespace functions
} // namespace apps
} // namespace wxm
 
#endif