I did a course at university on mathematical stability in floating point algorithms which mainly taught me there's no way I was smart enough to write correct code.

Plus then making them work optimally on N different microarchitectures.

I actually wonder what language/system they were using that lacked this in the stdlib or at least in very widely available libraries.