Method postprocessIntegrateData

src/IB/IBMethod.cpp:416–483

} // preprocessIntegrateData

void
IBMethod::postprocessIntegrateData(double current_time, double new_time, int /*num_cycles*/)
{
    int ierr;
    const int coarsest_ln = 0;
    const int finest_ln = d_hierarchy->getFinestLevelNumber();
    const double dt = new_time - current_time;
    const int integrator_step = d_ib_solver->getIntegratorStep();

    // Update the instrumentation data.
    updateIBInstrumentationData(integrator_step + 1, new_time);
    if (d_instrument_panel->isInstrumented())
    {
        const std::vector<std::string>& instrument_name = d_instrument_panel->getInstrumentNames();
        const std::vector<double>& flow_data = d_instrument_panel->getFlowValues();
        for (unsigned int m = 0; m < flow_data.size(); ++m)
        {
            // NOTE: Flow volume is calculated in default units.
            d_total_flow_volume[m] += flow_data[m] * dt;
            if (d_do_log)
                plog << "flow volume through " << instrument_name[m] << ":\t " << d_total_flow_volume[m] << "\n";
        }
    }

    // Reset time-dependent Lagrangian data.
    for (int ln = coarsest_ln; ln <= finest_ln; ++ln)
    {
        if (!d_l_data_manager->levelContainsLagrangianData(ln)) continue;
        ierr = VecSwap(d_X_current_data[ln]->getVec(), d_X_new_data[ln]->getVec());
        IBTK_CHKERRQ(ierr);
        ierr = VecSwap(d_U_current_data[ln]->getVec(), d_U_new_data[ln]->getVec());
        IBTK_CHKERRQ(ierr);
        if (d_F_new_data[ln])
        {
            ierr = VecSwap(d_F_current_data[ln]->getVec(), d_F_new_data[ln]->getVec());
            IBTK_CHKERRQ(ierr);
        }
        else if (d_F_half_data[ln])
        {
            ierr = VecSwap(d_F_current_data[ln]->getVec(), d_F_half_data[ln]->getVec());
            IBTK_CHKERRQ(ierr);
        }
    }
    d_X_current_needs_ghost_fill = true;
    d_F_current_needs_ghost_fill = true;

    // Deallocate Lagrangian scratch data.
    d_X_current_data.clear();
    d_X_new_data.clear();
    d_X_half_data.clear();
    d_X_jac_data.clear();
    d_X_LE_new_data.clear();
    d_X_LE_half_data.clear();
    d_U_current_data.clear();
    d_U_new_data.clear();
    d_U_half_data.clear();
    d_U_jac_data.clear();
    d_F_current_data.clear();