This made me wonder if RLibm has been extended to doubles yet[0]. Sadly that doesn't seem to be the case. It's a floating point math library that always produce

This made me wonder if RLibm has been extended to doubles yet[0]. Sadly that doesn't seem to be the case. It's a floating point math library that always produces correct rounding results, but it is limited to 32-bit floats (and posits, for that matter). Interestingly, despite being the only library that always rounds correctly, it also claims to be faster than other libraries (at same bit-depth). I can recommend the papers, they're quite accessible and there is something very satisfying about the math tricks they came up with to build the polynomials. They do seem to have added more rounding modes though, which presumably makes the library more generally useful.[1] [0] https://people.cs.rutgers.edu/~sn349/rlibm/ https://people.cs.rutgers.edu/~sn349/rlibm/ [1] https://blog.sigplan.org/2022/04/28/one-polynomial-approximation-to-produce-correctly-rounded-results-for-multiple-representations-and-rounding-modes/ https://blog.sigplan.org/2022/04/28/one-polynomial-approxima...