Re: [Scheme-reports] Proposed language for 'eqv?' applied to inexact real numbers Mark H Weaver 12 Nov 2012 17:53 UTC

I wrote:
> Noah Lavine <noah.b.lavine@gmail.com> writes:
>>     The 'eqv?' predicate on two elements of S must return #t if its
>> arguments are the same member of S, and #f otherwise. Note that a
>> single member of S may have different representations, but arithmetic
>> operations are defined on the abstract set S and not on the
>> representations."
>>
>> The goal is basically to push parts of the definition of eqv? down to
>> implementations, but do it in a structured way. This would require
>> that (eqv? 1.0 1.0) => #t,
>
> It is a mistake to try to ensure that (eqv? x y) => #true when x and y
> represent the same numeric value in different representations.  You are
> trying to make 'eqv?' act like '=', and that is a fundamental error, and
> a very common one.

I meant to include a specific example of why (eqv? 1.0 1.0) should
sometimes return #false, if those two arguments are represented
differently.

Suppose those two arguments are represented with different precisions.
For example, one might be a 64-bit IEEE double and the other might be a
32-bit float, or perhaps an MPFR big float.  Now consider the following
predicate:

  (define (foo? x) (zero? (/ x (expt 10 400))))

If one uses the usual precision-propagation rules, then (foo? 1.0) will
return #true if its argument is an IEEE double or float, but might
return #false if its argument is a MPFR big float.

Therefore, 'eqv?' must distinguish between inexact numbers of different
precisions, otherwise memoization does not work.

Consider a square root procedure, where the precision of the result
depends on the precision of its argument.  So (sqrt 2.0) will return a
more precise answer if its argument if an MPFR big float.  This
procedure can only be properly memoized if 'eqv?' distinguishes numbers
of different precisions.

     Mark

_______________________________________________
Scheme-reports mailing list
Scheme-reports@scheme-reports.org
http://lists.scheme-reports.org/cgi-bin/mailman/listinfo/scheme-reports