Method initializePatchHierarchy

src/IB/CIBMethod.cpp:506–623

} // initializeLevelData

void
CIBMethod::initializePatchHierarchy(Pointer<PatchHierarchy<NDIM>> hierarchy,
                                    Pointer<GriddingAlgorithm<NDIM>> gridding_alg,
                                    int u_data_idx,
                                    const std::vector<Pointer<CoarsenSchedule<NDIM>>>& u_synch_scheds,
                                    const std::vector<Pointer<RefineSchedule<NDIM>>>& u_ghost_fill_scheds,
                                    int integrator_step,
                                    double init_data_time,
                                    bool initial_time)
{
    // Initialize various Lagrangian data objects required by the conventional
    // IB method.
    IBMethod::initializePatchHierarchy(hierarchy,
                                       gridding_alg,
                                       u_data_idx,
                                       u_synch_scheds,
                                       u_ghost_fill_scheds,
                                       integrator_step,
                                       init_data_time,
                                       initial_time);

    // Set structure index info.
    const int struct_ln = getStructuresLevelNumber();
    std::vector<int> structIDs = d_l_data_manager->getLagrangianStructureIDs(struct_ln);
    std::sort(structIDs.begin(), structIDs.end());
    const auto structs_on_this_ln = static_cast<unsigned>(structIDs.size());

    for (unsigned struct_no = 0; struct_no < structs_on_this_ln; ++struct_no)
    {
        d_struct_lag_idx_range[struct_no] =
            d_l_data_manager->getLagrangianStructureIndexRange(structIDs[struct_no], struct_ln);
    }

    // Initialize Lagrangian and Eulerian lambda at initial time.
    const int coarsest_ln = 0;
    const int finest_ln = d_hierarchy->getFinestLevelNumber();
    if (initial_time)
    {
        // Initialize Lagrangian lambda.
        setInitialLambda(finest_ln);

        // Initialize Eulerian lambda (S[lambda]) variable.
        for (int ln = coarsest_ln; ln <= finest_ln; ++ln)
        {
            Pointer<PatchLevel<NDIM>> level = d_hierarchy->getPatchLevel(ln);
            for (PatchLevel<NDIM>::Iterator p(level); p; p++)
            {
                Pointer<Patch<NDIM>> patch = level->getPatch(p());
                Pointer<CellData<NDIM, double>> lambda_data = patch->getPatchData(d_eul_lambda_idx);
                lambda_data->fillAll(0.0);
            }
        }
    }

    // Register plot quantities.
    if (d_silo_writer)
    {
        for (int ln = coarsest_ln; ln <= finest_ln; ++ln)