From 47f91938af6c92f36d234738d267ac18e00842a6 Mon Sep 17 00:00:00 2001 From: Remi Lehe Date: Fri, 8 May 2020 09:41:22 -0700 Subject: Fix bug for PSATD momentum-conserving runs (#965) * Relax PSATD check on boxes, so that PSATD momentum-conserving can be used * Increase tolerance of test * Add automated tests for momentum-conserving PSATD --- Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp') diff --git a/Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp b/Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp index 68bc95f53..76299b7de 100644 --- a/Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp +++ b/Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp @@ -233,10 +233,10 @@ SpectralFieldData::ForwardTransform( const MultiFab& mf, realspace_bx = mf[mfi].box(); // Keep guard cells } realspace_bx.enclosedCells(); // Discard last point in nodal direction - AMREX_ALWAYS_ASSERT( realspace_bx == tmpRealField[mfi].box() ); + AMREX_ALWAYS_ASSERT( realspace_bx.contains(tmpRealField[mfi].box()) ); Array4 mf_arr = mf[mfi].array(); Array4 tmp_arr = tmpRealField[mfi].array(); - ParallelFor( realspace_bx, + ParallelFor( tmpRealField[mfi].box(), [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept { tmp_arr(i,j,k) = mf_arr(i,j,k,i_comp); }); -- cgit v1.2.3