Method rot

ibtk/src/math/PatchMathOps.cpp:2146–2256

} // curl

void
PatchMathOps::rot(Pointer<SideData<NDIM, double>> dst,
                  const Pointer<NodeData<NDIM, double>> src,
                  const Pointer<Patch<NDIM>> patch,
                  CartSideRobinPhysBdryOp* bc_op,
                  const double fill_time) const
{
#if (NDIM != 2)
    TBOX_ERROR("PatchMathOps::rot():\n"
               << "  not implemented for NDIM != 2" << std::endl);
    NULL_USE(dst);
    NULL_USE(src);
    NULL_USE(patch);
    NULL_USE(bc_op);
    NULL_USE(fill_time);
#endif
#if (NDIM == 2)
    const Pointer<CartesianPatchGeometry<NDIM>> pgeom = patch->getPatchGeometry();
    const double* const dx = pgeom->getDx();

    double* const w0 = dst->getPointer(0);
    double* const w1 = dst->getPointer(1);
    const int w_ghosts = (dst->getGhostCellWidth()).max();

    const double* const u0 = src->getPointer(0);
    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::rot():\n"
                   << "  dst does not have uniform ghost cell widths" << std::endl);
    }

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

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

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

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