Modify initialization in ECT solver

ECP-WarpX · Jan 18, 2025 · ace399b · ace399b
1 parent 0980516
commit ace399b
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 9 deletions.
diff --git a/Regression/Checksum/benchmarks_json/test_3d_embedded_boundary_rotated_cube.json b/Regression/Checksum/benchmarks_json/test_3d_embedded_boundary_rotated_cube.json
@@ -1,10 +1,10 @@
 {
   "lev=0": {
-    "Bx": 1.252616939910365e-05,
-    "By": 0.02473895628331097,
-    "Bz": 0.024738956316621142,
-    "Ex": 10253221.850298548,
-    "Ey": 10387.334582977643,
-    "Ez": 10387.532806510022
+    "Bx": 1.280747509243305e-05,
+    "By": 2.473900144296397e-02,
+    "Bz": 2.473890786894079e-02,
+    "Ex": 1.025322901921306e+07,
+    "Ey": 1.042254197269831e+04,
+    "Ez": 1.040011664019071e+04
   }
 }
diff --git a/Source/Initialization/WarpXInitData.cpp b/Source/Initialization/WarpXInitData.cpp
@@ -1261,13 +1261,13 @@ void WarpX::InitializeEBGridData (int lev)
             ScaleAreas(face_areas_lev, CellSize(lev));
 
             if (WarpX::electromagnetic_solver_id == ElectromagneticSolverAlgo::ECT) {
+                // Compute additional quantities required for the ECT solver
+                MarkExtensionCells();
+                ComputeFaceExtensions();
                 // Mark on which grid points E should be updated
                 MarkUpdateECellsECT( m_eb_update_E[lev], edge_lengths_lev );
                 // Mark on which grid points B should be updated
                 MarkUpdateBCellsECT( m_eb_update_B[lev], face_areas_lev, edge_lengths_lev);
-                // Compute additional quantities required for the ECT solver
-                MarkExtensionCells();
-                ComputeFaceExtensions();
             } else {
                 // Mark on which grid points E should be updated (stair-case approximation)
                 MarkUpdateCellsStairCase(