Fix bug with ES solver and MR: E_aux=E_fp in UpdateAuxilaryData

Examples/Tests/electrostatic_sphere_eb/analysis_rz.py

@@ -31,39 +31,52 @@
 fn = sys.argv[1]
 ds = yt.load( fn )
 
-all_data_level_0 = ds.covering_grid(level=0,left_edge=ds.domain_left_edge, dims=ds.domain_dimensions)
-phi = all_data_level_0['boxlib', 'phi'].v.squeeze()
-Er = all_data_level_0['boxlib', 'Er'].v.squeeze()
+def get_error_per_lev(ds, level):
+    if level == 1:
+        dims = ds.domain_dimensions * ds.refine_by**level
+        all_data_level = ds.covering_grid(level=level, left_edge=ds.domain_left_edge, dims=dims)
+        phi = all_data_level['boxlib', 'phi'].v.squeeze()
+        phi = phi[:,:,0]
+        Er = all_data_level['boxlib', 'Er'].v.squeeze()
+        Er = Er[:,:,0]
+    elif level==0:
+        dims = ds.domain_dimensions * ds.refine_by**level
+        all_data_level = ds.covering_grid(level=level, left_edge=ds.domain_left_edge, dims=dims)
+        phi = all_data_level['boxlib', 'phi'].v.squeeze()
+        Er = all_data_level['boxlib', 'Er'].v.squeeze()
 
-Dx = ds.domain_width/ds.domain_dimensions
-dr = Dx[0]
-rmin = ds.domain_left_edge[0]
-rmax = ds.domain_right_edge[0]
-nr = phi.shape[0]
-r = np.linspace(rmin+dr/2.,rmax-dr/2.,nr)
-B = 1.0/np.log(0.1/0.5)
-A = -B*np.log(0.5)
+    Dx = ds.domain_width/dims
+    dr = Dx[0]
+    rmin = ds.domain_left_edge[0]
+    rmax = ds.domain_right_edge[0]
+    nr = phi.shape[0]
+    r = np.linspace(rmin+dr/2.,rmax-dr/2.,nr)
+    B = 1.0/np.log(0.1/0.5)
+    A = -B*np.log(0.5)
+    err = 0.0
+    errmax_phi = 0.0
+    errmax_Er = 0.0
+    for i in range(len(r)):
+        # outside EB and last cutcell
+        if r[i] > 0.1*m + dr:
+            phi_theory = A+B*np.log(r[i])
+            Er_theory = -B/float(r[i])
+            err = abs( phi_theory - phi[i,:] ).max() / phi_theory
+            if err>errmax_phi:
+                errmax_phi = err
+            err = abs( Er_theory - Er[i,:] ).max() / Er_theory
+            # Exclude the last inaccurate interpolation.
+            if err>errmax_Er and i<len(r)-1:
+                errmax_Er = err
 
-err = 0.0
-errmax_phi = 0.0
-errmax_Er = 0.0
-for i in range(len(r)):
-    # outside EB and last cutcell
-    if r[i] > 0.1*m + dr:
-        phi_theory = A+B*np.log(r[i])
-        Er_theory = -B/float(r[i])
-        err = abs( phi_theory - phi[i,:] ).max() / phi_theory
-        if err>errmax_phi:
-            errmax_phi = err
-        err = abs( Er_theory - Er[i,:] ).max() / Er_theory
-        # Exclude the last inaccurate interpolation.
-        if err>errmax_Er and i<len(r)-1:
-            errmax_Er = err
+    print('max error of phi = ', errmax_phi)
+    print('max error of Er = ', errmax_Er)
+    print('tolerance = ', tolerance)
+    assert(errmax_phi<tolerance and errmax_Er<tolerance)
 
-print('max error of phi = ', errmax_phi)
-print('max error of Er = ', errmax_Er)
-print('tolerance = ', tolerance)
-assert(errmax_phi<tolerance and errmax_Er<tolerance)
+levels = [0,1]
+for level in levels:
+    get_error_per_lev(ds,level)
 
 test_name = os.path.split(os.getcwd())[1]
 checksumAPI.evaluate_checksum(test_name, fn, do_particles=False)

Regression/Checksum/benchmarks_json/ElectrostaticSphereEB_RZ_MR.json

@@ -4,7 +4,7 @@
     "phi": 2036.0428085225362
   },
   "lev=1": {
-    "Er": 19519.172334977942,
-    "phi": 3291.0262856782897
+    "Er": 9563.742949788306,
+    "phi": 3290.9609834276607
   }
-}
+}

Regression/Checksum/benchmarks_json/ElectrostaticSphereLabFrame_MR_emass_10.json

@@ -6,15 +6,15 @@
     "rho": 0.0
   },
   "lev=1": {
-    "Ex": 14.281015560380963,
-    "Ey": 14.281015560380965,
-    "Ez": 14.281015560380965,
+    "Ex": 7.170105936287823,
+    "Ey": 7.17010593628782,
+    "Ez": 7.170105936287821,
     "rho": 2.6092568008333786e-10
   },
   "electron": {
-    "particle_momentum_x": 1.80842228672388e-24,
-    "particle_momentum_y": 1.8084222867238806e-24,
-    "particle_momentum_z": 1.7598771525647628e-24,
+    "particle_momentum_x": 9.257577597262615e-25,
+    "particle_momentum_y": 9.257577597262618e-25,
+    "particle_momentum_z": 9.257577597262624e-25,
     "particle_position_x": 327.46875,
     "particle_position_y": 327.46875,
     "particle_position_z": 327.46875,

Regression/WarpX-tests.ini

@@ -554,7 +554,7 @@ analysisRoutine = Examples/Tests/electrostatic_sphere_eb/analysis_rz.py
 [ElectrostaticSphereEB_RZ_MR]
 buildDir = .
 inputFile = Examples/Tests/electrostatic_sphere_eb/inputs_rz_mr
-runtime_params =  warpx.abort_on_warning_threshold = medium
+runtime_params =  warpx.abort_on_warning_threshold = medium warpx.fine_tag_lo = 0.0  -0.5 warpx.fine_tag_hi = 0.5  0.5 amr.ref_ratio_vect = 2 2 2
 dim = 2
 addToCompileString = USE_EB=TRUE USE_RZ=TRUE
 cmakeSetupOpts = -DWarpX_DIMS=RZ -DWarpX_EB=ON