Я просмотрел некоторую документацию о том, как работают функции numpy/scipy с точки зрения численной стабильности, например. любые средства, используемые для улучшения численной устойчивости или существуют альтернативные стабильные реализации.
Меня особенно интересует оператор (+
) массивов с плавающей запятой, numpy.sum()
, numpy.cumsum()
и numpy.dot()
. Во всех случаях я суммирую очень большое количество чисел с плавающей запятой, и меня беспокоит точность таких вычислений.
Кто-нибудь знает какую-либо ссылку на такие проблемы в документации numpy/scipy или на какой-то другой источник?