diff options
Diffstat (limited to 'Examples')
41 files changed, 108 insertions, 148 deletions
diff --git a/Examples/Modules/RigidInjection/inputs b/Examples/Modules/RigidInjection/inputs index d7e3b77e8..3c0d52ed0 100644 --- a/Examples/Modules/RigidInjection/inputs +++ b/Examples/Modules/RigidInjection/inputs @@ -23,10 +23,6 @@ geometry.prob_hi = 2. 2. 4. warpx.verbose = 1 # Algorithms -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 # interpolation interpolation.nox = 3 diff --git a/Examples/Modules/boosted_diags/inputs.2d b/Examples/Modules/boosted_diags/inputs.2d index 33eee5392..528eb6cd9 100644 --- a/Examples/Modules/boosted_diags/inputs.2d +++ b/Examples/Modules/boosted_diags/inputs.2d @@ -24,10 +24,7 @@ geometry.prob_hi = 150.e-6 150.e-6 0. warpx.verbose = 1 # Algorithms -algo.current_deposition = 3 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 +algo.current_deposition = direct # Numerics interpolation.nox = 3 diff --git a/Examples/Modules/boosted_diags/inputs.3d b/Examples/Modules/boosted_diags/inputs.3d index 33eee5392..528eb6cd9 100644 --- a/Examples/Modules/boosted_diags/inputs.3d +++ b/Examples/Modules/boosted_diags/inputs.3d @@ -24,10 +24,7 @@ geometry.prob_hi = 150.e-6 150.e-6 0. warpx.verbose = 1 # Algorithms -algo.current_deposition = 3 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 +algo.current_deposition = direct # Numerics interpolation.nox = 3 diff --git a/Examples/Modules/charged_beam/inputs b/Examples/Modules/charged_beam/inputs index 36b991c3e..18b645281 100644 --- a/Examples/Modules/charged_beam/inputs +++ b/Examples/Modules/charged_beam/inputs @@ -23,10 +23,7 @@ geometry.prob_hi = 20.e-6 20.e-6 20.e-6 warpx.verbose = 1 # Algorithms -algo.current_deposition = 3 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 +algo.current_deposition = direct # CFL warpx.cfl = 1.0 diff --git a/Examples/Modules/gaussian_beam/gaussian_beam_PICMI.py b/Examples/Modules/gaussian_beam/gaussian_beam_PICMI.py index b9bfc2864..a3a2330d7 100644 --- a/Examples/Modules/gaussian_beam/gaussian_beam_PICMI.py +++ b/Examples/Modules/gaussian_beam/gaussian_beam_PICMI.py @@ -17,7 +17,7 @@ number_sim_particles = 32768 total_charge = 8.010883097437485e-07 beam_rms_size = 0.25 -electron_beam_divergence = -0.04 +electron_beam_divergence = -0.04*picmi.c em_order = 3 diff --git a/Examples/Modules/gaussian_beam/inputs b/Examples/Modules/gaussian_beam/inputs index 0cab8d140..46cd785f2 100644 --- a/Examples/Modules/gaussian_beam/inputs +++ b/Examples/Modules/gaussian_beam/inputs @@ -23,10 +23,6 @@ geometry.prob_hi = 2. 2. 2. warpx.verbose = 1 # Algorithms -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 # interpolation interpolation.nox = 3 diff --git a/Examples/Modules/laser_injection/inputs b/Examples/Modules/laser_injection/inputs index 4186f2682..91ea135d1 100644 --- a/Examples/Modules/laser_injection/inputs +++ b/Examples/Modules/laser_injection/inputs @@ -21,10 +21,6 @@ geometry.prob_hi = 20.e-6 20.e-6 12.e-6 warpx.verbose = 1 # Algorithms -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 # CFL warpx.cfl = 1.0 diff --git a/Examples/Modules/laser_injection/inputs.2d.rt b/Examples/Modules/laser_injection/inputs.2d.rt index 1c40525f8..db726a98c 100644 --- a/Examples/Modules/laser_injection/inputs.2d.rt +++ b/Examples/Modules/laser_injection/inputs.2d.rt @@ -23,10 +23,7 @@ warpx.serialize_ics = 1 warpx.verbose = 1 # Algorithms -algo.current_deposition = 1 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 +algo.current_deposition = esirkepov # CFL warpx.cfl = 1.0 diff --git a/Examples/Modules/laser_injection/inputs.rt b/Examples/Modules/laser_injection/inputs.rt index b6f49d661..cec710403 100644 --- a/Examples/Modules/laser_injection/inputs.rt +++ b/Examples/Modules/laser_injection/inputs.rt @@ -23,10 +23,6 @@ warpx.serialize_ics = 1 warpx.verbose = 1 # Algorithms -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 # CFL warpx.cfl = 1.0 diff --git a/Examples/Modules/nci_corrector/inputs2d b/Examples/Modules/nci_corrector/inputs2d index 2f022f4cf..8db312fc6 100644 --- a/Examples/Modules/nci_corrector/inputs2d +++ b/Examples/Modules/nci_corrector/inputs2d @@ -25,10 +25,7 @@ warpx.verbose = 1 warpx.use_filter = 1 # Algorithms -algo.current_deposition = 1 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 +algo.current_deposition = esirkepov interpolation.nox = 3 interpolation.noy = 3 interpolation.noz = 3 diff --git a/Examples/Physics_applications/laser_acceleration/inputs.2d b/Examples/Physics_applications/laser_acceleration/inputs.2d index 11ddc6b0b..1a76feacf 100644 --- a/Examples/Physics_applications/laser_acceleration/inputs.2d +++ b/Examples/Physics_applications/laser_acceleration/inputs.2d @@ -17,10 +17,6 @@ warpx.fine_tag_hi = 5.e-6 -25.e-6 ################################# ############ NUMERICS ########### ################################# -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 interpolation.nox = 3 # Particle interpolation order. Must be the same in x, y, and z interpolation.noy = 3 interpolation.noz = 3 diff --git a/Examples/Physics_applications/laser_acceleration/inputs.2d.boost b/Examples/Physics_applications/laser_acceleration/inputs.2d.boost index 18c354634..fb03ae3f6 100644 --- a/Examples/Physics_applications/laser_acceleration/inputs.2d.boost +++ b/Examples/Physics_applications/laser_acceleration/inputs.2d.boost @@ -1,8 +1,8 @@ ################################# ######### BOX PARAMETERS ######## ################################# -# max_step = 2700 -stop_time = 1.9e-12 +max_step = 1000 +# stop_time = 1.9e-12 amr.n_cell = 128 1024 amr.max_grid_size = 64 amr.blocking_factor = 32 @@ -19,10 +19,8 @@ geometry.prob_hi = 128.e-6 0.96e-6 ################################# warpx.verbose = 1 amrex.v = 1 -algo.current_deposition = 2 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 1 +algo.current_deposition = direct-vectorized +algo.particle_pusher = vay algo.maxwell_fdtd_solver = ckc interpolation.nox = 3 interpolation.noy = 3 @@ -39,7 +37,7 @@ warpx.serialize_ics = 1 ################################# ####### BOOST PARAMETERS ######## ################################# -warpx.gamma_boost = 10. +warpx.gamma_boost = 30. warpx.boost_direction = z warpx.do_boosted_frame_diagnostic = 1 warpx.num_snapshots_lab = 7 @@ -62,11 +60,11 @@ electrons.momentum_distribution_type = "gaussian" electrons.xmin = -120.e-6 electrons.xmax = 120.e-6 electrons.zmin = 0.5e-3 -electrons.zmax = .0035 +electrons.zmax = 1. electrons.profile = "predefined" electrons.predefined_profile_name = "parabolic_channel" # predefined_profile_params = z_start ramp_up plateau ramp_down rc n0 -electrons.predefined_profile_params = .5e-3 .5e-3 2.e-3 .5e-3 50.e-6 3.5e24 +electrons.predefined_profile_params = .5e-3 .5e-3 2.e-3 .5e-3 50.e-6 3.5e25 electrons.do_continuous_injection = 1 ions.charge = q_e @@ -77,11 +75,11 @@ ions.momentum_distribution_type = "gaussian" ions.xmin = -120.e-6 ions.xmax = 120.e-6 ions.zmin = 0.5e-3 -ions.zmax = .0035 +ions.zmax = 1. ions.profile = "predefined" ions.predefined_profile_name = "parabolic_channel" # predefined_profile_params = z_start ramp_up plateau ramp_down rc n0 -ions.predefined_profile_params = .5e-3 .5e-3 2.e-3 .5e-3 50.e-6 3.5e24 +ions.predefined_profile_params = .5e-3 .5e-3 2.e-3 .5e-3 50.e-6 3.5e25 ions.do_continuous_injection = 1 beam.charge = -q_e diff --git a/Examples/Physics_applications/laser_acceleration/inputs.3d b/Examples/Physics_applications/laser_acceleration/inputs.3d index 8e854cbf4..f055708d1 100644 --- a/Examples/Physics_applications/laser_acceleration/inputs.3d +++ b/Examples/Physics_applications/laser_acceleration/inputs.3d @@ -17,10 +17,6 @@ warpx.fine_tag_hi = 5.e-6 5.e-6 -30.e-6 ################################# ############ NUMERICS ########### ################################# -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 interpolation.nox = 3 # Particle interpolation order. Must be the same in x, y, and z interpolation.noy = 3 interpolation.noz = 3 diff --git a/Examples/Physics_applications/plasma_acceleration/inputs.2d b/Examples/Physics_applications/plasma_acceleration/inputs.2d index e092cf932..9d0edbc7e 100644 --- a/Examples/Physics_applications/plasma_acceleration/inputs.2d +++ b/Examples/Physics_applications/plasma_acceleration/inputs.2d @@ -18,10 +18,6 @@ warpx.fine_tag_hi = 12.e-6 -100.e-6 ################################# ############ NUMERICS ########### ################################# -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 algo.maxwell_fdtd_solver = ckc warpx.use_filter = 1 warpx.do_pml = 1 diff --git a/Examples/Physics_applications/plasma_acceleration/inputs.2d.boost b/Examples/Physics_applications/plasma_acceleration/inputs.2d.boost index 59e7c432b..ba9166dee 100644 --- a/Examples/Physics_applications/plasma_acceleration/inputs.2d.boost +++ b/Examples/Physics_applications/plasma_acceleration/inputs.2d.boost @@ -16,10 +16,6 @@ geometry.prob_hi = 125.e-6 1.e-6 ################################# ############ NUMERICS ########### ################################# -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 algo.maxwell_fdtd_solver = ckc warpx.verbose = 1 warpx.plot_raw_fields = 0 diff --git a/Examples/Physics_applications/plasma_acceleration/inputs.3d.boost b/Examples/Physics_applications/plasma_acceleration/inputs.3d.boost index 1d58ef940..1af4fac86 100644 --- a/Examples/Physics_applications/plasma_acceleration/inputs.3d.boost +++ b/Examples/Physics_applications/plasma_acceleration/inputs.3d.boost @@ -16,10 +16,6 @@ geometry.prob_hi = 0.00015 0.00015 1.e-06 ################################# ############ NUMERICS ########### ################################# -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 algo.maxwell_fdtd_solver = ckc warpx.verbose = 1 warpx.plot_raw_fields = 0 diff --git a/Examples/Physics_applications/plasma_mirror/inputs.2d b/Examples/Physics_applications/plasma_mirror/inputs.2d index e16313b42..dc0a48ebf 100644 --- a/Examples/Physics_applications/plasma_mirror/inputs.2d +++ b/Examples/Physics_applications/plasma_mirror/inputs.2d @@ -18,10 +18,6 @@ warpx.serialize_ics = 1 ################################# ############ NUMERICS ########### ################################# -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 interpolation.nox = 3 interpolation.noy = 3 interpolation.noz = 3 diff --git a/Examples/Physics_applications/uniform_plasma/inputs.2d b/Examples/Physics_applications/uniform_plasma/inputs.2d index f3eb2decb..4304638dd 100644 --- a/Examples/Physics_applications/uniform_plasma/inputs.2d +++ b/Examples/Physics_applications/uniform_plasma/inputs.2d @@ -16,10 +16,6 @@ geometry.prob_hi = 20.e-6 20.e-6 ################################# warpx.serialize_ics = 1 warpx.verbose = 1 -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 warpx.cfl = 1.0 amr.plot_int = 10 diff --git a/Examples/Physics_applications/uniform_plasma/inputs.3d b/Examples/Physics_applications/uniform_plasma/inputs.3d index 245f3204c..6792dac2c 100644 --- a/Examples/Physics_applications/uniform_plasma/inputs.3d +++ b/Examples/Physics_applications/uniform_plasma/inputs.3d @@ -16,10 +16,6 @@ geometry.prob_hi = 20.e-6 20.e-6 20.e-6 ################################# warpx.serialize_ics = 1 warpx.verbose = 1 -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 warpx.cfl = 1.0 amr.plot_int = 10 diff --git a/Examples/Tests/Langmuir/inputs b/Examples/Tests/Langmuir/inputs index cd13aebbd..a989f2249 100644 --- a/Examples/Tests/Langmuir/inputs +++ b/Examples/Tests/Langmuir/inputs @@ -26,10 +26,7 @@ warpx.moving_window_dir = x warpx.moving_window_v = 0.0 # in units of the speed of light # Algorithms -algo.current_deposition = 3 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 +algo.current_deposition = direct # CFL warpx.cfl = 1.0 diff --git a/Examples/Tests/Langmuir/inputs.2d.rz b/Examples/Tests/Langmuir/inputs.2d.rz index c4132d493..a4464fe1d 100644 --- a/Examples/Tests/Langmuir/inputs.2d.rz +++ b/Examples/Tests/Langmuir/inputs.2d.rz @@ -27,10 +27,7 @@ warpx.moving_window_dir = z warpx.moving_window_v = 0.0 # in units of the speed of light # Algorithms -algo.current_deposition = 3 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 +algo.current_deposition = direct interpolation.nox = 1 interpolation.noy = 1 interpolation.noz = 1 diff --git a/Examples/Tests/Langmuir/inputs.lb b/Examples/Tests/Langmuir/inputs.lb index 5900cb1bb..73ededb23 100644 --- a/Examples/Tests/Langmuir/inputs.lb +++ b/Examples/Tests/Langmuir/inputs.lb @@ -32,10 +32,7 @@ warpx.moving_window_v = 0.0 # in units of the speed of light warpx.regrid_int = 1 # Algorithms -algo.current_deposition = 3 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 +algo.current_deposition = direct # CFL warpx.cfl = 1.0 diff --git a/Examples/Tests/Langmuir/inputs.multi.2d.rt b/Examples/Tests/Langmuir/inputs.multi.2d.rt index 79bc2b383..54d4e2c16 100644 --- a/Examples/Tests/Langmuir/inputs.multi.2d.rt +++ b/Examples/Tests/Langmuir/inputs.multi.2d.rt @@ -25,10 +25,7 @@ warpx.serialize_ics = 1 warpx.verbose = 1 # Algorithms -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 1 -algo.particle_pusher = 0 +algo.field_gathering = standard # Interpolation interpolation.nox = 1 diff --git a/Examples/Tests/Langmuir/inputs.multi.rt b/Examples/Tests/Langmuir/inputs.multi.rt index c9969d39a..8e6d6d8e1 100644 --- a/Examples/Tests/Langmuir/inputs.multi.rt +++ b/Examples/Tests/Langmuir/inputs.multi.rt @@ -26,10 +26,8 @@ warpx.serialize_ics = 1 warpx.verbose = 1 # Algorithms -algo.current_deposition = 3 -algo.charge_deposition = 0 -algo.field_gathering = 1 -algo.particle_pusher = 0 +algo.current_deposition = direct +algo.field_gathering = standard # Interpolation interpolation.nox = 1 diff --git a/Examples/Tests/Langmuir/inputs.multi.rz.rt b/Examples/Tests/Langmuir/inputs.multi.rz.rt index 6bc84a434..e4099f9c5 100644 --- a/Examples/Tests/Langmuir/inputs.multi.rz.rt +++ b/Examples/Tests/Langmuir/inputs.multi.rz.rt @@ -25,10 +25,7 @@ warpx.serialize_ics = 1 warpx.verbose = 1 # Algorithms -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 1 -algo.particle_pusher = 0 +algo.field_gathering = standard # Interpolation interpolation.nox = 1 diff --git a/Examples/Tests/Langmuir/inputs.nolb b/Examples/Tests/Langmuir/inputs.nolb index cd2900b29..6ac0629cb 100644 --- a/Examples/Tests/Langmuir/inputs.nolb +++ b/Examples/Tests/Langmuir/inputs.nolb @@ -32,10 +32,7 @@ warpx.moving_window_v = 0.0 # in units of the speed of light warpx.regrid_int = -1 # Algorithms -algo.current_deposition = 3 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 +algo.current_deposition = direct # CFL warpx.cfl = 1.0 diff --git a/Examples/Tests/Langmuir/inputs.rt b/Examples/Tests/Langmuir/inputs.rt index 3ec8fc47d..be5a730c2 100644 --- a/Examples/Tests/Langmuir/inputs.rt +++ b/Examples/Tests/Langmuir/inputs.rt @@ -25,10 +25,7 @@ warpx.serialize_ics = 1 warpx.verbose = 1 # Algorithms -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 1 -algo.particle_pusher = 0 +algo.field_gathering = standard # Interpolation interpolation.nox = 1 diff --git a/Examples/Tests/Langmuir/langmuir2d_PICMI.py b/Examples/Tests/Langmuir/langmuir2d_PICMI.py index c8c5a9559..8108a39ba 100644 --- a/Examples/Tests/Langmuir/langmuir2d_PICMI.py +++ b/Examples/Tests/Langmuir/langmuir2d_PICMI.py @@ -9,7 +9,9 @@ ymin = -20.e-6 xmax = +20.e-6 ymax = +20.e-6 -uniform_plasma = picmi.UniformDistribution(density=1.e25, upper_bound=[0., None, None], directed_velocity=[0.1, 0., 0.]) +uniform_plasma = picmi.UniformDistribution(density = 1.e25, + upper_bound = [0., None, None], + directed_velocity = [0.1*picmi.c, 0., 0.]) electrons = picmi.Species(particle_type='electron', name='electrons', initial_distribution=uniform_plasma) diff --git a/Examples/Tests/Langmuir/langmuir_PICMI.py b/Examples/Tests/Langmuir/langmuir_PICMI.py index 67f37ec65..d09c90714 100644 --- a/Examples/Tests/Langmuir/langmuir_PICMI.py +++ b/Examples/Tests/Langmuir/langmuir_PICMI.py @@ -14,7 +14,9 @@ xmax = +20.e-6 ymax = +20.e-6 zmax = +20.e-6 -uniform_plasma = picmi.UniformDistribution(density=1.e25, upper_bound=[0., None, None], directed_velocity=[0.1, 0., 0.]) +uniform_plasma = picmi.UniformDistribution(density = 1.e25, + upper_bound = [0., None, None], + directed_velocity = [0.1*picmi.c, 0., 0.]) electrons = picmi.Species(particle_type='electron', name='electrons', initial_distribution=uniform_plasma) diff --git a/Examples/Tests/Langmuir/langmuir_PICMI_rt.py b/Examples/Tests/Langmuir/langmuir_PICMI_rt.py index 30540313e..e629e268d 100644 --- a/Examples/Tests/Langmuir/langmuir_PICMI_rt.py +++ b/Examples/Tests/Langmuir/langmuir_PICMI_rt.py @@ -14,7 +14,9 @@ xmax = +20.e-6 ymax = +20.e-6 zmax = +20.e-6 -uniform_plasma = picmi.UniformDistribution(density=1.e25, upper_bound=[0., None, None], directed_velocity=[0.1, 0., 0.]) +uniform_plasma = picmi.UniformDistribution(density = 1.e25, + upper_bound = [0., None, None], + directed_velocity = [0.1*picmi.c, 0., 0.]) electrons = picmi.Species(particle_type='electron', name='electrons', initial_distribution=uniform_plasma) diff --git a/Examples/Tests/Langmuir/langmuir_PICMI_rz.py b/Examples/Tests/Langmuir/langmuir_PICMI_rz.py index 4eb998177..88b3947a6 100644 --- a/Examples/Tests/Langmuir/langmuir_PICMI_rz.py +++ b/Examples/Tests/Langmuir/langmuir_PICMI_rz.py @@ -9,7 +9,9 @@ zmin = -20.e-6 rmax = +20.e-6 zmax = +20.e-6 -uniform_plasma = picmi.UniformDistribution(density=1.e25, upper_bound=[None, None, 0.], directed_velocity=[0., 0., 0.1]) +uniform_plasma = picmi.UniformDistribution(density = 1.e25, + upper_bound = [None, None, 0.], + directed_velocity = [0., 0., 0.1*picmi.c]) electrons = picmi.Species(particle_type='electron', name='electrons', initial_distribution=uniform_plasma) diff --git a/Examples/Tests/Larmor/inputs b/Examples/Tests/Larmor/inputs index 729d315e9..da8fb5bae 100644 --- a/Examples/Tests/Larmor/inputs +++ b/Examples/Tests/Larmor/inputs @@ -41,10 +41,6 @@ warpx.B_external = 0.0 0.00078110417851950768 0.0 warpx.verbose = 1 # Algorithms -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 # CFL warpx.cfl = 1.0 diff --git a/Examples/Tests/Larmor/inputs.ml b/Examples/Tests/Larmor/inputs.ml index 9a0c9eb3d..3cfb978d7 100644 --- a/Examples/Tests/Larmor/inputs.ml +++ b/Examples/Tests/Larmor/inputs.ml @@ -41,10 +41,6 @@ warpx.B_external = 0.0 0.00078110417851950768 0.0 warpx.verbose = 1 # Algorithms -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 # CFL warpx.cfl = 1.0 diff --git a/Examples/Tests/PML/inputs2d b/Examples/Tests/PML/inputs2d index c84cf0555..5b936a333 100644 --- a/Examples/Tests/PML/inputs2d +++ b/Examples/Tests/PML/inputs2d @@ -23,10 +23,6 @@ geometry.prob_hi = 30.e-6 120.e-6 warpx.verbose = 1 # Algorithms -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 warpx.cfl = 1.0 warpx.do_pml = 1 diff --git a/Examples/Tests/SingleParticle/bilinear_filter_analysis.py b/Examples/Tests/SingleParticle/bilinear_filter_analysis.py new file mode 100755 index 000000000..1a47796fb --- /dev/null +++ b/Examples/Tests/SingleParticle/bilinear_filter_analysis.py @@ -0,0 +1,38 @@ +#! /usr/bin/env python + +import sys +import yt ; yt.funcs.mylog.setLevel(0) +import numpy as np +from scipy import signal + +# Build Jx without filter (from other simulation) +my_F_nofilter = np.zeros([16,16]) +my_F_nofilter[8,8] = -1.601068065642412e-11 +my_F_nofilter[8,7] = -1.601068065642412e-11 + +# Build 2D filter +filter0 = np.array([.25,.5,.25]) +my_order = [1,5] +my_filterx = filter0 +my_filtery = filter0 +while my_order[0]>1: + my_filterx = np.convolve(my_filterx,filter0) + my_order[0] -= 1 +while my_order[1]>1: + my_filtery = np.convolve(my_filtery,filter0) + my_order[1] -= 1 +my_filter = my_filterx[:,None]*my_filtery + +# Apply filter. my_F_filetered is the theoretical value for filtered field +my_F_filtered = signal.convolve2d(my_F_nofilter, my_filter, boundary='symm', mode='same') + +# Get simulation result for F_filtered +filename = sys.argv[1] +ds = yt.load( filename ) +sl = yt.SlicePlot(ds, 2, 'jx', aspect=1) +all_data_level_0 = ds.covering_grid(level=0,left_edge=ds.domain_left_edge, dims=ds.domain_dimensions) +F_filtered = all_data_level_0['boxlib', 'jx'].v.squeeze() + +# Compare theory and PIC for filtered value +error = np.sum( np.abs(F_filtered - my_F_filtered) ) / np.sum( np.abs(my_F_filtered) ) +assert( error < 1.e-14 ) diff --git a/Examples/Tests/SingleParticle/inputs b/Examples/Tests/SingleParticle/inputs new file mode 100644 index 000000000..548848d79 --- /dev/null +++ b/Examples/Tests/SingleParticle/inputs @@ -0,0 +1,23 @@ +max_step = 1 +amr.n_cell = 16 16 +amr.max_level = 0 +amr.blocking_factor = 8 +amr.max_grid_size = 8 +amr.plot_int = 1 +geometry.coord_sys = 0 +geometry.is_periodic = 0 0 +geometry.prob_lo = -8 -12 +geometry.prob_hi = 8 12 +warpx.do_pml = 0 +algo.charge_deposition = standard +algo.field_gathering = standard +warpx.cfl = 1.0 + +particles.nspecies = 1 +particles.species_names = electron +electron.charge = -q_e +electron.mass = m_e +electron.injection_style = "SingleParticle" +electron.single_particle_pos = 0.0 0.0 0.0 +electron.single_particle_vel = 1.e20 0.0 0.0 # gamma*beta +electron.single_particle_weight = 1.0 diff --git a/Examples/Tests/gpu_test/inputs b/Examples/Tests/gpu_test/inputs index 2b719655d..e4ae27469 100644 --- a/Examples/Tests/gpu_test/inputs +++ b/Examples/Tests/gpu_test/inputs @@ -24,10 +24,8 @@ warpx.do_pml = 0 warpx.verbose = 1 # Algorithms -algo.current_deposition = 3 -algo.charge_deposition = 0 -algo.field_gathering = 1 -algo.particle_pusher = 0 +algo.current_deposition = direct +algo.field_gathering = standard interpolation.nox = 1 interpolation.noy = 1 diff --git a/Examples/Tests/laser_on_fine/inputs b/Examples/Tests/laser_on_fine/inputs index 8fc8876cb..9d3133dd7 100644 --- a/Examples/Tests/laser_on_fine/inputs +++ b/Examples/Tests/laser_on_fine/inputs @@ -38,10 +38,9 @@ warpx.pml_ncell = 10 warpx.verbose = 1 # Algorithms -algo.current_deposition = 1 -algo.charge_deposition = 1 -algo.field_gathering = 1 -algo.particle_pusher = 0 +algo.current_deposition = esirkepov +algo.charge_deposition = standard +algo.field_gathering = standard # CFL warpx.cfl = 1.0 diff --git a/Examples/Tests/laser_on_fine/inputs.2d b/Examples/Tests/laser_on_fine/inputs.2d index c4af2df56..d011ff4ee 100644 --- a/Examples/Tests/laser_on_fine/inputs.2d +++ b/Examples/Tests/laser_on_fine/inputs.2d @@ -38,10 +38,9 @@ warpx.pml_ncell = 10 warpx.verbose = 1 # Algorithms -algo.current_deposition = 1 -algo.charge_deposition = 1 -algo.field_gathering = 1 -algo.particle_pusher = 0 +algo.current_deposition = esirkepov +algo.charge_deposition = standard +algo.field_gathering = standard # CFL warpx.cfl = 1.0 diff --git a/Examples/Tests/self_force_test/inputs b/Examples/Tests/self_force_test/inputs index 1e55c1174..5e7bcc4f8 100644 --- a/Examples/Tests/self_force_test/inputs +++ b/Examples/Tests/self_force_test/inputs @@ -42,10 +42,9 @@ warpx.pml_ncell = 10 warpx.verbose = 1 # Algorithms -algo.current_deposition = 1 -algo.charge_deposition = 1 -algo.field_gathering = 1 -algo.particle_pusher = 0 +algo.current_deposition = esirkepov +algo.charge_deposition = standard +algo.field_gathering = standard # CFL warpx.cfl = 1.0 diff --git a/Examples/Tests/subcycling/inputs.2d b/Examples/Tests/subcycling/inputs.2d index 9790ee973..15034be52 100644 --- a/Examples/Tests/subcycling/inputs.2d +++ b/Examples/Tests/subcycling/inputs.2d @@ -37,10 +37,6 @@ warpx.n_current_deposition_buffer = 0 warpx.n_field_gather_buffer = 0 # Algorithms -algo.current_deposition = 0 -algo.charge_deposition = 0 -algo.field_gathering = 0 -algo.particle_pusher = 0 algo.maxwell_fdtd_solver = "ckc" # CFL |