StrangImplicitSpectralEM.H

Go to the documentation of this file.

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

 
class StrangImplicitSpectralEM : public ImplicitSolver
{
public:
 
    StrangImplicitSpectralEM() = default;
 
    ~StrangImplicitSpectralEM() override = default;
 
    // Prohibit Move and Copy operations
    StrangImplicitSpectralEM(const StrangImplicitSpectralEM&) = delete;
    StrangImplicitSpectralEM& operator=(const StrangImplicitSpectralEM&) = delete;
    StrangImplicitSpectralEM(StrangImplicitSpectralEM&&) = delete;
    StrangImplicitSpectralEM& operator=(StrangImplicitSpectralEM&&) = delete;
 
    void Define ( WarpX*  a_WarpX ) override;
 
    void PrintParameters () const override;
 
    void 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;

    amrex::Vector<std::array< std::unique_ptr<amrex::MultiFab>, 3 > > m_Bold;

    void UpdateWarpXFields ( WarpXSolverVec const&  a_E,
                             amrex::Real            half_time );

    void FinishFieldUpdate ( amrex::Real  a_new_time );
 
};
 
#endif