Method copyFluidVariable

src/IB/ConstraintIBMethod.cpp:2347–2386

} // midpointStep

void
ConstraintIBMethod::copyFluidVariable(int copy_from_idx, int copy_to_idx)
{
    const int coarsest_ln = 0;
    const int finest_ln = d_hierarchy->getFinestLevelNumber();

    for (int ln = coarsest_ln; ln <= finest_ln; ++ln)
    {
        Pointer<PatchLevel<NDIM>> level = d_hierarchy->getPatchLevel(ln);
        if (!level->checkAllocated(copy_to_idx)) level->allocatePatchData(copy_to_idx);
    }

    SAMRAIVectorReal<NDIM, double> u_from(d_object_name + "from", d_hierarchy, coarsest_ln, finest_ln);
    SAMRAIVectorReal<NDIM, double> u_to(d_object_name + "to", d_hierarchy, coarsest_ln, finest_ln);

    u_from.addComponent(d_u_fluidSolve_var, copy_from_idx, d_wgt_sc_idx);
    u_to.addComponent(d_u_fluidSolve_var, copy_to_idx, d_wgt_sc_idx);

    u_to.copyVector(Pointer<SAMRAIVectorReal<NDIM, double>>(&u_from, false));

    using InterpolationTransactionComponent = IBTK::HierarchyGhostCellInterpolation::InterpolationTransactionComponent;
    std::vector<InterpolationTransactionComponent> transaction_comps;
    InterpolationTransactionComponent component(copy_to_idx,
                                                DATA_REFINE_TYPE,
                                                USE_CF_INTERPOLATION,
                                                SIDE_DATA_COARSEN_TYPE,
                                                BDRY_EXTRAP_TYPE,
                                                CONSISTENT_TYPE_2_BDRY,
                                                std::vector<SAMRAI::solv::RobinBcCoefStrategy<NDIM>*>(NDIM, nullptr),
                                                nullptr);
    transaction_comps.push_back(component);

    Pointer<HierarchyGhostCellInterpolation> hier_bdry_fill = new HierarchyGhostCellInterpolation();
    hier_bdry_fill->initializeOperatorState(transaction_comps, d_hierarchy, coarsest_ln, finest_ln);
    const bool homogeneous_bc = true;
    hier_bdry_fill->setHomogeneousBc(homogeneous_bc);
    hier_bdry_fill->fillData(0.0);

    return;
} // copyFluidVariable

void
ConstraintIBMethod::copyDensityVariable(int copy_from_idx, int copy_to_idx)