Добрый день! Вопрос к тем, кто имеет опыт работы (или разбирался с исходниками) с б-кой gsl (gnu scientific library,
http://www.gnu.org/software/gsl), в частности с решателем систем нелинейных алг. ур-й. Мы тут использовали решатель gsl_multiroot_fsolver_dnewton (discrete Newton algorithm без использования матрицы производных). В решаемых задачах сходимость наблюдалась не во всех случаях: у некоторых задач вообще всё время расходилось (хотя решение существует, и довольно близко), у других удалось добиться сходимости изменением параметров задачи.
Вопрос в том, является ли этот решатель чистым методом Ньютона, или в него ещё встроены какие-то улучшатели сходимости? В овервью написано, что "Additional strategies can be used to enlarge the region of convergence" (
http://www.gnu.org/software/gsl/manual/ ... nding.html), но в самом описании пользовательского интерфейса (
http://www.gnu.org/software/gsl/manual/ ... tives.html) никаких настроек по этому поводу не упомянуто.