current_density.h

Go to the documentation of this file.

#ifndef WXM_APPS_MHD_CURRENT_DENSITY_H
#define WXM_APPS_MHD_CURRENT_DENSITY_H
 
// General includes
#include "apps/wmapplication.h"
 
namespace wxm
{
namespace apps
{
namespace rmhd
{
 
class CurrentDensity : public WmApplication
{
public:
    CurrentDensity();
 
    ~CurrentDensity() override;
 
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getInputVariableIndexes(int flag) const override
    {
        return _grad_B;
    }
 
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return _J;
    }
 
    // const std::vector<int>& getAuxiliaryVariableIndexes(int flag) const override
    // {
    //     return _auxVariables;
    // }
 
    void evaluate_function(const real* q,
                           const real* aux,
                           const solverVariables_t* pSV,
                           real* result) const override;
 
protected:
    std::vector<int> _grad_B; // conserved variable vector to output to
    std::vector<int> _J;      // conserved variable vector to output to
    real _skin_depth_norm;
 
    // if cylindrical curlB has to be calculated differently
    std::string _coordinate_system;
    int _radial_index;
    int _ir, _ith, _iz; // r, theta, and z indices
 
private:
    CurrentDensity& operator=(const CurrentDensity& var);
    CurrentDensity(const CurrentDensity& var);
};
 
} // namespace rmhd
} // namespace apps
} // namespace wxm
 
#endif // WXM_APPS_MHD_CURRENT_DENSITY_H