diff options
| author |  Dave Grote <grote1@llnl.gov>
| 2019-04-26 21:30:35 -0700 |
|---|---|---|
| committer | Dave Grote <grote1@llnl.gov>
| 2019-04-26 21:30:35 -0700 |
| commit | d528e233a7710dc497174c825620f1de7323faa5 (patch) | |
| tree | 9e8d02d4bc84b3b1cfcc8403f7feddcac3a7a352 /Source/WarpX.cpp | |
| parent | 1d4e9797cb022491378034253c3019babb822e83 (diff) | |
| download | WarpX-d528e233a7710dc497174c825620f1de7323faa5.tar.gz WarpX-d528e233a7710dc497174c825620f1de7323faa5.tar.zst WarpX-d528e233a7710dc497174c825620f1de7323faa5.zip | |
Implemented dive cleaning for RZ geometry
Diffstat (limited to 'Source/WarpX.cpp')
| -rw-r--r-- | Source/WarpX.cpp | 18 |
1 files changed, 16 insertions, 2 deletions
diff --git a/Source/WarpX.cpp b/Source/WarpX.cpp index 51d9133be..c0212d4b3 100644 --- a/Source/WarpX.cpp +++ b/Source/WarpX.cpp @@ -957,12 +957,19 @@ WarpX::ComputeDivE (MultiFab& divE, int dcomp, for (MFIter mfi(divE, true); mfi.isValid(); ++mfi) { const Box& bx = mfi.tilebox(); +#ifdef WARPX_RZ + const Real xmin = bx.smallEnd(0)*dx[0]; +#endif WRPX_COMPUTE_DIVE(bx.loVect(), bx.hiVect(), BL_TO_FORTRAN_N_ANYD(divE[mfi],dcomp), BL_TO_FORTRAN_ANYD((*E[0])[mfi]), BL_TO_FORTRAN_ANYD((*E[1])[mfi]), BL_TO_FORTRAN_ANYD((*E[2])[mfi]), - dx.data()); + dx.data() +#ifdef WARPX_RZ + ,&xmin +#endif + ); } } @@ -977,12 +984,19 @@ WarpX::ComputeDivE (MultiFab& divE, int dcomp, for (MFIter mfi(divE, true); mfi.isValid(); ++mfi) { Box bx = mfi.growntilebox(ngrow); +#ifdef WARPX_RZ + const Real xmin = bx.smallEnd(0)*dx[0]; +#endif WRPX_COMPUTE_DIVE(bx.loVect(), bx.hiVect(), BL_TO_FORTRAN_N_ANYD(divE[mfi],dcomp), BL_TO_FORTRAN_ANYD((*E[0])[mfi]), BL_TO_FORTRAN_ANYD((*E[1])[mfi]), BL_TO_FORTRAN_ANYD((*E[2])[mfi]), - dx.data()); + dx.data() +#ifdef WARPX_RZ + ,&xmin +#endif + ); } } |