WarpXSumGuardCells.H

Go to the documentation of this file.

/* Copyright 2019 Maxence Thevenet, Remi Lehe, Weiqun Zhang
 *
 *
 * This file is part of WarpX.
 *
 * License: BSD-3-Clause-LBNL
 */
#ifndef WARPX_SUM_GUARD_CELLS_H_
#define WARPX_SUM_GUARD_CELLS_H_

#include <AMReX_MultiFab.H>

inline void
WarpXSumGuardCells(amrex::MultiFab& mf, const amrex::Periodicity& period,
                   const int icomp=0, const int ncomp=1){
#ifdef WARPX_USE_PSATD
   // Update both valid cells and guard cells
   const amrex::IntVect n_updated_guards = mf.nGrowVect();
#else
   // Update only the valid cells
   const amrex::IntVect n_updated_guards = amrex::IntVect::TheZeroVector();
#endif
    mf.SumBoundary(icomp, ncomp, n_updated_guards, period);
}

inline void
WarpXSumGuardCells(amrex::MultiFab& dst, amrex::MultiFab& src,
                   const amrex::Periodicity& period,
                   const int icomp=0, const int ncomp=1){
#ifdef WARPX_USE_PSATD
    // Update both valid cells and guard cells
    const amrex::IntVect n_updated_guards = dst.nGrowVect();
#else
    // Update only the valid cells
    const amrex::IntVect n_updated_guards = amrex::IntVect::TheZeroVector();
#endif
    src.SumBoundary(0, ncomp, n_updated_guards, period);
    amrex::Copy( dst, src, 0, icomp, ncomp, n_updated_guards );
}

#endif // WARPX_SUM_GUARD_CELLS_H_