aboutsummaryrefslogtreecommitdiff
path: root/Examples/Modules/embedded_boundary_cube/analysis_fields_2d.py
diff options
context:
space:
mode:
Diffstat (limited to 'Examples/Modules/embedded_boundary_cube/analysis_fields_2d.py')
-rwxr-xr-xExamples/Modules/embedded_boundary_cube/analysis_fields_2d.py5
1 files changed, 4 insertions, 1 deletions
diff --git a/Examples/Modules/embedded_boundary_cube/analysis_fields_2d.py b/Examples/Modules/embedded_boundary_cube/analysis_fields_2d.py
index 227048be9..67e893b19 100755
--- a/Examples/Modules/embedded_boundary_cube/analysis_fields_2d.py
+++ b/Examples/Modules/embedded_boundary_cube/analysis_fields_2d.py
@@ -48,7 +48,6 @@ for i in range(ncells[0]):
(-Lx / 2 <= x < Lx / 2) *
(-Lz / 2 <= z < Lz / 2) *
np.cos(np.pi / Lx * c * t))
-
rel_tol_err = 1e-3
# Compute relative l^2 error on By
@@ -56,6 +55,10 @@ By_sim = data['By'].to_ndarray()
rel_err_y = np.sqrt(np.sum(np.square(By_sim - By_th)) / np.sum(np.square(By_th)))
assert (rel_err_y < rel_tol_err)
+# Compute relative l^2 error on Ey
+Ey_sim = data['Ey'].to_ndarray()
+rel_err_y = np.sqrt(np.sum(np.square(Ey_sim/c - By_th)) / np.sum(np.square(By_th)))
+
test_name = os.path.split(os.getcwd())[1]
checksumAPI.evaluate_checksum(test_name, filename)