I have been asked to test a library provided by a 3rd party. The library is known to be accurate to n significant figures. Any less-significant errors can safely be ignored. I want to write a function to help me compare the results:
def nearlyequal(a, b, sigfig=5):
The purpose of this function is to determine if two floating-point numbers (a and b) are approximately equal. The function will return True if a==b (exact match) or if a and b have the same value when rounded to sigfig significant-figures when written in decimal.
Can anybody suggest a good implementation? I've written a mini unit-test. Unless you can see a bug in my tests then a good implementation should pass the following:
assert nearlyequal(1, 1, 5)
assert nearlyequal(1.0, 1.0, 5)
assert nearlyequal(1.0, 1.0, 5)
assert nearlyequal(-1e-9, 1e-9, 5)
assert nearlyequal(1e9, 1e9 + 1 , 5)
assert not nearlyequal(1e4, 1e4 + 1, 5)
assert nearlyequal(0.0, 1e-15, 5)
assert not nearlyequal(0.0, 1e-4, 6)
Additional notes:
- Values a and b might be of type int, float or numpy.float64. Values a and b will always be of the same type. It's vital that conversion does not introduce additional error into the function.
- Lets keep this numerical, so functions that convert to strings or use non-mathematical tricks are not ideal. This program will be audited by somebody who is a mathematician who will want to be able to prove that the function does what it is supposed to do.
- Speed... I've got to compare a lot of numbers so the faster the better.
- I've got numpy, scipy and the standard-library. Anything else will be hard for me to get, especially for such a small part of the project.
0.0001
is 4 digits. It can also be used to mean digits in the mantissa when written in scientific notation, i.e.1.001e-10
. I believe the latter is more common usage. These are not the same. It's easier to specify if you are talking about digits to the right of the decimal point. In that case, you meanfabs(x-y)<0.0001
means x is close to y. Otherwise, comparing to 0 is difficult.