Method curl

ibtk/src/math/PatchMathOps.cpp:1497–1592

/////////////////////////////// PUBLIC ///////////////////////////////////////

void
PatchMathOps::curl(Pointer<CellData<NDIM, double>> dst,
                   const Pointer<CellData<NDIM, double>> src,
                   const Pointer<Patch<NDIM>> patch) const
{
    const Pointer<CartesianPatchGeometry<NDIM>> pgeom = patch->getPatchGeometry();
    const double* const dx = pgeom->getDx();

    double* const W = dst->getPointer();
    const int W_ghosts = (dst->getGhostCellWidth()).max();

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

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

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

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

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

    const Box<NDIM>& U_box = src->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::curl():\n"
                   << "  src has insufficient ghost cell width" << std::endl);
    }

    const int W_depth = dst->getDepth();

    if (
#if (NDIM == 2)
        (W_depth != 1)
#elif (NDIM == 3)
        (W_depth != NDIM)
#else
        false
#endif
    )
    {
        TBOX_ERROR("PatchMathOps::curl():\n"
                   << "  dst has incorrect depth" << std::endl);