Refactor IEEECK to avoid generating IEEE exceptions by using IEEE_ARITHMETIC

[foxtran]

Description

I've replaced division-by-zero and other invalid floating-point operations with IEEE_ARITHMETIC module from modern Fortran. It should simplify debugging of numerical instabilities of code that also uses LAPACK.

Fixes #1172

Checklist

• [x] The documentation has been updated: no need to update documentation. Keep old documentation as is for compability with previous behaviour.
• [x] If the PR solves a specific issue, it is set to be closed on merge.