spfpm is a pure-Python toolkit for performing binary fixed-point arithmetic, including trigonometric and exponential functions.
The package provides:
On a modern desktop PC, spfpm is typically capable of hundreds of thousands of arithmetic operations per second, i.e. over 100 kilo-FLOPS, even for a few hundred bits of resolution. As a pure-Python library SPFPM offers portability and interactivity, but will never compete with the performance of lower-level libraries such as mpmath, Boost multiprecision or GMP for heavyweight calculations.
Development currently targets Python versions 3.3 and later, although the library may also be usable with python-2.7. The latest version of spfpm can be found on GitHub.
After installation there are two main classes that you need from the FixedPoint module:
from FixedPoint import FXfamily, FXnum
you can create fixed-point numbers with the default (64-bit) resolution as follows:
x = FXnum(22) / FXnum(7)
y = FXnum(3.1415)
print(x - y)
Creating numbers with a specific precision requires use
of the FXfamily
class:
fam100 = FXfamily(100)
z = FXnum(1, fam100)
z2 = fam100(2)
One can then apply various computations such as:
print(z.atan() * 4)
print(z2.sqrt())
The FXfamily
class also provides access to pre-computed constants
which should be accurate to about 1/2 of the least significant bit (LSB):
print('pi = ', fam100.pi)
print('log2 = ', fam100.log2)
print('sqrt2 = ', fam100.sqrt2)
This produces the following printed values of those constants:
The struck-through digits show where the values computed at 100-bit precision differ from the true digits in the decimal representations of these constants. Alternatively, one could print these values in base-16:
print(FXfamily(400).pi.toBinaryString(logBase=4))
giving a hexadecimal value of Pi as:
which agrees exactly with the accepted result.
All files are released under the Python PSF License and are Copyright 2006-2024 RW Penney.