Replace i1mach and d1mach calls

This pull request replaces the calls to `i1mach` and `d1mach` in the recently added polygamma code. I didn't see a cleaner way to replace the code to determine the number of digits in the significand of a Float64 value.