diff options
Diffstat (limited to 'Examples/Tests/PML/analysis_pml_psatd.py')
| -rwxr-xr-x | Examples/Tests/PML/analysis_pml_psatd.py | 5 |
1 files changed, 5 insertions, 0 deletions
diff --git a/Examples/Tests/PML/analysis_pml_psatd.py b/Examples/Tests/PML/analysis_pml_psatd.py index 164c7a19e..3e5bcef04 100755 --- a/Examples/Tests/PML/analysis_pml_psatd.py +++ b/Examples/Tests/PML/analysis_pml_psatd.py @@ -31,6 +31,8 @@ Bz = all_data_level_0['boxlib', 'Bz'].v.squeeze() Ex = all_data_level_0['boxlib', 'Ex'].v.squeeze() Ey = all_data_level_0['boxlib', 'Ey'].v.squeeze() Ez = all_data_level_0['boxlib', 'Ez'].v.squeeze() +rho = all_data_level_0['boxlib','rho'].v.squeeze() +divE = all_data_level_0['boxlib','divE'].v.squeeze() energyE = np.sum(scc.epsilon_0/2*(Ex**2+Ey**2+Ez**2)) energyB = np.sum(1./scc.mu_0/2*(Bx**2+By**2+Bz**2)) energy_end = energyE + energyB @@ -43,3 +45,6 @@ print("Reflectivity_theory: %s" %Reflectivity_theory) assert( abs(Reflectivity-Reflectivity_theory) < 5./100 * Reflectivity_theory ) +# Check relative L-infinity spatial norm of rho/epsilon_0 - div(E) +Linf_norm = np.amax( np.abs( rho/scc.epsilon_0 - divE ) ) / np.amax( np.abs( rho/scc.epsilon_0 ) ) +assert( Linf_norm < 2.e-2 ) |