Точность превыше всего!

[green_fr]

На работе использую программу для генерирования экономических сценариев (B+H ESG), при определённой конфигурации она вылетает с нечитаемой ошибкой. Hotline объясняет, что у меня данные неправильные, моя матрица должна быть положительно полуопределённой, бла-бла-бла.

Я проверяю — одно из собственных значений выходит −2E-16.
Формально да, отрицательное. Но мы все понимаем, что это ноль, проблема округления (более того, я могу с ручкой на бумажке доказать, что это ноль, я намеренно сделал вырожденную матрицу, мне так надо).

Весь день бодаюсь теперь с hotline’ом по почте, доказывая, что это не «extremely high degree of accuracy in calculations» (цитата из их ответа), а баг, требующий патча. Ищу понятные аргументы...

Update: победили, ошибку признали, постараются когда-нибудь починить, но уже точно не успеют к следующей (7.3.0) версии.

Comments

[birdwatcher.livejournal.com]

Ну сделай им вырожденную с положительным собственным числом.

[green-fr.livejournal.com]

В смысле? Она уже вырожденная, кроме этих "нулей" все остальные уже числа положительны.

[birdwatcher.livejournal.com]

Я имел в виду то же самое, что french_man. Думаю, достаточно сдвинуть любое число.

[green-fr.livejournal.com]

Ну да, пока оставляю это на крайний случай. Потому что мне потом же надо будет ещё обосновывать эту матрицу, у нас и так официальная документация пестрит «ввиду технических ограничений», «допущение для большей скорости вычислений» и пр...