Fix stair-case solver in 2D Cartesian (#2942)

* Fix stair-case solver in 2D Cartesian

* Generalization for the macroscropic solver

* Update automated test

* [pre-commit.ci] auto fixes from pre-commit.com hooks

for more information, see https://pre-commit.ci

* Update benchmark

Co-authored-by: pre-commit-ci[bot] <66853113+pre-commit-ci[bot]@users.noreply.github.com>

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)