[postgis-devel] PostGIS release

[Michael Fuhr]

On Thu, Jul 12, 2007 at 07:49:31PM -0600, Michael Fuhr wrote:
> ./pgis-test 16384.0001 16384.0001 16385
> FP_CONTAINS_INCL(16384.000100, 16384.000100, 16385.000000) = 0
> 
> When A gets large enough, A - PGIS_EPSILON requires more precision
> than available to be distinguishable from A, so A - PGIS_EPSILON
> equals A and is no longer less than B.

I get the above results on the various platforms that I've already
mentioned.  Here's what I get on a different system (Linux i686):

./pgis-test 16384.0001 16384.0001 16385
FP_CONTAINS_INCL(16384.000100, 16384.000100, 16385.000000) = 1

./pgis-test 33554432 33554432 40000000
FP_CONTAINS_INCL(33554432.000000, 33554432.000000, 40000000.000000) = 1

./pgis-test 33554433 33554433 40000000
FP_CONTAINS_INCL(33554433.000000, 33554433.000000, 40000000.000000) = 0

-- 
Michael Fuhr