Method laplace

ibtk/src/math/PatchMathOps.cpp:5039–5208

} // harmonic_interp

void
PatchMathOps::laplace(Pointer<CellData<NDIM, double>> dst,
                      const double alpha,
                      const double beta,
                      const Pointer<CellData<NDIM, double>> src1,
                      const double gamma,
                      const Pointer<CellData<NDIM, double>> src2,
                      const Pointer<Patch<NDIM>> patch,
                      const int l,
                      const int m,
                      const int n) const
{
    const Pointer<CartesianPatchGeometry<NDIM>> pgeom = patch->getPatchGeometry();
    const double* const dx = pgeom->getDx();

    double* const F = dst->getPointer(l);
    const int F_ghosts = (dst->getGhostCellWidth()).max();

    const double* const U = src1->getPointer(m);
    const int U_ghosts = (src1->getGhostCellWidth()).max();

    const Box<NDIM>& patch_box = patch->getBox();

#if !defined(NDEBUG)
    if (F_ghosts != (dst->getGhostCellWidth()).min())
    {
        TBOX_ERROR("PatchMathOps::laplace():\n"
                   << "  dst does not have uniform ghost cell widths" << std::endl);
    }

    if (U_ghosts != (src1->getGhostCellWidth()).min())
    {
        TBOX_ERROR("PatchMathOps::laplace():\n"
                   << "  src1 does not have uniform ghost cell widths" << std::endl);
    }

    if (src1 == dst)
    {
        TBOX_ERROR("PatchMathOps::laplace():\n"
                   << "  src1 == dst." << std::endl);
    }

    if ((src1 == src2) && (gamma != 0.0))
    {
        TBOX_ERROR("PatchMathOps::laplace():\n"
                   << "  src1 == src2 but gamma is nonzero." << std::endl);
    }

    const Box<NDIM>& U_box = src1->getGhostBox();
    const Box<NDIM> U_box_shrunk = Box<NDIM>::grow(U_box, -1);

    if ((!U_box_shrunk.contains(patch_box.lower())) || (!U_box_shrunk.contains(patch_box.upper())))
    {
        TBOX_ERROR("PatchMathOps::laplace():\n"
                   << "  src1 has insufficient ghost cell width" << std::endl);
    }

    if (patch_box != dst->getBox())