Fast_math: A collection of functions for fast number crunching using Fortran

[FedericoPerini]

@rwmsu

I can only guess that the old “backward compatability”

Not an expert, but I would bet that ensuring IEEE arithmetics compliance for all edge cases is also part of the answer.