From ca67ae0700d1f56e7921da9283d42184149cb15b Mon Sep 17 00:00:00 2001 From: Lorenzo Giacomel <47607756+lgiacome@users.noreply.github.com> Date: Thu, 13 Jan 2022 20:05:19 +0100 Subject: Fixing staircased EM solver (#2739) * Fixing the staircase consistency * Removed the face_areas multifabs everywhere they're not needed * Bug fix * More fixes * Another fix * Another fix * Initialize areas anyways for the initialization * [pre-commit.ci] auto fixes from pre-commit.com hooks for more information, see https://pre-commit.ci Co-authored-by: lgiacome Co-authored-by: pre-commit-ci[bot] <66853113+pre-commit-ci[bot]@users.noreply.github.com> --- Source/Initialization/WarpXInitData.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'Source/Initialization/WarpXInitData.cpp') diff --git a/Source/Initialization/WarpXInitData.cpp b/Source/Initialization/WarpXInitData.cpp index be810d21e..95efa48aa 100644 --- a/Source/Initialization/WarpXInitData.cpp +++ b/Source/Initialization/WarpXInitData.cpp @@ -888,8 +888,8 @@ void WarpX::InitializeEBGridData (int lev) auto const eb_fact = fieldEBFactory(lev); ComputeEdgeLengths(m_edge_lengths[lev], eb_fact); - ComputeFaceAreas(m_face_areas[lev], eb_fact); ScaleEdges(m_edge_lengths[lev], CellSize(lev)); + ComputeFaceAreas(m_face_areas[lev], eb_fact); ScaleAreas(m_face_areas[lev], CellSize(lev)); if (WarpX::maxwell_solver_id == MaxwellSolverAlgo::ECT) { -- cgit v1.2.3