code: static private boolean PE_EQ(double a, double b) {
  		return (a == b)
  				|| PE_ABS(a - b) <= PE_EPS * (1 + (PE_ABS(a) + PE_ABS(b)) / 2);
  	}
docs: Determine if two doubles are equal within a default tolerance