Method Init

Source/Parallelization/GuardCellManager.cpp:35–393

using namespace amrex;

void
guardCellManager::Init (
    const amrex::Real dt,
    const amrex::Real *dx,
    const bool do_subcycling,
    const bool do_fdtd_nci_corr,
    ablastr::utils::enums::GridType grid_type,
    const bool do_moving_window,
    const int moving_window_dir,
    const int particle_max_grid_crossings,
    const int nox,
    const int nox_fft, const int noy_fft, const int noz_fft,
    const int nci_corr_stencil,
    const ElectromagneticSolverAlgo electromagnetic_solver_id,
    const EvolveScheme evolve_scheme,
    const int max_level,
    const amrex::Vector<amrex::Real>& v_galilean,
    const amrex::Vector<amrex::Real>& v_comoving,
    const bool safe_guard_cells,
    const int do_psatd_JRhom,
    const bool fft_do_time_averaging,
    const bool do_pml,
    const int do_pml_in_domain,
    const int pml_ncell,
    const amrex::Vector<amrex::IntVect>& ref_ratios,
    const bool use_filter,
    const amrex::IntVect& bilinear_filter_stencil_length)
{
    // When using subcycling, the particles on the fine level perform two pushes
    // before being redistributed ; therefore, we need one extra guard cell
    // (the particles may move by 2*c*dt)
    int ngx_tmp = (max_level > 0 && do_subcycling) ? nox+1 : nox;
    int ngy_tmp = (max_level > 0 && do_subcycling) ? nox+1 : nox;
    int ngz_tmp = (max_level > 0 && do_subcycling) ? nox+1 : nox;

    const bool galilean = (v_galilean[0] != 0. || v_galilean[1] != 0. || v_galilean[2] != 0.);
    const bool comoving = (v_comoving[0] != 0. || v_comoving[1] != 0. || v_comoving[2] != 0.);

    // Add one guard cell in the case of the Galilean or comoving algorithms
    if (galilean || comoving)
    {
      ngx_tmp += 1;
      ngy_tmp += 1;
      ngz_tmp += 1;
    }

    // Ex, Ey, Ez, Bx, By, and Bz have the same number of ghost cells.
    // jx, jy, jz and rho have the same number of ghost cells.
    // E and B have the same number of ghost cells as j and rho if NCI filter is not used,
    // but different number of ghost cells in z-direction if NCI filter is used.
    // The number of cells should be even, in order to easily perform the
    // interpolation from coarse grid to fine grid.
    int ngx = (ngx_tmp % 2) ? ngx_tmp+1 : ngx_tmp;  // Always even number
    int ngy = (ngy_tmp % 2) ? ngy_tmp+1 : ngy_tmp;  // Always even number
    const int ngz_nonci = (ngz_tmp % 2) ? ngz_tmp+1 : ngz_tmp;  // Always even number
    int ngz;
    if (do_fdtd_nci_corr) {
        const int ng = ngz_tmp + nci_corr_stencil;