From 81dd5a670799b72d0e99f04cb5028566ad767b18 Mon Sep 17 00:00:00 2001 From: Edoardo Zoni <59625522+EZoni@users.noreply.github.com> Date: Thu, 2 Jul 2020 16:02:38 -0700 Subject: PSATD: current correction works only with global FFTs (#1131) * Clean up implementation of current correction: - abort when current correction runs with psatd.periodic_single_box_fft=0 - rename input parameter from psatd.do_current_correction to psatd.current_correction - add nodal tests in 2D and 3D (using direct current deposition) - change 2D test: use Esirkepov deposition, instead of direct, as in 3D test - add/update relevant checksum benchmarks * Small clean-up --- Python/pywarpx/picmi.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'Python/pywarpx/picmi.py') diff --git a/Python/pywarpx/picmi.py b/Python/pywarpx/picmi.py index 6f4d7a4ba..df700a124 100644 --- a/Python/pywarpx/picmi.py +++ b/Python/pywarpx/picmi.py @@ -453,7 +453,7 @@ class ElectromagneticSolver(picmistandard.PICMI_ElectromagneticSolver): if self.method == 'PSATD': self.periodic_single_box_fft = kw.pop('warpx_periodic_single_box_fft', None) self.fftw_plan_measure = kw.pop('warpx_fftw_plan_measure', None) - self.do_current_correction = kw.pop('warpx_do_current_correction', None) + self.current_correction = kw.pop('warpx_current_correction', None) def initialize_inputs(self): @@ -466,7 +466,7 @@ class ElectromagneticSolver(picmistandard.PICMI_ElectromagneticSolver): if self.method == 'PSATD': pywarpx.psatd.periodic_single_box_fft = self.periodic_single_box_fft pywarpx.psatd.fftw_plan_measure = self.fftw_plan_measure - pywarpx.psatd.do_current_correction = self.do_current_correction + pywarpx.psatd.current_correction = self.current_correction if self.stencil_order is not None: pywarpx.psatd.nox = self.stencil_order[0] -- cgit v1.2.3