[Libmesh-devel] Total Linear Its During Nonlinear Solve

[Derek G. <fri...@gm...>]

A while back we changed NonlinearSolver::solve() to return a std::pair that included the total number of nonlinear iterations and the final nonlinear residual.

I'm now also wanting the total number of _linear_ iterations.

What are we going to do here long term?

This used to not be a problem because we didn't call clear() after the solve (like we now do) so you could just get the snes object out and inspect it after a solve… but that is no longer possible.

Should I just modify the interface so it returns a struct?  Or should we use Boost::tuple?  Or should we save off this data as class members that you can query with something like NonlinearSolver::total_nonlinear_iterations(), NonlinearSolver::total_linear_iterations() ???  Kind of like we have get_converged_reason in PetscNonlinearSolver ?

For now I think I'm going to add a get_total_linear_iterations() just to PetscNonlinearSolver so I can get my work done…. it will follow the path of get_converged_reason()….

Derek