[tlb]

I wrote an interval arithmetic package once too. It was slow, because it had to change the FP rounding flags multiple times for some operations.

In the end, it seemed like any substantial computation ended up having extremely wide bounds, much wider than they deserved. Trying to invert a matrix often resulted in [-Inf .. +Inf] bounds.