use POSIX;
use constant epsilon => 100 * DBL_EPSILON;