SemiImplicitEM.H

Go to the documentation of this file.

/* Copyright 2024 Justin Angus
 *
 * This file is part of WarpX.
 *
 * License: BSD-3-Clause-LBNL
 */
#ifndef SEMI_IMPLICIT_EM_H_
#define SEMI_IMPLICIT_EM_H_
 
#include "FieldSolver/ImplicitSolvers/WarpXSolverVec.H"
 
#include <AMReX_Array.H>
#include <AMReX_MultiFab.H>
#include <AMReX_REAL.H>
 
#include "ImplicitSolver.H"

 
class SemiImplicitEM : public ImplicitSolver
{
public:
 
    SemiImplicitEM() = default;
 
    ~SemiImplicitEM() override = default;
 
    // Prohibit Move and Copy operations
    SemiImplicitEM(const SemiImplicitEM&) = delete;
    SemiImplicitEM& operator=(const SemiImplicitEM&) = delete;
    SemiImplicitEM(SemiImplicitEM&&) = delete;
    SemiImplicitEM& operator=(SemiImplicitEM&&) = delete;
 
    void Define (WarpX* a_WarpX, bool a_from_restart) override;
 
    void PrintParameters () const override;
 
    int OneStep (amrex::Real  start_time,
                 amrex::Real  a_dt,
                 int          a_step) override;
 
    void ComputeRHS ( WarpXSolverVec&  a_RHS,
                const WarpXSolverVec&  a_E,
                      amrex::Real      start_time,
                      int              a_nl_iter,
                      bool             a_from_jacobian ) override;
 
    // This parameter is used for the time-step fraction in the PC for implicit
    // treatment of light waves in the curl-curl MLMG solver.
    // This function should return zero if light waves are not treated implicitly
    [[nodiscard]] amrex::Real  GetThetaForPC () const override { return 0.; }
 
private:

    WarpXSolverVec m_E, m_Eold;
 
};
 
#endif