Method vc_laplace

ibtk/src/math/PatchMathOps.cpp:5847–6023

} // laplace

void
PatchMathOps::vc_laplace(Pointer<SideData<NDIM, double>> dst,
                         const double alpha,
                         const double beta,
                         const Pointer<NodeData<NDIM, double>> coef1,
                         const Pointer<SideData<NDIM, double>> coef2,
                         const Pointer<SideData<NDIM, double>> src1,
                         const double gamma_in,
                         const Pointer<SideData<NDIM, double>> src2_in,
                         const Pointer<Patch<NDIM>> patch,
                         const bool use_harmonic_interp,
                         const int l,
                         const int m,
                         const int n) const
{
#if (NDIM == 2)
    const Pointer<CartesianPatchGeometry<NDIM>> pgeom = patch->getPatchGeometry();
    const double* const dx = pgeom->getDx();

    double* const f0 = dst->getPointer(0, l);
    double* const f1 = dst->getPointer(1, l);

    const int f_ghosts = (dst->getGhostCellWidth()).max();

    const double* const mu = coef1->getPointer();
    const int mu_ghosts = (coef1->getGhostCellWidth()).max();

    const int rho_varying = !(coef2.isNull());
    const double* rho0 = nullptr;
    const double* rho1 = nullptr;
    const int rho_ghosts = rho_varying ? (coef2->getGhostCellWidth()).max() : 0;
    if (rho_varying)
    {
        rho0 = coef2->getPointer(0, m);
        rho1 = coef2->getPointer(1, m);
    }

    const double* const u0 = src1->getPointer(0, m);
    const double* const u1 = src1->getPointer(1, m);
    const int u_ghosts = (src1->getGhostCellWidth()).max();

    const double gamma = (src2_in ? gamma_in : 0.0);
    const Pointer<SideData<NDIM, double>> src2 = (src2_in ? src2_in : src1);
    const double* const v0 = src2->getPointer(0, n);
    const double* const v1 = src2->getPointer(1, n);
    const int v_ghosts = (src2->getGhostCellWidth()).max();

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

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

    if (mu_ghosts != (coef1->getGhostCellWidth()).min())
    {