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 index d45887ca932..118214948a5 100644 --- 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 } } \ No newline at end of file diff --git a/Source/Initialization/WarpXInitData.cpp b/Source/Initialization/WarpXInitData.cpp index 4db7096eb05..248448e817e 100644 --- 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(