Some time ago I've used a simple algorithm by using 128bit floats for doubles (and 64bit for floats) which seems to work very nicely and is straightforward to implement. It passed the full tests for conversion for 32bit floats and the sparse tests for doubles.

I did a blog post about it: https://www.fixscript.org/blog/math-library (includes interactive demos)

Any opinions?