Some updates

Author: WeiqunZhang

Docs/sphinx_documentation/source/FFT.rst

@@ -113,16 +113,16 @@ boundaries. Its implementation utilizes :cpp:`FFT::OpenBCSolver`, which can
 be used for implementing convolution based solvers with a user provided
 Green's function. If users want to extend the open BC solver to 2D or other
 types of Green's function, they could use :cpp:`FFT::PoissonOpenBC` as an
-example. Below is an example of solving the Poisson equation with open
+example. Below is an example of solving Poisson's equation with open
 boundaries.
 
 .. highlight:: c++
 
 ::
 
     Geometry geom(...);
-    MultiFab soln(...); // soln & rhs can be either nodal or cell-centered
-    MultiFab rhs(...);
+    MultiFab soln(...); // soln can be either nodal or cell-centered.
+    MultiFab rhs(...);  // rhs must have the same index type as soln.
 
     int ng = ...; // ng can be non-zero, if we want to compute potential
                   // outside the domain.

Src/Base/AMReX_BoxArray.cpp

@@ -2052,15 +2052,15 @@ BoxArray decompose (Box const& domain, int nboxes,
                 if (b.ok()) {
                     if (no_overlap) {
                         for (int idim = 0; idim < AMREX_SPACEDIM; ++idim) {
-                            if (ixtyp.nodeCentered() &&
+                            if (ixtyp.nodeCentered(idim) &&
                                 b.bigEnd(idim) == ccdomain.bigEnd(idim))
                             {
                                 b.growHi(idim, 1);
                             }
                         }
                     } else {
                         for (int idim = 0; idim < AMREX_SPACEDIM; ++idim) {
-                            if (ixtyp.nodeCentered()) {
+                            if (ixtyp.nodeCentered(idim)) {
                                 b.growHi(idim, 1);
                             }
                         }

Tests/FFT/OpenBC/main.cpp

@@ -5,8 +5,6 @@
 #include <AMReX_MultiFabUtil.H>
 #include <AMReX_ParmParse.H>
 
-#include <AMReX_VisMF.H>
-
 using namespace amrex;
 
 int main (int argc, char* argv[])