aboutsummaryrefslogtreecommitdiff
path: root/Examples/Modules
diff options
context:
space:
mode:
authorGravatar Luca Fedeli <luca.fedeli.88@gmail.com> 2019-11-15 02:09:49 +0100
committerGravatar Luca Fedeli <luca.fedeli.88@gmail.com> 2019-11-15 02:09:49 +0100
commitd19a19a9f9cd0c2e37cd76a27276b6243aea5b13 (patch)
tree916595639944499c238714d068528b1c13f6a9e5 /Examples/Modules
parentf02985414079dd492446a70c24c36e92cfac3889 (diff)
downloadWarpX-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.py36
-rw-r--r--Examples/Modules/qed/quantum_synchrotron/inputs.3d_test_optical_depth_evolution114
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
-####