Maxwell.h

//
// Class Maxwell
//   Base class for solvers for Maxwell's equations
//
 
#ifndef IPPL_MAXWELL_H
#define IPPL_MAXWELL_H
 
#include "Types/Vector.h"
 
#include "Field/Field.h"
 
#include "Utility/ParameterList.h"
#include "FieldLayout/FieldLayout.h"
#include "Meshes/UniformCartesian.h"
 
namespace ippl {
 
    template <typename EMField, typename SourceField>
    class Maxwell {
    public:
        constexpr static unsigned Dim = EMField::dim;
 
        Maxwell() {}
 
        Maxwell(SourceField& four_current, EMField& E, EMField& B) {
            setSources(four_current);
            setEMFields(E, B);
        }
 
        virtual void setSources(SourceField& four_current) { JN_mp = &four_current; }
 
        void setEMFields(EMField& E, EMField& B) {
            En_mp = &E;
            Bn_mp = &B;
        }
 
 
        void mergeParameters(const ParameterList& params) { params_m.merge(params); }
 
        virtual void solve() = 0;
 
        virtual ~Maxwell() {}
 
    protected:
        // Parameters
        ParameterList params_m;
 
        // Field for four-current (rho, J)
        SourceField* JN_mp = nullptr;
 
        // E and B fields
        EMField* En_mp = nullptr;
        EMField* Bn_mp = nullptr;
    };
}  // namespace ippl
 
#endif