What an interesting piece!
Every developer should know that floating point numbers implement approximations.
During my master thesys in Civil Enginineering (twenty years ago!) I learned that professor Jonathan R Shewchuk wrote almost perfect C code about this, and used in its Delaunay triangulation. At that time I hadn’t the time to really study it as it seemed to me a little like magic.