2016-09-03 9:48 GMT+02:00 Nicolas Cellier <nicolas.cellier.aka.nice@gmail.com>:

For floating points it would be good to also have something related to unit of least precision like I think it exists in google tests:

assert: aFloat isWithin: anInteger ulpFrom: anotherFloat

������ ^(aFloat - anotherFloat) abs <= (anotherFloat ulp max: aFloat ulp)

It's testing if the result are the same with a tolerance on least significant bits.
(anotherFloat ulp max: aFloat ulp) is not really necessary, it's only for making the comparison symmetric.

indeed, 2 and 2 predecessor do not have the same exponent, so do not have the same ulp.

If we don't have symmetry, then

������ self assert: 2.0 predecessor predecessor isWithin: 1 ulpFrom: 2.0.
might differ from:
������ self assert: 2.0 isWithin: 1 ulpFrom: 2.0 predecessor predecessor .

Note that the this message cannot be used to test for zero result (because 0.0 ulp is the smallest denormal).

For example:

������ self assert: 2.0 - 2.0 predecessor isWithin: 4 ulpFrom: 0.0.

would anwer false, but it makes no sense!

One should write

������ self assert: 2.0 isWithin: 4 ulpFrom: 2.0 predecessor.

Bad result will be in other cases too (not only zero):

expected := 8085- 8084.9999.

computed := (80.85 * 100) - 8084.9999.

(computed - expected) abs <= (computed ulp max: expected ulp) �� ��"==> false"