From ef7e2ae56cb62da09c2602488398137c3d8a90f4 Mon Sep 17 00:00:00 2001 From: Edoardo Zoni <59625522+EZoni@users.noreply.github.com> Date: Thu, 3 Mar 2022 15:14:54 -0800 Subject: Refactor Vay Deposition Functions (#2885) * Refactor Vay Deposition Functions * Use rho Index Type for FFT of D * Improve Doxygen Strings --- .../SpectralAlgorithms/PsatdAlgorithm.cpp | 19 +------------------ 1 file changed, 1 insertion(+), 18 deletions(-) (limited to 'Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithm.cpp') diff --git a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithm.cpp b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithm.cpp index 52813cb13..35429109a 100644 --- a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithm.cpp +++ b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithm.cpp @@ -832,25 +832,13 @@ void PsatdAlgorithm::CurrentCorrection (SpectralFieldData& field_data) } void -PsatdAlgorithm::VayDeposition ( - const int lev, - SpectralFieldData& field_data, - std::array,3>& current) +PsatdAlgorithm::VayDeposition (SpectralFieldData& field_data) { // Profiling BL_PROFILE("PsatdAlgorithm::VayDeposition()"); const SpectralFieldIndex& Idx = m_spectral_index; - // Forward Fourier transform of D (temporarily stored in current): - // D is nodal and does not match the staggering of J, therefore we pass the - // actual staggering of D (IntVect(1)) to the ForwardTransform function - field_data.ForwardTransform(lev, *current[0], Idx.Jx, 0, IntVect(1)); - field_data.ForwardTransform(lev, *current[1], Idx.Jy, 0, IntVect(1)); - field_data.ForwardTransform(lev, *current[2], Idx.Jz, 0, IntVect(1)); - - const amrex::IntVect& fill_guards = m_fill_guards; - // Loop over boxes for (amrex::MFIter mfi(field_data.fields); mfi.isValid(); ++mfi) { @@ -903,11 +891,6 @@ PsatdAlgorithm::VayDeposition ( else fields(i,j,k,Idx.Jz) = 0._rt; }); } - - // Backward Fourier transform of J - field_data.BackwardTransform(lev, *current[0], Idx.Jx, 0, fill_guards); - field_data.BackwardTransform(lev, *current[1], Idx.Jy, 0, fill_guards); - field_data.BackwardTransform(lev, *current[2], Idx.Jz, 0, fill_guards); } #endif // WARPX_USE_PSATD -- cgit v1.2.3