diff options
author | 2020-11-30 09:01:34 -0800 | |
---|---|---|
committer | 2020-11-30 09:01:34 -0800 | |
commit | 8f3ccb17b31d0b26fedebb3c095700f2014f315c (patch) | |
tree | f228ca33f506065252a95a9c3e5c8265425acf67 /Examples/Tests/galilean/analysis_2d.py | |
parent | bedec348b35ac28cf85cfe588e0ffd618fb2a1f9 (diff) | |
download | WarpX-8f3ccb17b31d0b26fedebb3c095700f2014f315c.tar.gz WarpX-8f3ccb17b31d0b26fedebb3c095700f2014f315c.tar.zst WarpX-8f3ccb17b31d0b26fedebb3c095700f2014f315c.zip |
Fix bug in yt workaround for Galilean 2D/RZ tests (#1536)
* Use negative velocity in 2D/RZ Galilean tests, remove yt workaround
* Remove unused import
Diffstat (limited to 'Examples/Tests/galilean/analysis_2d.py')
-rwxr-xr-x | Examples/Tests/galilean/analysis_2d.py | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/Examples/Tests/galilean/analysis_2d.py b/Examples/Tests/galilean/analysis_2d.py index ea423a71d..d6486a4e2 100755 --- a/Examples/Tests/galilean/analysis_2d.py +++ b/Examples/Tests/galilean/analysis_2d.py @@ -37,23 +37,23 @@ Ez= ds.index.grids[0]['boxlib', 'Ez'].squeeze().v if (averaged): # energyE_ref was calculated with Galilean PSATD method (v_galilean = (0,0,0.99498743710662)) - energyE_ref = 26913.546573259937 - tolerance_rel = 1e-5 + energyE_ref = 32532.00882239954 + tolerance_rel = 1e-6 elif (not dims_RZ and not current_correction): # energyE_ref was calculated with standard PSATD method (v_galilean = (0.,0.,0.)) - energyE_ref = 38362.88743899688 + energyE_ref = 35657.99361677053 tolerance_rel = 1e-8 elif (not dims_RZ and current_correction): # energyE_ref was calculated with standard PSATD method (v_galilean = (0.,0.,0.)): - energyE_ref = 33331.13857367177 + energyE_ref = 35024.02751955393 tolerance_rel = 2e-8 elif (dims_RZ and not current_correction): # energyE_ref was calculated with standard PSATD method (v_galilean = (0.,0.,0.)) - energyE_ref = 178013.54481470847 + energyE_ref = 239019.10670780553 tolerance_rel = 1e-8 elif (dims_RZ and current_correction): # energyE_ref was calculated with standard PSATD method (v_galilean = (0.,0.,0.)) - energyE_ref = 455119.0695414118 + energyE_ref = 471730.0524143545 tolerance_rel = 1e-9 energyE = np.sum(scc.epsilon_0/2*(Ex**2+Ey**2+Ez**2)) |