EmbeddedBoundaryInit.H

Go to the documentation of this file.

/* Copyright 2021-2025 Lorenzo Giacomel, Luca Fedeli
 *
 * This file is part of WarpX.
 *
 * License: BSD-3-Clause-LBNL
 */
 
#ifndef WARPX_EMBEDDED_BOUNDARY_EMBEDDED_BOUNDARY_INIT_H_
#define WARPX_EMBEDDED_BOUNDARY_EMBEDDED_BOUNDARY_INIT_H_
 
#include "Enabled.H"
 
#ifdef AMREX_USE_EB
 
#include <ablastr/fields/MultiFabRegister.H>
 
#include <AMReX_EBFabFactory.H>
#include <AMReX_iMultiFab.H>
#include <AMReX_Periodicity.H>
#include <AMReX_REAL.H>
 
#include <array>
 
namespace warpx::embedded_boundary
{

    void MarkReducedShapeCells (
        std::unique_ptr<amrex::iMultiFab> & eb_reduce_particle_shape,
        amrex::EBFArrayBoxFactory const & eb_fact,
        int particle_shape_order,
        const amrex::Periodicity& periodicity);

    void MarkUpdateCellsStairCase (
        std::array< std::unique_ptr<amrex::iMultiFab>,3> & eb_update,
        ablastr::fields::VectorField const & field,
        amrex::EBFArrayBoxFactory const & eb_fact );

    void MarkUpdateECellsECT (
        std::array< std::unique_ptr<amrex::iMultiFab>,3> & eb_update_E,
        ablastr::fields::VectorField const& edge_lengths );

    void MarkUpdateBCellsECT (
        std::array< std::unique_ptr<amrex::iMultiFab>,3> & eb_update_B,
        ablastr::fields::VectorField const& face_areas,
        ablastr::fields::VectorField const& edge_lengths );

    void MarkExtensionCells(
        const std::array<amrex::Real,3>& cell_size,
        std::array< std::unique_ptr<amrex::iMultiFab>, 3 > & flag_info_face,
        std::array< std::unique_ptr<amrex::iMultiFab>, 3 > & flag_ext_face,
        const ablastr::fields::VectorField& b_field,
        const ablastr::fields::VectorField& face_areas,
        const ablastr::fields::VectorField& edge_lengths,
        const ablastr::fields::VectorField& area_mod);

    void ComputeEdgeLengths (
        ablastr::fields::VectorField& edge_lengths,
        const amrex::EBFArrayBoxFactory& eb_fact);
    void ComputeFaceAreas (
        ablastr::fields::VectorField& face_areas,
        const amrex::EBFArrayBoxFactory& eb_fact);

    void ScaleEdges (
        ablastr::fields::VectorField& edge_lengths,
        const std::array<amrex::Real,3>& cell_size);
    void ScaleAreas (
        ablastr::fields::VectorField& face_areas,
        const std::array<amrex::Real,3>& cell_size);
}
 
#endif
 
#endif //WARPX_EMBEDDED_BOUNDARY_EMBEDDED_BOUNDARY_H_