Method initializeHierarchyDependentData

src/IB/IBInstrumentPanel.cpp:689–914

} // initializeHierarchyIndependentData

void
IBInstrumentPanel::initializeHierarchyDependentData(const Pointer<PatchHierarchy<NDIM>> hierarchy,
                                                    LDataManager* const l_data_manager,
                                                    const int timestep_num,
                                                    const double data_time)
{
    if (!d_initialized)
    {
        initializeHierarchyIndependentData(hierarchy, l_data_manager);
    }
    if (d_num_meters == 0) return;

    IBAMR_TIMER_START(t_initialize_hierarchy_dependent_data);

    const int coarsest_ln = 0;
    const int finest_ln = hierarchy->getFinestLevelNumber();

    // Keep track of the timestep number and time at which the instrument data
    // was collected.
    d_instrument_read_timestep_num = timestep_num;
    d_instrument_read_time = data_time;

    // Loop over all local nodes to determine the positions of the local
    // perimeter nodes.
    for (unsigned int m = 0; m < d_num_meters; ++m)
    {
        for (int n = 0; n < d_num_perimeter_nodes[m]; ++n)
        {
            d_X_perimeter[m][n] = Point::Zero();
        }
    }
    for (int ln = coarsest_ln; ln <= finest_ln; ++ln)
    {
        if (l_data_manager->levelContainsLagrangianData(ln))
        {
            // Extract the local position array.
            Pointer<LData> X_data = l_data_manager->getLData(LDataManager::POSN_DATA_NAME, ln);
            Vec X_vec = X_data->getVec();
            double* X_arr;
            int ierr = VecGetArray(X_vec, &X_arr);
            IBTK_CHKERRQ(ierr);

            // Store the local positions of the perimeter nodes.
            const Pointer<LMesh> mesh = l_data_manager->getLMesh(ln);
            const std::vector<LNode*>& local_nodes = mesh->getLocalNodes();
            for (const auto& node_idx : local_nodes)
            {
                const IBInstrumentationSpec* const spec = node_idx->getNodeDataItem<IBInstrumentationSpec>();
                if (spec)
                {
                    const int& petsc_idx = node_idx->getLocalPETScIndex();
                    const double* const X = &X_arr[NDIM * petsc_idx];
                    const int m = spec->getMeterIndex();
                    const int n = spec->getNodeIndex();
                    std::copy(X, X + NDIM, d_X_perimeter[m][n].data());
                }
            }