aboutsummaryrefslogtreecommitdiff
path: root/Examples/Tests/particles_in_PML/analysis_particles_in_pml.py
diff options
context:
space:
mode:
Diffstat (limited to 'Examples/Tests/particles_in_PML/analysis_particles_in_pml.py')
-rwxr-xr-xExamples/Tests/particles_in_PML/analysis_particles_in_pml.py60
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)