From c78b9ab3a6cbf390e0ca32e80e00ff1d589cc0a4 Mon Sep 17 00:00:00 2001 From: David Grote Date: Mon, 12 Jul 2021 12:17:11 -0700 Subject: Fixes to RZ PSATD (#1945) * For RZ, changed the sign of the density corrections near the axis * Further fixes for deposition correction near axis * Yet one more sign fix for charge density * For RZ spectral solver, filled in the guard cells below the radial axis * Fix white space at end of line * In RZ spectral backtransform, ensure box is valid * For RZ inverse volume scaling, fixed use of nGrow to use nGrowVect * Temporary fix adding damped cells in the domain interior * Bug fix for RZ PSATD scalar backward transform * Fixes for damping of the fields in the z-guards * Bug fix in DampFieldsInGuards * Bug fix in DampFieldsInGuards (for tiling) * Added warpx_amr_check_input input parameter * Removed unneeded damp and zero_in_domain input * Removed damping related code from picmi * Improved some comments in code copying field to the radial guard cells * Update Source/FieldSolver/SpectralSolver/SpectralFieldDataRZ.cpp Simplify the expression for the sign Co-authored-by: Edoardo Zoni <59625522+EZoni@users.noreply.github.com> * Updated benchmarks * Updated tolerance for Langmuir analysis script * Updated CI test galilean_rz_psatd_current_correction Co-authored-by: Edoardo Zoni <59625522+EZoni@users.noreply.github.com> --- Python/pywarpx/picmi.py | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) (limited to 'Python/pywarpx/picmi.py') diff --git a/Python/pywarpx/picmi.py b/Python/pywarpx/picmi.py index 257cadb9a..d73e954ef 100644 --- a/Python/pywarpx/picmi.py +++ b/Python/pywarpx/picmi.py @@ -533,7 +533,6 @@ class ElectromagneticSolver(picmistandard.PICMI_ElectromagneticSolver): self.psatd_current_correction = kw.pop('warpx_current_correction', None) self.psatd_update_with_rho = kw.pop('warpx_psatd_update_with_rho', None) self.psatd_do_time_averaging = kw.pop('warpx_psatd_do_time_averaging', None) - self.psatd_use_damp_fields_in_z_guard = kw.pop('warpx_use_damp_fields_in_z_guard', None) def initialize_inputs(self): @@ -549,7 +548,6 @@ class ElectromagneticSolver(picmistandard.PICMI_ElectromagneticSolver): pywarpx.psatd.current_correction = self.psatd_current_correction pywarpx.psatd.update_with_rho = self.psatd_update_with_rho pywarpx.psatd.do_time_averaging = self.psatd_do_time_averaging - pywarpx.psatd.use_damp_fields_in_z_guard = self.psatd_use_damp_fields_in_z_guard if self.grid.guard_cells is not None: pywarpx.psatd.nx_guard = self.grid.guard_cells[0] @@ -735,6 +733,7 @@ class Simulation(picmistandard.PICMI_Simulation): self.costs_heuristic_particles_wt = kw.pop('warpx_costs_heuristic_particles_wt', None) self.costs_heuristic_cells_wt = kw.pop('warpx_costs_heuristic_cells_wt', None) self.use_fdtd_nci_corr = kw.pop('warpx_use_fdtd_nci_corr', None) + self.amr_check_input = kw.pop('warpx_amr_check_input', None) self.inputs_initialized = False self.warpx_initialized = False @@ -772,6 +771,8 @@ class Simulation(picmistandard.PICMI_Simulation): pywarpx.particles.use_fdtd_nci_corr = self.use_fdtd_nci_corr + pywarpx.amr.check_input = self.amr_check_input + particle_shape = self.particle_shape for s in self.species: if s.particle_shape is not None: -- cgit v1.2.3