NonStandardFDTDSolver.h

Go to the documentation of this file.

 
#ifndef IPPL_NON_STANDARD_FDTD_SOLVER_H
#define IPPL_NON_STANDARD_FDTD_SOLVER_H
 
#include <cstddef>
using std::size_t;
#include "Types/Vector.h"
 
#include "FieldLayout/FieldLayout.h"
#include "MaxwellSolvers/AbsorbingBC.h"
#include "MaxwellSolvers/FDTDSolverBase.h"
#include "MaxwellSolvers/Maxwell.h"
#include "Meshes/UniformCartesian.h"
#include "Particle/ParticleBase.h"
 
namespace ippl {
 
    template <typename EMField, typename SourceField, fdtd_bc boundary_conditions = periodic>
    class NonStandardFDTDSolver : public FDTDSolverBase<EMField, SourceField, boundary_conditions> {
    public:
        NonStandardFDTDSolver(SourceField& source, EMField& E, EMField& B);
 
        constexpr static unsigned Dim = EMField::dim;  // Dimension of the electromagnetic field.
        using scalar = typename EMField::value_type::value_type;  // Scalar type used in the
                                                                  // electromagnetic field.
        using Vector_t = Vector<typename EMField::value_type::value_type,
                                Dim>;  // Vector type used in the electromagnetic field.
        using SourceVector_t =
            typename SourceField::value_type;  // Vector type used in the source field.
 
        template <typename scalar>
        struct nondispersive {
            scalar a1;
            scalar a2;
            scalar a4;
            scalar a6;
            scalar a8;
        };
 
        void step() override;
        void initialize() override;
    };
}  // namespace ippl
 
#include "NonStandardFDTDSolver.hpp"
#endif