Method interp

ibtk/src/math/PatchMathOps.cpp:3936–4009

} // grad

void
PatchMathOps::interp(Pointer<CellData<NDIM, double>> dst,
                     const Pointer<FaceData<NDIM, double>> src,
                     const Pointer<Patch<NDIM>> patch) const
{
    const int U_ghosts = (dst->getGhostCellWidth()).max();
    const int v_ghosts = (src->getGhostCellWidth()).max();

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

#if !defined(NDEBUG)
    if (dst->getDepth() != NDIM * src->getDepth())
    {
        TBOX_ERROR("PatchMathOps::interp():\n"
                   << "  src and dst have incompatible depths" << std::endl);
    }

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

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

    if (patch_box != dst->getBox())
    {
        TBOX_ERROR("PatchMathOps::interp():\n"
                   << "  dst and src must live on the same patch" << std::endl);
    }

    if (patch_box != src->getBox())
    {
        TBOX_ERROR("PatchMathOps::interp():\n"
                   << "  dst and src must live on the same patch" << std::endl);
    }
#endif

    for (int depth = 0; depth < src->getDepth(); ++depth)
    {
        // Interpolate.
        double* const U = dst->getPointer(NDIM * depth);

        const double* const v0 = src->getPointer(0, depth);
        const double* const v1 = src->getPointer(1, depth);
#if (NDIM == 3)
        const double* const v2 = src->getPointer(2, depth);
#endif

        F_TO_C_INTERP_FC(U,
                         U_ghosts,
                         v0,
                         v1,
#if (NDIM == 3)