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/BoundaryConditions/PML.cpp | 13 ------------- 1 file changed, 13 deletions(-) (limited to 'Source/BoundaryConditions/PML.cpp') diff --git a/Source/BoundaryConditions/PML.cpp b/Source/BoundaryConditions/PML.cpp index 6be5dbc1b..1241c864e 100644 --- a/Source/BoundaryConditions/PML.cpp +++ b/Source/BoundaryConditions/PML.cpp @@ -681,12 +681,6 @@ PML::PML (const int lev, const BoxArray& grid_ba, const DistributionMapping& gri WarpX::GetInstance().getEfield_fp(0,1).ixType().toIntVect() ), dm, WarpX::ncomps, max_guard_EB ); pml_edge_lengths[2] = std::make_unique(amrex::convert( ba, WarpX::GetInstance().getEfield_fp(0,2).ixType().toIntVect() ), dm, WarpX::ncomps, max_guard_EB ); - pml_face_areas[0] = std::make_unique(amrex::convert( ba, - WarpX::GetInstance().getBfield_fp(0,0).ixType().toIntVect() ), dm, WarpX::ncomps, max_guard_EB ); - pml_face_areas[1] = std::make_unique(amrex::convert( ba, - WarpX::GetInstance().getBfield_fp(0,1).ixType().toIntVect() ), dm, WarpX::ncomps, max_guard_EB ); - pml_face_areas[2] = std::make_unique(amrex::convert( ba, - WarpX::GetInstance().getBfield_fp(0,2).ixType().toIntVect() ), dm, WarpX::ncomps, max_guard_EB ); if (WarpX::maxwell_solver_id == MaxwellSolverAlgo::Yee || WarpX::maxwell_solver_id == MaxwellSolverAlgo::CKC || @@ -695,10 +689,8 @@ PML::PML (const int lev, const BoxArray& grid_ba, const DistributionMapping& gri auto const eb_fact = fieldEBFactory(); WarpX::ComputeEdgeLengths(pml_edge_lengths, eb_fact); - WarpX::ComputeFaceAreas(pml_face_areas, eb_fact); std::array cellsize = {AMREX_D_DECL(geom->CellSize()[0],geom->CellSize()[1],geom->CellSize()[2])}; WarpX::ScaleEdges(pml_edge_lengths, cellsize); - WarpX::ScaleAreas(pml_face_areas, cellsize); } #endif @@ -1056,11 +1048,6 @@ PML::Get_edge_lengths() return {pml_edge_lengths[0].get(), pml_edge_lengths[1].get(), pml_edge_lengths[2].get()}; } -std::array -PML::Get_face_areas() -{ - return {pml_face_areas[0].get(), pml_face_areas[1].get(), pml_face_areas[2].get()}; -} MultiFab* PML::GetF_fp () -- cgit v1.2.3