From 050e12f5293d42aabd812ff314baa2199060cdc3 Mon Sep 17 00:00:00 2001 From: Olga Shapoval <30510597+oshapoval@users.noreply.github.com> Date: Tue, 18 Feb 2020 22:25:34 -0800 Subject: Galilean PSATD with shift (#704) * Read Galilean velocity * Prepare structures for Galilean solver * Started implementing Galilean equations * Analytical limits for X1, X2, X3, X4 coefficients added * Slight changes added * Added Galilean position pusher * Scale galilean velocity * Remove unneeded Abort * Fix Galilean pusher * Allocate Theta2 array * Fix definition of coefficients * Increase guard cells for Galilean * Add guard cell in particle exchange * Type corrected * v_gal added to warpx_current_deposition * v_gal added to WarpXParticleContainer.H * Bug fixed - update particle x-position over one time step * Fix issues with merge from dev * Preparation for merging dev into galilean. * Adding galilean shift * Implemented galilean shift * Changed method's name from GalileanShift to ShiftGalileanBoundary * Added doxygen string for ShiftGalileanBoundary * Removed never used method LowerCornerWithCentering * Removed temporary comments * Removed dt as a variable from DepositCharge method and its dependencies * Converted tab to spaces * Removed EOL white space * Add documentation and automated tests * Fix compilation error * Add automated test * Update automated test * Removed temporary used galilean shift * Removed temporary used particle's push for Galilean PSATD * Removed unused statement * Remove EOL white space. * Added zero shift for LowerCorner in RZ geometry * Minor changes to Galilean implementation * Modifications for GPU * Fix typo Co-authored-by: Remi Lehe --- Source/Parallelization/GuardCellManager.cpp | 16 +++++++++++++--- 1 file changed, 13 insertions(+), 3 deletions(-) (limited to 'Source/Parallelization/GuardCellManager.cpp') diff --git a/Source/Parallelization/GuardCellManager.cpp b/Source/Parallelization/GuardCellManager.cpp index d845a7447..d7d2bce7c 100644 --- a/Source/Parallelization/GuardCellManager.cpp +++ b/Source/Parallelization/GuardCellManager.cpp @@ -24,14 +24,24 @@ guardCellManager::Init( const int nci_corr_stencil, const int maxwell_fdtd_solver_id, const int max_level, + const amrex::Array v_galilean, const bool exchange_all_guard_cells) { // 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) - const int ngx_tmp = (max_level > 0 && do_subcycling == 1) ? nox+1 : nox; - const int ngy_tmp = (max_level > 0 && do_subcycling == 1) ? nox+1 : nox; - const int ngz_tmp = (max_level > 0 && do_subcycling == 1) ? nox+1 : nox; + int ngx_tmp = (max_level > 0 && do_subcycling == 1) ? nox+1 : nox; + int ngy_tmp = (max_level > 0 && do_subcycling == 1) ? nox+1 : nox; + int ngz_tmp = (max_level > 0 && do_subcycling == 1) ? nox+1 : nox; + + if ((v_galilean[0]!=0) or + (v_galilean[1]!=0) or + (v_galilean[2]!=0)){ + // Add one guard cell in the case of the galilean algorithm + 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. -- cgit v1.2.3