Method initializeHierarchyIndependentData

src/IB/IBInstrumentPanel.cpp:574–687

} // isInstrumented

void
IBInstrumentPanel::initializeHierarchyIndependentData(const Pointer<PatchHierarchy<NDIM>> hierarchy,
                                                      LDataManager* const l_data_manager)
{
    IBAMR_TIMER_START(t_initialize_hierarchy_independent_data);

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

    // Determine how many flow meters/pressure gauges are present in the local
    // data.
    int max_meter_index = -1;
    std::vector<int> max_node_index;
    for (int ln = coarsest_ln; ln <= finest_ln; ++ln)
    {
        if (l_data_manager->levelContainsLagrangianData(ln))
        {
            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 m = spec->getMeterIndex();
                    max_meter_index = std::max(m, max_meter_index);

                    const int n = spec->getNodeIndex();
                    max_node_index.resize(max_meter_index + 1, -1);
                    max_node_index[m] = std::max(n, max_node_index[m]);
                }
            }
        }
    }

    // Communicate local data to all processes.
    d_num_meters = IBTK_MPI::maxReduction(max_meter_index) + 1;
    max_node_index.resize(d_num_meters, -1);
    d_num_perimeter_nodes.clear();
    d_num_perimeter_nodes.resize(d_num_meters, -1);
    for (unsigned int m = 0; m < d_num_meters; ++m)
    {
        d_num_perimeter_nodes[m] = max_node_index[m] + 1;
    }
    IBTK_MPI::maxReduction(d_num_meters > 0 ? &d_num_perimeter_nodes[0] : nullptr, d_num_meters);
#if !defined(NDEBUG)
    for (unsigned int m = 0; m < d_num_meters; ++m)
    {
        TBOX_ASSERT(d_num_perimeter_nodes[m] > 0);
    }
#endif

    // Resize arrays.
    d_X_centroid.resize(d_num_meters);
    d_X_perimeter.resize(d_num_meters);
    for (unsigned int m = 0; m < d_num_meters; ++m)
    {
        d_X_perimeter[m].resize(boost::extents[d_num_perimeter_nodes[m]]);