diff options
author | 2019-11-15 02:09:49 +0100 | |
---|---|---|
committer | 2019-11-15 02:09:49 +0100 | |
commit | d19a19a9f9cd0c2e37cd76a27276b6243aea5b13 (patch) | |
tree | 916595639944499c238714d068528b1c13f6a9e5 /Examples/Modules | |
parent | f02985414079dd492446a70c24c36e92cfac3889 (diff) | |
download | WarpX-d19a19a9f9cd0c2e37cd76a27276b6243aea5b13.tar.gz WarpX-d19a19a9f9cd0c2e37cd76a27276b6243aea5b13.tar.zst WarpX-d19a19a9f9cd0c2e37cd76a27276b6243aea5b13.zip |
removed unwanted test
Diffstat (limited to 'Examples/Modules')
-rw-r--r-- | Examples/Modules/qed/quantum_synchrotron/analysis_3d_optical_depth_evolution.py | 36 | ||||
-rw-r--r-- | Examples/Modules/qed/quantum_synchrotron/inputs.3d_test_optical_depth_evolution | 114 |
2 files changed, 0 insertions, 150 deletions
diff --git a/Examples/Modules/qed/quantum_synchrotron/analysis_3d_optical_depth_evolution.py b/Examples/Modules/qed/quantum_synchrotron/analysis_3d_optical_depth_evolution.py deleted file mode 100644 index d38a2d47e..000000000 --- a/Examples/Modules/qed/quantum_synchrotron/analysis_3d_optical_depth_evolution.py +++ /dev/null @@ -1,36 +0,0 @@ -#! /usr/bin/env python3 -import yt -import numpy as np -import scipy.stats as st -import sys - -# This script checks if electrons and positrons initialized with -# Quantum Synchrotron process enabled -# do actually have an exponentially distributed optical depth - -# Tolerance -tol = 1e-2 - -def check(): - return - filename = sys.argv[1] - data_set = yt.load(filename) - - all_data = data_set.all_data() - res_ele_tau = all_data["electrons", 'particle_tau'] - res_pos_tau = all_data["positrons", 'particle_tau'] - - loc_ele, scale_ele = st.expon.fit(res_ele_tau) - loc_pos, scale_pos = st.expon.fit(res_pos_tau) - - # loc should be very close to 0, scale should be very close to 1 - assert(np.abs(loc_ele - 0) < tol) - assert(np.abs(loc_pos - 0) < tol) - assert(np.abs(scale_ele - 1) < tol) - assert(np.abs(scale_pos - 1) < tol) - -def main(): - check() - -if __name__ == "__main__": - main() diff --git a/Examples/Modules/qed/quantum_synchrotron/inputs.3d_test_optical_depth_evolution b/Examples/Modules/qed/quantum_synchrotron/inputs.3d_test_optical_depth_evolution deleted file mode 100644 index 496838e4b..000000000 --- a/Examples/Modules/qed/quantum_synchrotron/inputs.3d_test_optical_depth_evolution +++ /dev/null @@ -1,114 +0,0 @@ -#An inputfile to demonstrate lookup table generation for the Breit Wheeler engine - -################################# -####### GENERAL PARAMETERS ###### -################################# -max_step = 10 -amr.n_cell = 64 64 64 -amr.max_grid_size = 32 # maximum size of each AMReX box, used to decompose the domain -amr.blocking_factor = 8 # minimum size of each AMReX box, used to decompose the domain -amr.plot_int = 10 -geometry.coord_sys = 0 # 0: Cartesian -geometry.is_periodic = 1 1 1 # Is periodic? -geometry.prob_lo = -1.e-6 -1.e-6 -1e-6 # physical domain -geometry.prob_hi = 1.e-6 1.e-6 1e-6 -amr.max_level = 0 # Maximum level in hierarchy (1 might be unstable, >1 is not supported) - - -################################# -############ NUMERICS ########### -################################# -algo.current_deposition = esirkepov -algo.charge_deposition = standard -algo.field_gathering = energy-conserving -algo.particle_pusher = boris -interpolation.nox = 3 # Particle interpolation order. Must be the same in x, y, and z -interpolation.noy = 3 -interpolation.noz = 3 -warpx.verbose = 1 -warpx.do_dive_cleaning = 0 -warpx.plot_raw_fields = 0 -warpx.plot_raw_fields_guards = 0 -warpx.plot_finepatch = 0 -warpx.plot_crsepatch = 0 -warpx.use_filter = 1 -warpx.cfl = 1. # if 1., the time step is set to its CFL limit -warpx.do_pml = 1 # use Perfectly Matched Layer as boundary condition -warpx.serialize_ics = 1 - -################################# -############ PLASMA ############# -################################# -particles.nspecies = 2 # number of species -particles.species_names = electrons positrons -################################# - -electrons.charge = -q_e -electrons.mass = m_e -electrons.injection_style = "NUniformPerCell" -electrons.profile = "constant" -electrons.xmin = -0.5e-6 -electrons.ymin = -0.5e-6 -electrons.zmin = -0.5e-6 -electrons.xmax = 0.5e-6 -electrons.ymax = 0.5e-6 -electrons.zmax = 0.5e-6 -electrons.num_particles_per_cell_each_dim = 2 2 -electrons.density = 1e19 -electrons.profile = "constant" -electrons.momentum_distribution_type = "gaussian" -electrons.ux_m = 0.0 -electrons.uy_m = 0.0 -electrons.uz_m = 0.0 -electrons.ux_th = 100. -electrons.uy_th = 100. -electrons.uz_th = 100. -##########QED#################### -electrons.do_qed = 1 -electrons.do_qed_quantum_sync = 1 -electrons.do_classical_radiation_reaction = 1 -################################# - -positrons.charge = q_e -positrons.mass = m_e -positrons.injection_style = "NUniformPerCell" -positrons.profile = "constant" -positrons.xmin = -0.5e-6 -positrons.ymin = -0.5e-6 -positrons.zmin = -0.5e-6 -positrons.xmax = 0.5e-6 -positrons.ymax = 0.5e-6 -positrons.zmax = 0.5e-6 -positrons.num_particles_per_cell_each_dim = 2 2 -positrons.density = 1e19 -positrons.profile = "constant" -positrons.momentum_distribution_type = "gaussian" -positrons.ux_m = 0.0 -positrons.uy_m = 0.0 -positrons.uz_m = 0.0 -positrons.ux_th = 100. -positrons.uy_th = 100. -positrons.uz_th = 100. -##########QED#################### -positrons.do_qed = 1 -positrons.do_qed_quantum_sync = 1 -positrons.do_classical_radiation_reaction = 1 -################################# - -##########QED TABLES#################### -qed_qs.chi_min = 0.001 -qed_qs.generate_table = 1 -qed_qs.tab_dndt_chi_min = 0.001 -qed_qs.tab_dndt_chi_max = 200 -qed_qs.tab_dndt_how_many = 64 -qed_qs.tab_em_chi_min = 0.001 -qed_qs.tab_em_chi_max = 200 -qed_qs.tab_em_chi_how_many = 2 -qed_qs.tab_em_prob_how_many = 2 -qed_qs.save_table_in = "qs_table" -#qed_qs.load_table_from = "qs_table" -################################# - -### EXTERNAL FIELD ### -warpx.B_external_particle = 0 0 1.0e7 -#### |