Method solveSystem

src/IB/CIBSaddlePointSolver.cpp:433–473

} // getCIBMobilitySolver

bool
CIBSaddlePointSolver::solveSystem(Vec x, Vec b)
{
    IBTK_TIMER_START(t_solve_system);

    // Initialize the solver, when necessary.
    const bool deallocate_after_solve = !d_is_initialized;
    if (deallocate_after_solve) initializeSolverState(x, b);

    d_petsc_x = x;
    VecCopy(b, d_petsc_b);

    // Modify RHS for inhomogeneous Bcs.
    d_A->setHomogeneousBc(false);
    d_A->modifyRhsForBcs(d_petsc_b);
    d_A->setHomogeneousBc(true);

    // Solve the system.
    KSPSolve(d_petsc_ksp, d_petsc_b, d_petsc_x);
    KSPGetIterationNumber(d_petsc_ksp, &d_current_iterations);
    KSPGetResidualNorm(d_petsc_ksp, &d_current_residual_norm);

    // Impose Solution Bcs.
    d_A->setHomogeneousBc(false);
    d_A->imposeSolBcs(d_petsc_x);

    // Determine the convergence reason.
    KSPConvergedReason reason;
    KSPGetConvergedReason(d_petsc_ksp, &reason);
    const bool converged = (static_cast<int>(reason) > 0);
    if (d_enable_logging) reportPETScKSPConvergedReason(d_object_name, reason, plog);

    // Invalidate d_petsc_x Vec.
    d_petsc_x = nullptr;

    // Deallocate the solver, when necessary.
    if (deallocate_after_solve) deallocateSolverState();

    IBTK_TIMER_STOP(t_solve_system);
    return converged;
} // solveSystem

void
CIBSaddlePointSolver::initializeSolverState(Vec x, Vec b)