Replace diff with epsilon test

[Jacob Williams]

The diff function is a trick to check against the machine epsilon. Replace it with an actual test, and then we can eliminate this trick.

See: https://fortran-lang.discourse.group/t/looking-at-some-old-code/4580