Я использую порт minpack hybrd1
в Eigen , который использует метод Пауэлла , чтобы найти решение f (x) = 0; в этом случае якобиан вычисляется численно.
Решатель принимает параметр xtol
, относительный допуск, который по умолчанию sqrt(eps)
для данного числового типа (double
в моем случае)
В большинстве случаев я получаю удовлетворительное решение с векторной нормой остатков, например 1.446575e-9
. Иногда, однако, решение дает норму, например 523.543, что заставляет меня задуматься о том, что на самом деле означает xtol
. Какова контрольная величина, с которой сравнивается текущая норма остатков?