diff options
Diffstat (limited to 'Examples/Tests/particles_in_PML/analysis_particles_in_pml.py')
-rwxr-xr-x | Examples/Tests/particles_in_PML/analysis_particles_in_pml.py | 60 |
1 files changed, 0 insertions, 60 deletions
diff --git a/Examples/Tests/particles_in_PML/analysis_particles_in_pml.py b/Examples/Tests/particles_in_PML/analysis_particles_in_pml.py deleted file mode 100755 index 85878d3ad..000000000 --- a/Examples/Tests/particles_in_PML/analysis_particles_in_pml.py +++ /dev/null @@ -1,60 +0,0 @@ -#!/usr/bin/env python3 - -# Copyright 2019-2020 Luca Fedeli, Maxence Thevenet, Remi Lehe -# -# -# This file is part of WarpX. -# -# License: BSD-3-Clause-LBNL - -""" -This script tests the absorption of particles in the PML. - -The input file inputs_2d/inputs is used: it features a positive and a -negative particle, going in opposite direction and eventually -leaving the box. This script tests that the field in the box -is close to 0 once the particles have left. With regular -PML, this test fails, since the particles leave a spurious -charge, with associated fields, behind them. -""" -import os -import sys - -import yt - -yt.funcs.mylog.setLevel(0) -sys.path.insert(1, '../../../../warpx/Regression/Checksum/') -import checksumAPI - -# Open plotfile specified in command line -filename = sys.argv[1] -ds = yt.load( filename ) - -# Check that the field is low enough -ad0 = ds.covering_grid(level=0, left_edge=ds.domain_left_edge, dims=ds.domain_dimensions) -Ex_array = ad0[('mesh','Ex')].to_ndarray() -Ey_array = ad0[('mesh','Ey')].to_ndarray() -Ez_array = ad0[('mesh','Ez')].to_ndarray() -max_Efield = max(Ex_array.max(), Ey_array.max(), Ez_array.max()) -print( "max_Efield = %s" %max_Efield ) - -# The field associated with the particle does not have -# the same amplitude in 2d and 3d -if ds.dimensionality == 2: - if ds.max_level == 0: - tolerance_abs = 0.0003 - elif ds.max_level == 1: - tolerance_abs = 0.0006 -elif ds.dimensionality == 3: - if ds.max_level == 0: - tolerance_abs = 10 - elif ds.max_level == 1: - tolerance_abs = 110 -else: - raise ValueError("Unknown dimensionality") - -print("tolerance_abs: " + str(tolerance_abs)) -assert max_Efield < tolerance_abs - -test_name = os.path.split(os.getcwd())[1] -checksumAPI.evaluate_checksum(test_name, filename) |