diff options
| author |  Lorenzo Giacomel <47607756+lgiacome@users.noreply.github.com>
| 2021-11-16 18:22:20 +0100 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2021-11-16 17:22:20 +0000 |
| commit | f515975a6ec2680fd9aada32aee00857b25588c9 (patch) | |
| tree | eeea5ff3e8d0b3c48f56cd52008ffeafcb7e42a9 /Source/Initialization/WarpXInitData.cpp | |
| parent | 0e058ef2644f5420dc3e12cccd52782945e00975 (diff) | |
| download | WarpX-f515975a6ec2680fd9aada32aee00857b25588c9.tar.gz WarpX-f515975a6ec2680fd9aada32aee00857b25588c9.tar.zst WarpX-f515975a6ec2680fd9aada32aee00857b25588c9.zip | |
Fixes to the EB init (#2565)
Diffstat (limited to 'Source/Initialization/WarpXInitData.cpp')
| -rw-r--r-- | Source/Initialization/WarpXInitData.cpp | 10 |
1 files changed, 7 insertions, 3 deletions
diff --git a/Source/Initialization/WarpXInitData.cpp b/Source/Initialization/WarpXInitData.cpp index be8225c81..59b7d1f25 100644 --- a/Source/Initialization/WarpXInitData.cpp +++ b/Source/Initialization/WarpXInitData.cpp @@ -486,13 +486,14 @@ WarpX::InitLevelData (int lev, Real /*time*/) #ifdef AMREX_USE_EB if(lev==maxLevel()) { - if (WarpX::maxwell_solver_id == MaxwellSolverAlgo::Yee - || WarpX::maxwell_solver_id == MaxwellSolverAlgo::ECT) { + if (WarpX::maxwell_solver_id == MaxwellSolverAlgo::Yee || + WarpX::maxwell_solver_id == MaxwellSolverAlgo::CKC || + WarpX::maxwell_solver_id == MaxwellSolverAlgo::ECT) { + ComputeEdgeLengths(); ComputeFaceAreas(); ScaleEdges(); ScaleAreas(); - ComputeDistanceToEB(); const auto &period = Geom(lev).periodicity(); WarpXCommUtil::FillBoundary(*m_edge_lengths[lev][0], guard_cells.ng_alloc_EB, period); @@ -516,6 +517,9 @@ WarpX::InitLevelData (int lev, Real /*time*/) ComputeFaceExtensions(); } } + + ComputeDistanceToEB(); + } #endif |