diff options
Diffstat (limited to 'Examples')
| -rwxr-xr-x | Examples/Modules/embedded_boundary_cube/analysis_fields_2d.py | 5 | ||||
| -rw-r--r-- | Examples/Modules/embedded_boundary_cube/inputs_2d | 5 |
2 files changed, 9 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) diff --git a/Examples/Modules/embedded_boundary_cube/inputs_2d b/Examples/Modules/embedded_boundary_cube/inputs_2d index 2b3348582..476992b36 100644 --- a/Examples/Modules/embedded_boundary_cube/inputs_2d +++ b/Examples/Modules/embedded_boundary_cube/inputs_2d @@ -20,6 +20,7 @@ my_constants.zmax = 0.5 warpx.eb_implicit_function = "max(max(x+xmin,-(x+xmax)), max(z+zmin,-(z+zmax)))" warpx.B_ext_grid_init_style = parse_B_ext_grid_function +warpx.E_ext_grid_init_style = parse_E_ext_grid_function my_constants.m = 0 my_constants.p = 1 @@ -30,6 +31,10 @@ warpx.Bz_external_grid_function(x,y,z) = 0 warpx.Bx_external_grid_function(x,y,z) = 0 warpx.By_external_grid_function(x,y,z) = cos(m * pi / Lx * (x - Lx / 2)) * cos(p * pi / Lz * (z - Lz / 2))*mu0 +warpx.Ez_external_grid_function(x,y,z) = 0 +warpx.Ex_external_grid_function(x,y,z) = 0 +warpx.Ey_external_grid_function(x,y,z) = cos(m * pi / Lx * (x - Lx / 2)) * cos(p * pi / Lz * (z - Lz / 2))*mu0*clight + diagnostics.diags_names = diag1 diag1.intervals = 1 diag1.diag_type = Full |