Method interpolateVelocity

src/IB/IBMethod.cpp:671–694

} // updateFixedLEOperators

void
IBMethod::interpolateVelocity(const int u_data_idx,
                              const std::vector<Pointer<CoarsenSchedule<NDIM>>>& u_synch_scheds,
                              const std::vector<Pointer<RefineSchedule<NDIM>>>& u_ghost_fill_scheds,
                              const double data_time)
{
    std::vector<Pointer<LData>>*U_data, *X_LE_data;
    bool* X_LE_needs_ghost_fill;
    getVelocityData(&U_data, data_time);
    getLECouplingPositionData(&X_LE_data, &X_LE_needs_ghost_fill, data_time);
    d_l_data_manager->interp(u_data_idx, *U_data, *X_LE_data, u_synch_scheds, u_ghost_fill_scheds, data_time);
    resetAnchorPointValues(*U_data,
                           /*coarsest_ln*/ 0,
                           /*finest_ln*/ d_hierarchy->getFinestLevelNumber());

    if (!IBTK::rel_equal_eps(data_time, d_half_time))
    {
        std::vector<Pointer<LData>>* U_half_data;
        getVelocityData(&U_half_data, d_half_time);
        reinitMidpointData(d_U_current_data, d_U_new_data, *U_half_data);
    }

    return;
} // interpolateVelocity

void
IBMethod::interpolateLinearizedVelocity(const int u_data_idx,