diff options
Diffstat (limited to 'Source/Utils/CoarsenMR.cpp')
-rw-r--r-- | Source/Utils/CoarsenMR.cpp | 104 |
1 files changed, 0 insertions, 104 deletions
diff --git a/Source/Utils/CoarsenMR.cpp b/Source/Utils/CoarsenMR.cpp deleted file mode 100644 index 549ff6ea2..000000000 --- a/Source/Utils/CoarsenMR.cpp +++ /dev/null @@ -1,104 +0,0 @@ -#include "CoarsenMR.H" - -#include "Utils/TextMsg.H" - -#include <AMReX_BLProfiler.H> -#include <AMReX_BLassert.H> -#include <AMReX_BoxArray.H> -#include <AMReX_Config.H> -#include <AMReX_GpuControl.H> -#include <AMReX_GpuLaunch.H> -#include <AMReX_IndexType.H> -#include <AMReX_IntVect.H> -#include <AMReX_MFIter.H> -#include <AMReX_MultiFab.H> - -using namespace amrex; - -void -CoarsenMR::Loop ( MultiFab& mf_dst, - const MultiFab& mf_src, - const int ncomp, - const IntVect ngrow, - const IntVect crse_ratio ) -{ - // Staggering of source fine MultiFab and destination coarse MultiFab - const IntVect stag_src = mf_src.boxArray().ixType().toIntVect(); - const IntVect stag_dst = mf_dst.boxArray().ixType().toIntVect(); - - // Auxiliary integer arrays (always 3D) - GpuArray<int,3> sf; // staggering of source fine MultiFab - GpuArray<int,3> sc; // staggering of destination coarse MultiFab - GpuArray<int,3> cr; // coarsening ratio - - sf[0] = stag_src[0]; -#if defined(WARPX_DIM_1D_Z) - sf[1] = 0; -#else - sf[1] = stag_src[1]; -#endif -#if (AMREX_SPACEDIM <= 2) - sf[2] = 0; -#elif defined(WARPX_DIM_3D) - sf[2] = stag_src[2]; -#endif - - sc[0] = stag_dst[0]; -#if defined(WARPX_DIM_1D_Z) - sc[1] = 0; -#else - sc[1] = stag_dst[1]; -#endif -#if (AMREX_SPACEDIM <= 2) - sc[2] = 0; -#elif defined(WARPX_DIM_3D) - sc[2] = stag_dst[2]; -#endif - - cr[0] = crse_ratio[0]; -#if defined(WARPX_DIM_1D_Z) - cr[1] = 1; -#else - cr[1] = crse_ratio[1]; -#endif -#if (AMREX_SPACEDIM <= 2) - cr[2] = 1; -#elif defined(WARPX_DIM_3D) - cr[2] = crse_ratio[2]; -#endif - -#ifdef AMREX_USE_OMP -#pragma omp parallel if (Gpu::notInLaunchRegion()) -#endif - // Loop over boxes (or tiles if not on GPU) - for (MFIter mfi( mf_dst, TilingIfNotGPU() ); mfi.isValid(); ++mfi) - { - // Tiles defined at the coarse level - const Box& bx = mfi.growntilebox( ngrow ); - Array4<Real> const& arr_dst = mf_dst.array( mfi ); - Array4<Real const> const& arr_src = mf_src.const_array( mfi ); - ParallelFor( bx, ncomp, - [=] AMREX_GPU_DEVICE( int i, int j, int k, int n ) - { - arr_dst(i,j,k,n) = CoarsenMR::Interp( - arr_src, sf, sc, cr, i, j, k, n ); - } ); - } -} - -void -CoarsenMR::Coarsen ( MultiFab& mf_dst, - const MultiFab& mf_src, - const IntVect crse_ratio ) -{ - BL_PROFILE("CoarsenMR::Coarsen()"); - - WARPX_ALWAYS_ASSERT_WITH_MESSAGE( mf_src.ixType() == mf_dst.ixType(), - "source MultiFab and destination MultiFab have different IndexType" ); - - // Number of guard cells to fill on coarse patch and number of components - const IntVect ngrow = ( mf_src.nGrowVect() + 1 ) / crse_ratio; - const int ncomp = mf_src.nComp(); - - CoarsenMR::Loop( mf_dst, mf_src, ncomp, ngrow, crse_ratio ); -} |