divergence_error_calculation.h

Go to the documentation of this file.

#pragma once
 
// WARPXM includes
#include "apps/wmapplication.h"
 
namespace wxm
{
namespace apps
{
namespace functions
{
namespace maxwell
{
 
class DivergenceError : public WmApplication
{
public:
    DivergenceError();
 
    ~DivergenceError();
 
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getInputVariableIndexes(int flag) const override
    {
        return _input_indices;
    }
 
    const std::vector<int>& getAuxiliaryVariableIndexes(int flag) const override
    {
        return _auxiliary_fluid_indices;
    }
 
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return _output_divergence_error_indices;
    }
 
    void evaluate_function(const real* q,
                           const real* aux,
                           const solverVariables_t* pSV,
                           real* result) const override;
 
protected:
    real cartesian_divergence(const real* field_gradients) const;
    real cylindrical_divergence(const real* fields,
                                const real* field_gradients,
                                const solverVariables_t* pSV) const;
 
    std::vector<int> _input_indices;
    std::vector<int> _auxiliary_fluid_indices;
    std::vector<int> _output_divergence_error_indices;
 
    real _omega_p_tau;
    real _skin_depth_norm;
 
    int _num_fluids;
    std::vector<real> _charges;
    std::vector<real> _masses;
    std::vector<real> _com;
 
    real _factor;
 
    int _coordinate_system;
    int _radial_index;
    int _ir, _ith, _iz; // r, theta, and z indices
 
private:
    DivergenceError& operator=(const DivergenceError& var);
    DivergenceError(const DivergenceError& var);
};
 
} // namespace maxwell
} // namespace functions
} // namespace apps
} // namespace wxm