From 6fa77bd5f5a1ac20ad7eec70dbeea08afeb796fa Mon Sep 17 00:00:00 2001 From: Remi Lehe Date: Fri, 17 Jan 2025 06:16:41 -0800 Subject: [PATCH] Use previous criterion to update EB --- Source/Initialization/WarpXInitData.cpp | 21 ++++++--------------- 1 file changed, 6 insertions(+), 15 deletions(-) diff --git a/Source/Initialization/WarpXInitData.cpp b/Source/Initialization/WarpXInitData.cpp index 248448e817e..147d2bfd0de 100644 --- 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();