ohms_law_cyl_source.h

Go to the documentation of this file.

#pragma once
 
#include "apps/wmapplication.h"
 
namespace wxm
{
namespace apps
{
namespace hallmhd
{
namespace hyperresistivity
{
 
class OhmsLawCylSource : public WmApplication
{
public:
    OhmsLawCylSource();
 
    ~OhmsLawCylSource() override;
 
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getInputVariableIndexes(int flag) const override
    {
        return _input_variable_idcs;
    }
 
    const std::vector<int>& getAuxiliaryVariableIndexes(int flag) const override
    {
        return _aux_variable_idcs;
    }
 
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return _output_variable_idcs;
    }
 
    real source(const real* q,
                const real* aux,
                const elementGeometry_t* pEG,
                real* source) const override;
 
protected:
    real _hyperresistive_nu;
 
    int _radial_index;
    int _ir, _ith, _iz; // r, theta, and z indices
 
    std::vector<int> _input_variable_idcs;
    std::vector<int> _aux_variable_idcs;
    std::vector<int> _output_variable_idcs;
 
private:
    OhmsLawCylSource& operator=(const OhmsLawCylSource& var) = delete;
    OhmsLawCylSource(const OhmsLawCylSource& var) = delete;
};
} // namespace hyperresistivity
} // namespace hallmhd
} // namespace apps
} // namespace wxm