Method initializeDirectorDataOnPatchLevel

src/IB/IBRedundantInitializer.cpp:1172–1216

} // initializeMassOnPatchLevel

unsigned int
IBRedundantInitializer::initializeDirectorDataOnPatchLevel(const unsigned int /*global_index_offset*/,
                                                           const unsigned int local_index_offset,
                                                           Pointer<LData> D_data,
                                                           const Pointer<PatchHierarchy<NDIM>> hierarchy,
                                                           const int level_number,
                                                           const double /*init_data_time*/,
                                                           const bool /*can_be_refined*/,
                                                           const bool /*initial_time*/,
                                                           LDataManager* const /*l_data_manager*/)
{
#if !defined(NDEBUG)
    TBOX_ASSERT(d_data_processed);
#endif

    // Loop over all patches in the specified level of the patch level and
    // initialize the local vertices.
    boost::multi_array_ref<double, 2>& D_array = *D_data->getLocalFormVecArray();
    int local_idx = invalid_index;
    int local_node_count = 0;
    Pointer<PatchLevel<NDIM>> level = hierarchy->getPatchLevel(level_number);
    for (PatchLevel<NDIM>::Iterator p(level); p; p++)
    {
        Pointer<Patch<NDIM>> patch = level->getPatch(p());

        // Initialize the vertices whose initial locations will be within the
        // given patch.
        std::vector<std::pair<int, int>> patch_vertices;
        getPatchVertices(patch_vertices, patch, hierarchy);
        local_node_count += patch_vertices.size();
        for (const auto& point_idx : patch_vertices)
        {
            const int local_petsc_idx = ++local_idx + local_index_offset;

            // Initialize the director corresponding to the present vertex.
            const std::vector<double>& D = getVertexDirectors(point_idx, level_number);
            for (int d = 0; d < 3 * 3; ++d)
            {
                D_array[local_petsc_idx][d] = D[d];
            }
        }
    }
    D_data->restoreArrays();
    return local_node_count;
} // initializeDirectorOnPatchLevel

void
IBRedundantInitializer::tagCellsForInitialRefinement(const Pointer<PatchHierarchy<NDIM>> hierarchy,