Hi, I have a remark on "lines_search()" method in "NewtonSolver" class.
The method doesn't update "current_residual" since it is passed to it by value. The consequence is that the second call of "test_convergence()" in "NewtonSolver::solve()" is done using wrong value of "current_residual". Regards, Mladen ------------------------------------------------------------------------------ Open Source Business Conference (OSBC), March 24-25, 2009, San Francisco, CA -OSBC tackles the biggest issue in open source: Open Sourcing the Enterprise -Strategies to boost innovation and cut costs with open source participation -Receive a $600 discount off the registration fee with the source code: SFAD http://p.sf.net/sfu/XcvMzF8H _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
