Значение относительной толерантности в процедуре минимизации (minpack's hybrd1) - PullRequest
0 голосов
/ 29 июля 2011

Я использую порт minpack hybrd1 в Eigen , который использует метод Пауэлла , чтобы найти решение f (x) = 0; в этом случае якобиан вычисляется численно.

Решатель принимает параметр xtol, относительный допуск, который по умолчанию sqrt(eps) для данного числового типа (double в моем случае)

В большинстве случаев я получаю удовлетворительное решение с векторной нормой остатков, например 1.446575e-9. Иногда, однако, решение дает норму, например 523.543, что заставляет меня задуматься о том, что на самом деле означает xtol. Какова контрольная величина, с которой сравнивается текущая норма остатков?

1 Ответ

0 голосов
/ 03 августа 2011

Как я выяснил при проверке кода, xtol относится к норме начального предположения (предоставленного пользователем), где запускается решатель.

...