[Python] How do I access rounding modes defined in IEEE 754 in python?

I have the need to compute exactly in single precision floating points in Python.

The options I've tried are decimal.Decimal and numpy.float32. However Decimal is not based on IEEE 754, and float32 cannot allow the use of rounding modes. It is surprising rounding modes are standard features of IEEE 754 but there exists no built-in implementation of it in Python.

Continue reading...