diff options
| author |  Lorenzo Giacomel <47607756+lgiacome@users.noreply.github.com>
| 2022-01-13 20:05:19 +0100 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2022-01-13 11:05:19 -0800 |
| commit | ca67ae0700d1f56e7921da9283d42184149cb15b (patch) | |
| tree | 6ee96802136c95bf7fbbbebbf19a871516ea697a /Source/BoundaryConditions/PML.cpp | |
| parent | 726a9d85b951cf26f936ef0dd600a342657f4106 (diff) | |
| download | WarpX-ca67ae0700d1f56e7921da9283d42184149cb15b.tar.gz WarpX-ca67ae0700d1f56e7921da9283d42184149cb15b.tar.zst WarpX-ca67ae0700d1f56e7921da9283d42184149cb15b.zip | |
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 <lorenzo.giacome@cern.ch>
Co-authored-by: pre-commit-ci[bot] <66853113+pre-commit-ci[bot]@users.noreply.github.com>
Diffstat (limited to '')
| -rw-r--r-- | Source/BoundaryConditions/PML.cpp | 13 |
1 files changed, 0 insertions, 13 deletions
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<MultiFab>(amrex::convert( ba, WarpX::GetInstance().getEfield_fp(0,2).ixType().toIntVect() ), dm, WarpX::ncomps, max_guard_EB ); - pml_face_areas[0] = std::make_unique<MultiFab>(amrex::convert( ba, - WarpX::GetInstance().getBfield_fp(0,0).ixType().toIntVect() ), dm, WarpX::ncomps, max_guard_EB ); - pml_face_areas[1] = std::make_unique<MultiFab>(amrex::convert( ba, - WarpX::GetInstance().getBfield_fp(0,1).ixType().toIntVect() ), dm, WarpX::ncomps, max_guard_EB ); - pml_face_areas[2] = std::make_unique<MultiFab>(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<amrex::Real,3> 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<MultiFab*,3> -PML::Get_face_areas() -{ - return {pml_face_areas[0].get(), pml_face_areas[1].get(), pml_face_areas[2].get()}; -} MultiFab* PML::GetF_fp () |