Точность текста с плавающей запятой Python / numpy

Python 2.7 и более поздние repr уже имеют интеллектуальную реализацию для float, которая печатает 0,1 как 0.1 . Краткий вывод выбирается предпочтительнее других кандидатов, таких как 0.10000000000000001 потому что это кратчайшее представление этого конкретного числа, которое округляет до того же значения с плавающей запятой при чтении обратно на Python. Чтобы использовать этот алгоритм, конвертируйте свои 64-битные поплавки в реальные поплавки Python, прежде чем передавать их в repr :

 >>> map(repr, map(float, v64)) ['5.0', '0.1', '2.4', '4.555555555555555', '12345678.923456786'] 

Удивительно, но результат естественный и численно правильный. Более подробную информацию о версии 2.7 / 3.2 можно найти в « Что нового» и увлекательной лекции Марка Дикинсона.

К сожалению, этот трюк не будет работать для 32-битных поплавков, по крайней мере, не переосмысливая алгоритм, используемый в repr Python 2.7.

Чтобы однозначно определить число с плавающей запятой с одной точностью (32-разрядной) в формате IEEE-754, может потребоваться использовать 9 (значимое, т.е. не начинающееся с 0, если значение равно 0) десятичные цифры, а 9 цифр всегда достаточно.

Для чисел с плавающей запятой с двойной точностью (64 бит) может понадобиться 17 (значащих) десятичных цифр и всегда достаточны.

Я не совсем уверен, как определяется формат %g , он может позволить представлению начинаться с 0 (0,1), поэтому безопасные значения для точности будут .9 и .17 .

Если вы хотите минимизировать размер файла, запись байтовых представлений приведет к созданию гораздо меньшего файла, поэтому, если вы можете это сделать, это путь.

Код C, который реализует фантастическое представление в версии 2.7, в основном состоит из Python / dtoa.c (с обертками в Python / pystrtod.c и Objects / floatobject.c). В частности, посмотрите на _Py_dg_dtoa. Должна быть возможность заимствовать этот код и модифицировать его для работы с float вместо double. Затем вы можете обернуть это в модуль расширения или просто создать его так и ctypes.

Кроме того, обратите внимание, что источник говорит, что реализация «вдохновлена« Как печатать цифры с плавающей точкой »Томаном Ги Слилем-младшим и Джоном Л. Уайтом [Proc. ACM SIGPLAN '90, стр. 112-126] «. Таким образом, вы можете реализовать что-то менее гибкое и более простое, прочитав эту бумагу (и какой бы ни был modfications, описанный в комментариях dtoa.c, кажется уместным).

Наконец, код является незначительным изменением кода, представленным Дэвидом Гей в AT & T, и используется в ряде других библиотек (NSPR и т. Д.), Один из которых может иметь более доступную версию.

Но прежде чем делать что-либо из этого, убедитесь, что на самом деле проблема с производительностью, попробовав функцию Python и измеряя, слишком ли она слишком медленная.

И если это действительно критически важная для вас область, вы, вероятно, не захотите перебирать список и сначала вызвать функцию (или собственную собственную функцию C); вам, вероятно, нужна функция, которая преобразует массив numpy с плавающей точкой или удваивается в строку сразу. (В идеале вы хотели бы построить это в numpy, конечно.)

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

Излишне говорить, что вы также должны проверить, достаточно ли %.9g и %.17g для ваших целей, поскольку это следующая самая легкая вещь, которая может работать.