Я столкнулся с предел точности mpmath с плавающей запятой?

Question

mpmath претендует на поддержку «арифметики с произвольной точностью с плавающей запятой».

И все же. , ,

>>> import mpmath 
>>> 1 + mpmath.erf(-5.921) 
mpf('1.1102230246251565e-16') 
>>> 1 + mpmath.erf(-5.922) # I expect a smaller positive number here. 
mpf('0.0') 

Я что-то не хватает? Или это фундаментальное ограничение mpmath?

@jonrsharpe предположил, что проблема в том, что я отправил float в erf. Тем не менее, приведенный ниже код показывает, что это не проблема:

>>> 1 + mpmath.erf(mpmath.mpf('-5.922')) 
mpf('0.0') 

Answer 1

Вопрос в данном случае имеет отношение к mpmath «s global precision setting слишком низко. Значение по умолчанию для prec является

>>> mpmath.mp.prec 
53 

Когда я установил его 100, я получаю результат я ожидал:

>>> 1 + mpmath.erf(-5.922) 
mpf('5.5236667058718205581661131647751e-17') 

В этом случае разница в скорости не заметно, но отмечают, что увеличение точность обычно увеличивает время, необходимое для вычисления результата.