Method initializeSolverState

src/IB/CIBMobilitySolver.cpp:258–305

} // getMobilitySolvers

void
CIBMobilitySolver::initializeSolverState(Vec x, Vec b)
{
    IBAMR_TIMER_START(t_initialize_solver_state);

    // Deallocate the solver state if the solver is already initialized.
    if (d_is_initialized)
    {
        d_reinitializing_solver = true;
        deallocateSolverState();
    }

    if (d_mobility_solver_type == KRYLOV)
    {
        d_krylov_mob_solver->initializeSolverState(x, b);
    }
    else if (d_mobility_solver_type == DIRECT)
    {
        d_direct_mob_solver->initializeSolverState(x, b);
    }
    else
    {
        TBOX_ERROR("CIBMobilitySolver::initializeSolverState() Unknown mobility solver type" << std::endl);
    }

    d_has_free_parts = false;
    for (unsigned part = 0; part < d_num_rigid_parts && !d_has_free_parts; ++part)
    {
        int num_free_dofs;
        d_cib_strategy->getSolveRigidBodyVelocity(part, num_free_dofs);
        if (num_free_dofs)
        {
            d_has_free_parts = true;
        }
    }
    if (d_has_free_parts)
    {
        d_krylov_freebody_mob_solver->initializeSolverState(x, b);
    }

    // Indicate that the solver is initialized.
    d_reinitializing_solver = false;
    d_is_initialized = true;

    IBAMR_TIMER_STOP(t_initialize_solver_state);

    return;
} // initializeSolverState

void
CIBMobilitySolver::deallocateSolverState()