Use previous criterion to update EB

ECP-WarpX · Jan 18, 2025 · 6fa77bd · 6fa77bd
1 parent 4969716
commit 6fa77bd
Showing 1 changed file with 6 additions and 15 deletions.
diff --git a/Source/Initialization/WarpXInitData.cpp b/Source/Initialization/WarpXInitData.cpp
@@ -1264,22 +1264,13 @@ void WarpX::InitializeEBGridData (int lev)
                 // 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);
-            } else {
-                // Mark on which grid points E should be updated (stair-case approximation)
-                MarkUpdateCellsStairCase(
-                    m_eb_update_E[lev],
-                    m_fields.get_alldirs(FieldType::Efield_fp, lev),
-                    eb_fact );
-                // Mark on which grid points B should be updated (stair-case approximation)
-                MarkUpdateCellsStairCase(
-                    m_eb_update_B[lev],
-                    m_fields.get_alldirs(FieldType::Bfield_fp, lev),
-                    eb_fact );
             }
+
+            // 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);
+
         }
 
         ComputeDistanceToEB();