diff options
Diffstat (limited to 'Source/Diagnostics/FullDiagnostics.cpp')
| -rw-r--r-- | Source/Diagnostics/FullDiagnostics.cpp | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/Source/Diagnostics/FullDiagnostics.cpp b/Source/Diagnostics/FullDiagnostics.cpp index c272643ad..26d30ebce 100644 --- a/Source/Diagnostics/FullDiagnostics.cpp +++ b/Source/Diagnostics/FullDiagnostics.cpp @@ -469,10 +469,10 @@ FullDiagnostics::InitializeBufferData (int i_buffer, int lev ) { // To ensure that the diagnostic lo and hi are within the domain defined at level, lev. diag_dom.setLo(idim, std::max(m_lo[idim],warpx.Geom(lev).ProbLo(idim)) ); diag_dom.setHi(idim, std::min(m_hi[idim],warpx.Geom(lev).ProbHi(idim)) ); - if ( fabs(warpx.Geom(lev).ProbLo(idim) - diag_dom.lo(idim)) + if ( std::fabs(warpx.Geom(lev).ProbLo(idim) - diag_dom.lo(idim)) > warpx.Geom(lev).CellSize(idim) ) use_warpxba = false; - if ( fabs(warpx.Geom(lev).ProbHi(idim) - diag_dom.hi(idim)) + if ( std::fabs(warpx.Geom(lev).ProbHi(idim) - diag_dom.hi(idim)) > warpx.Geom(lev).CellSize(idim) ) use_warpxba = false; @@ -493,11 +493,11 @@ FullDiagnostics::InitializeBufferData (int i_buffer, int lev ) { amrex::IntVect hi(1); for (int idim=0; idim < AMREX_SPACEDIM; ++idim) { // lo index with same cell-size as simulation at level, lev. - lo[idim] = std::max( static_cast<int>( floor ( + lo[idim] = std::max( static_cast<int>( std::floor ( ( diag_dom.lo(idim) - warpx.Geom(lev).ProbLo(idim)) / warpx.Geom(lev).CellSize(idim)) ), 0 ); // hi index with same cell-size as simulation at level, lev. - hi[idim] = std::max( static_cast<int> ( ceil ( + hi[idim] = std::max( static_cast<int> ( std::ceil ( ( diag_dom.hi(idim) - warpx.Geom(lev).ProbLo(idim)) / warpx.Geom(lev).CellSize(idim) ) ), 0) - 1 ; // if hi<=lo, then hi = lo + 1, to ensure one cell in that dimension |