MagnetostaticSolver.H

Go to the documentation of this file.

/* Copyright 2022 S. Eric Clark, LLNL
 *
 * This file is part of WarpX.
 *
 * License: BSD-3-Clause-LBNL
 */
#ifndef WARPX_MAGNETOSTATICSOLVER_H_
#define WARPX_MAGNETOSTATICSOLVER_H_
 
#include <ablastr/fields/MultiFabRegister.H>
 
#include <AMReX_Array.H>
#include <AMReX_MultiFab.H>
#include <AMReX_MLMG.H>
#include <AMReX_REAL.H>
 
namespace MagnetostaticSolver {

    class VectorPoissonBoundaryHandler {
    public:
        amrex::Array<amrex::Array<amrex::LinOpBCType, AMREX_SPACEDIM>, 3> lobc, hibc;
        bool bcs_set = false;
        std::array<std::array<bool, AMREX_SPACEDIM * 2>, 3> dirichlet_flag;
        bool has_non_periodic = false;
 
        void defineVectorPotentialBCs ();
    };

    class EBCalcBfromVectorPotentialPerLevel {
      private:
        ablastr::fields::MultiLevelVectorField m_b_field;
        ablastr::fields::MultiLevelVectorField m_grad_buf_e_stag;
        ablastr::fields::MultiLevelVectorField m_grad_buf_b_stag;
 
      public:
        EBCalcBfromVectorPotentialPerLevel (ablastr::fields::MultiLevelVectorField const & b_field,
                                            ablastr::fields::MultiLevelVectorField const & grad_buf_e_stag,
                                            ablastr::fields::MultiLevelVectorField const & grad_buf_b_stag)
                : m_b_field(b_field),
                  m_grad_buf_e_stag(grad_buf_e_stag),
                  m_grad_buf_b_stag(grad_buf_b_stag)
        {}
 
        void operator() (amrex::Array<std::unique_ptr<amrex::MLMG>, 3> & mlmg, int lev);
 
        // Function to perform interpolation from cell edges to cell faces
        void doInterp (amrex::MultiFab & src, amrex::MultiFab & dst);
    };
} // namespace MagnetostaticSolver
 
#endif //WARPX_MAGNETOSTATICSOLVER_H_