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
------------------------------------------------------------------------------
BlackBerry® DevCon Americas, Oct. 18-20, San Francisco, CA
Learn about the latest advances in developing for the 
BlackBerry® mobile platform with sessions, labs & more.
See new tools and technologies. Register for BlackBerry® DevCon today!
http://p.sf.net/sfu/rim-devcon-copy1 
_______________________________________________
Libmesh-devel mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-devel

Reply via email to