pcmaxwell.h

Go to the documentation of this file.

#ifndef WXM_APPS_PCMAXWELL_H
#define WXM_APPS_PCMAXWELL_H
 
// WarpM includes
#include "apps/wmapplication.h"
 
namespace wxm
{
namespace apps
{
namespace pcmaxwell
{
 
class pcmaxwell_t : public WmApplication
{
public:
    pcmaxwell_t();
 
    ~pcmaxwell_t() override;
 
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getInputVariableIndexes(int flag) const override
    {
        return _variables;
    }
    const std::vector<int>& getAuxiliaryVariableIndexes(int flag) const override
    {
        return _gradient_variables;
    }
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return _variables;
    }
 
    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;
 
protected:
    std::string _coordinate_system;
    std::string _gradient_method;
    real _chi_e;
    real _chi_b;
 
    int _radial_index;
    int _ir, _ith, _iz; // r, theta, and z indices
 
    std::vector<int> _variables;
    std::vector<int> _gradient_variables;
 
    // flux tensor is used to calculate left and right fluxes for numerical flux
    // mutable because this will continuously be changed upon call, but functions are
    // labled const this is just a helper vector this behavior should be justified.
    mutable std::vector<std::vector<real>> _flux_tensor;
 
    // specific calls
    real internal_flux_cartesian(const real* q,
                                 const real* aux,
                                 const solverVariables_t* pSV,
                                 std::vector<std::vector<real>>& internalFlux) const;
 
    real internal_flux_cylindrical(const real* q,
                                   const real* aux,
                                   const solverVariables_t* pSV,
                                   std::vector<std::vector<real>>& internalFlux) const;
 
    real numerical_flux_ldg(const real* q_l,
                            const real* q_r,
                            const real* aux_l,
                            const real* aux_r,
                            const solverVariables_t* pFV,
                            real* numericalFlux) const;
 
private:
    pcmaxwell_t& operator=(const pcmaxwell_t& var);
    pcmaxwell_t(const pcmaxwell_t& var);
};
} // namespace pcmaxwell
} // namespace apps
} // namespace wxm
 
#endif // WXM_APPS_PCMAXWELL_H