2000-01-05 07:06:23 +01:00
|
|
|
--
|
|
|
|
-- POINT
|
|
|
|
--
|
Change floating-point output format for improved performance.
Previously, floating-point output was done by rounding to a specific
decimal precision; by default, to 6 or 15 decimal digits (losing
information) or as requested using extra_float_digits. Drivers that
wanted exact float values, and applications like pg_dump that must
preserve values exactly, set extra_float_digits=3 (or sometimes 2 for
historical reasons, though this isn't enough for float4).
Unfortunately, decimal rounded output is slow enough to become a
noticable bottleneck when dealing with large result sets or COPY of
large tables when many floating-point values are involved.
Floating-point output can be done much faster when the output is not
rounded to a specific decimal length, but rather is chosen as the
shortest decimal representation that is closer to the original float
value than to any other value representable in the same precision. The
recently published Ryu algorithm by Ulf Adams is both relatively
simple and remarkably fast.
Accordingly, change float4out/float8out to output shortest decimal
representations if extra_float_digits is greater than 0, and make that
the new default. Applications that need rounded output can set
extra_float_digits back to 0 or below, and take the resulting
performance hit.
We make one concession to portability for systems with buggy
floating-point input: we do not output decimal values that fall
exactly halfway between adjacent representable binary values (which
would rely on the reader doing round-to-nearest-even correctly). This
is known to be a problem at least for VS2013 on Windows.
Our version of the Ryu code originates from
https://github.com/ulfjack/ryu/ at commit c9c3fb1979, but with the
following (significant) modifications:
- Output format is changed to use fixed-point notation for small
exponents, as printf would, and also to use lowercase 'e', a
minimum of 2 exponent digits, and a mandatory sign on the exponent,
to keep the formatting as close as possible to previous output.
- The output of exact midpoint values is disabled as noted above.
- The integer fast-path code is changed somewhat (since we have
fixed-point output and the upstream did not).
- Our project style has been largely applied to the code with the
exception of C99 declaration-after-statement, which has been
retained as an exception to our present policy.
- Most of upstream's debugging and conditionals are removed, and we
use our own configure tests to determine things like uint128
availability.
Changing the float output format obviously affects a number of
regression tests. This patch uses an explicit setting of
extra_float_digits=0 for test output that is not expected to be
exactly reproducible (e.g. due to numerical instability or differing
algorithms for transcendental functions).
Conversions from floats to numeric are unchanged by this patch. These
may appear in index expressions and it is not yet clear whether any
change should be made, so that can be left for another day.
This patch assumes that the only supported floating point format is
now IEEE format, and the documentation is updated to reflect that.
Code by me, adapting the work of Ulf Adams and other contributors.
References:
https://dl.acm.org/citation.cfm?id=3192369
Reviewed-by: Tom Lane, Andres Freund, Donald Dong
Discussion: https://postgr.es/m/87r2el1bx6.fsf@news-spur.riddles.org.uk
2019-02-13 16:20:33 +01:00
|
|
|
-- avoid bit-exact output here because operations may not be bit-exact.
|
|
|
|
SET extra_float_digits = 0;
|
2000-01-05 07:06:23 +01:00
|
|
|
CREATE TABLE POINT_TBL(f1 point);
|
|
|
|
INSERT INTO POINT_TBL(f1) VALUES ('(0.0,0.0)');
|
|
|
|
INSERT INTO POINT_TBL(f1) VALUES ('(-10.0,0.0)');
|
|
|
|
INSERT INTO POINT_TBL(f1) VALUES ('(-3.0,4.0)');
|
|
|
|
INSERT INTO POINT_TBL(f1) VALUES ('(5.1, 34.5)');
|
|
|
|
INSERT INTO POINT_TBL(f1) VALUES ('(-5.0,-12.0)');
|
2018-09-26 10:45:21 +02:00
|
|
|
INSERT INTO POINT_TBL(f1) VALUES ('(1e-300,-1e-300)'); -- To underflow
|
|
|
|
INSERT INTO POINT_TBL(f1) VALUES ('(1e+300,Inf)'); -- To overflow
|
|
|
|
INSERT INTO POINT_TBL(f1) VALUES (' ( Nan , NaN ) ');
|
2010-11-23 21:27:50 +01:00
|
|
|
-- bad format points
|
2000-01-05 07:06:23 +01:00
|
|
|
INSERT INTO POINT_TBL(f1) VALUES ('asdfasdf');
|
2003-09-25 08:58:07 +02:00
|
|
|
ERROR: invalid input syntax for type point: "asdfasdf"
|
2008-09-01 22:42:46 +02:00
|
|
|
LINE 1: INSERT INTO POINT_TBL(f1) VALUES ('asdfasdf');
|
|
|
|
^
|
2000-01-05 07:06:23 +01:00
|
|
|
INSERT INTO POINT_TBL(f1) VALUES ('10.0,10.0');
|
|
|
|
INSERT INTO POINT_TBL(f1) VALUES ('(10.0 10.0)');
|
2003-09-25 08:58:07 +02:00
|
|
|
ERROR: invalid input syntax for type point: "(10.0 10.0)"
|
2008-09-01 22:42:46 +02:00
|
|
|
LINE 1: INSERT INTO POINT_TBL(f1) VALUES ('(10.0 10.0)');
|
|
|
|
^
|
2018-09-26 10:45:21 +02:00
|
|
|
INSERT INTO POINT_TBL(f1) VALUES ('(10.0, 10.0) x');
|
|
|
|
ERROR: invalid input syntax for type point: "(10.0, 10.0) x"
|
|
|
|
LINE 1: INSERT INTO POINT_TBL(f1) VALUES ('(10.0, 10.0) x');
|
|
|
|
^
|
2000-01-05 07:06:23 +01:00
|
|
|
INSERT INTO POINT_TBL(f1) VALUES ('(10.0,10.0');
|
2003-09-25 08:58:07 +02:00
|
|
|
ERROR: invalid input syntax for type point: "(10.0,10.0"
|
2008-09-01 22:42:46 +02:00
|
|
|
LINE 1: INSERT INTO POINT_TBL(f1) VALUES ('(10.0,10.0');
|
|
|
|
^
|
2018-09-26 10:45:21 +02:00
|
|
|
INSERT INTO POINT_TBL(f1) VALUES ('(10.0, 1e+500)'); -- Out of range
|
|
|
|
ERROR: "1e+500" is out of range for type double precision
|
|
|
|
LINE 1: INSERT INTO POINT_TBL(f1) VALUES ('(10.0, 1e+500)');
|
|
|
|
^
|
2005-04-07 03:51:41 +02:00
|
|
|
SELECT '' AS six, * FROM POINT_TBL;
|
2018-09-26 10:45:21 +02:00
|
|
|
six | f1
|
|
|
|
-----+-------------------
|
2000-01-05 07:06:23 +01:00
|
|
|
| (0,0)
|
|
|
|
| (-10,0)
|
|
|
|
| (-3,4)
|
|
|
|
| (5.1,34.5)
|
|
|
|
| (-5,-12)
|
2018-09-26 10:45:21 +02:00
|
|
|
| (1e-300,-1e-300)
|
|
|
|
| (1e+300,Infinity)
|
|
|
|
| (NaN,NaN)
|
2000-01-05 07:06:23 +01:00
|
|
|
| (10,10)
|
2018-09-26 10:45:21 +02:00
|
|
|
(9 rows)
|
1997-04-05 13:24:54 +02:00
|
|
|
|
2010-11-23 21:27:50 +01:00
|
|
|
-- left of
|
2000-01-05 07:06:23 +01:00
|
|
|
SELECT '' AS three, p.* FROM POINT_TBL p WHERE p.f1 << '(0.0, 0.0)';
|
|
|
|
three | f1
|
|
|
|
-------+----------
|
|
|
|
| (-10,0)
|
|
|
|
| (-3,4)
|
|
|
|
| (-5,-12)
|
1997-04-05 13:24:54 +02:00
|
|
|
(3 rows)
|
|
|
|
|
2010-11-23 21:27:50 +01:00
|
|
|
-- right of
|
2000-01-05 07:06:23 +01:00
|
|
|
SELECT '' AS three, p.* FROM POINT_TBL p WHERE '(0.0,0.0)' >> p.f1;
|
|
|
|
three | f1
|
|
|
|
-------+----------
|
|
|
|
| (-10,0)
|
|
|
|
| (-3,4)
|
|
|
|
| (-5,-12)
|
1997-04-05 13:24:54 +02:00
|
|
|
(3 rows)
|
|
|
|
|
2010-11-23 21:27:50 +01:00
|
|
|
-- above
|
2000-01-05 07:06:23 +01:00
|
|
|
SELECT '' AS one, p.* FROM POINT_TBL p WHERE '(0.0,0.0)' >^ p.f1;
|
|
|
|
one | f1
|
|
|
|
-----+----------
|
|
|
|
| (-5,-12)
|
1997-04-05 13:24:54 +02:00
|
|
|
(1 row)
|
|
|
|
|
2010-11-23 21:27:50 +01:00
|
|
|
-- below
|
2000-01-05 07:06:23 +01:00
|
|
|
SELECT '' AS one, p.* FROM POINT_TBL p WHERE p.f1 <^ '(0.0, 0.0)';
|
|
|
|
one | f1
|
|
|
|
-----+----------
|
|
|
|
| (-5,-12)
|
1997-04-05 13:24:54 +02:00
|
|
|
(1 row)
|
|
|
|
|
2010-11-23 21:27:50 +01:00
|
|
|
-- equal
|
2000-01-05 07:06:23 +01:00
|
|
|
SELECT '' AS one, p.* FROM POINT_TBL p WHERE p.f1 ~= '(5.1, 34.5)';
|
|
|
|
one | f1
|
|
|
|
-----+------------
|
|
|
|
| (5.1,34.5)
|
1997-04-05 13:24:54 +02:00
|
|
|
(1 row)
|
|
|
|
|
2010-11-23 21:27:50 +01:00
|
|
|
-- point in box
|
2000-01-05 07:06:23 +01:00
|
|
|
SELECT '' AS three, p.* FROM POINT_TBL p
|
2006-09-10 02:29:35 +02:00
|
|
|
WHERE p.f1 <@ box '(0,0,100,100)';
|
2000-01-05 07:06:23 +01:00
|
|
|
three | f1
|
|
|
|
-------+------------
|
|
|
|
| (0,0)
|
|
|
|
| (5.1,34.5)
|
|
|
|
| (10,10)
|
1997-04-27 04:58:38 +02:00
|
|
|
(3 rows)
|
1997-04-05 13:24:54 +02:00
|
|
|
|
2010-01-14 17:31:09 +01:00
|
|
|
SELECT '' AS three, p.* FROM POINT_TBL p
|
|
|
|
WHERE box '(0,0,100,100)' @> p.f1;
|
|
|
|
three | f1
|
|
|
|
-------+------------
|
|
|
|
| (0,0)
|
|
|
|
| (5.1,34.5)
|
|
|
|
| (10,10)
|
|
|
|
(3 rows)
|
|
|
|
|
2000-01-05 07:06:23 +01:00
|
|
|
SELECT '' AS three, p.* FROM POINT_TBL p
|
2006-09-10 02:29:35 +02:00
|
|
|
WHERE not p.f1 <@ box '(0,0,100,100)';
|
2018-09-26 10:45:21 +02:00
|
|
|
three | f1
|
|
|
|
-------+-------------------
|
2000-01-05 07:06:23 +01:00
|
|
|
| (-10,0)
|
|
|
|
| (-3,4)
|
|
|
|
| (-5,-12)
|
2018-09-26 10:45:21 +02:00
|
|
|
| (1e-300,-1e-300)
|
|
|
|
| (1e+300,Infinity)
|
|
|
|
| (NaN,NaN)
|
|
|
|
(6 rows)
|
1997-04-05 13:24:54 +02:00
|
|
|
|
2000-01-05 07:06:23 +01:00
|
|
|
SELECT '' AS two, p.* FROM POINT_TBL p
|
2006-09-10 02:29:35 +02:00
|
|
|
WHERE p.f1 <@ path '[(0,0),(-10,0),(-10,10)]';
|
2018-09-26 10:45:21 +02:00
|
|
|
two | f1
|
|
|
|
-----+------------------
|
2000-01-05 07:06:23 +01:00
|
|
|
| (0,0)
|
|
|
|
| (-10,0)
|
2018-09-26 10:45:21 +02:00
|
|
|
| (1e-300,-1e-300)
|
|
|
|
(3 rows)
|
1997-04-05 13:24:54 +02:00
|
|
|
|
2010-01-14 17:31:09 +01:00
|
|
|
SELECT '' AS three, p.* FROM POINT_TBL p
|
|
|
|
WHERE not box '(0,0,100,100)' @> p.f1;
|
2018-09-26 10:45:21 +02:00
|
|
|
three | f1
|
|
|
|
-------+-------------------
|
2010-01-14 17:31:09 +01:00
|
|
|
| (-10,0)
|
|
|
|
| (-3,4)
|
|
|
|
| (-5,-12)
|
2018-09-26 10:45:21 +02:00
|
|
|
| (1e-300,-1e-300)
|
|
|
|
| (1e+300,Infinity)
|
|
|
|
| (NaN,NaN)
|
|
|
|
(6 rows)
|
2010-01-14 17:31:09 +01:00
|
|
|
|
2000-01-05 07:06:23 +01:00
|
|
|
SELECT '' AS six, p.f1, p.f1 <-> point '(0,0)' AS dist
|
1997-05-11 17:42:09 +02:00
|
|
|
FROM POINT_TBL p
|
|
|
|
ORDER BY dist;
|
2018-09-26 10:45:21 +02:00
|
|
|
six | f1 | dist
|
|
|
|
-----+-------------------+----------------------
|
|
|
|
| (0,0) | 0
|
|
|
|
| (1e-300,-1e-300) | 1.4142135623731e-300
|
|
|
|
| (-3,4) | 5
|
|
|
|
| (-10,0) | 10
|
|
|
|
| (-5,-12) | 13
|
|
|
|
| (10,10) | 14.142135623731
|
|
|
|
| (5.1,34.5) | 34.8749193547455
|
|
|
|
| (1e+300,Infinity) | Infinity
|
|
|
|
| (NaN,NaN) | NaN
|
|
|
|
(9 rows)
|
1997-04-05 13:24:54 +02:00
|
|
|
|
2000-01-05 07:06:23 +01:00
|
|
|
SELECT '' AS thirtysix, p1.f1 AS point1, p2.f1 AS point2, p1.f1 <-> p2.f1 AS dist
|
1997-05-11 17:42:09 +02:00
|
|
|
FROM POINT_TBL p1, POINT_TBL p2
|
2007-01-09 03:14:16 +01:00
|
|
|
ORDER BY dist, p1.f1[0], p2.f1[0];
|
2018-09-26 10:45:21 +02:00
|
|
|
thirtysix | point1 | point2 | dist
|
|
|
|
-----------+-------------------+-------------------+----------------------
|
|
|
|
| (-10,0) | (-10,0) | 0
|
|
|
|
| (-5,-12) | (-5,-12) | 0
|
|
|
|
| (-3,4) | (-3,4) | 0
|
|
|
|
| (0,0) | (0,0) | 0
|
|
|
|
| (1e-300,-1e-300) | (1e-300,-1e-300) | 0
|
|
|
|
| (5.1,34.5) | (5.1,34.5) | 0
|
|
|
|
| (10,10) | (10,10) | 0
|
|
|
|
| (0,0) | (1e-300,-1e-300) | 1.4142135623731e-300
|
|
|
|
| (1e-300,-1e-300) | (0,0) | 1.4142135623731e-300
|
|
|
|
| (-3,4) | (0,0) | 5
|
|
|
|
| (-3,4) | (1e-300,-1e-300) | 5
|
|
|
|
| (0,0) | (-3,4) | 5
|
|
|
|
| (1e-300,-1e-300) | (-3,4) | 5
|
|
|
|
| (-10,0) | (-3,4) | 8.06225774829855
|
|
|
|
| (-3,4) | (-10,0) | 8.06225774829855
|
|
|
|
| (-10,0) | (0,0) | 10
|
|
|
|
| (-10,0) | (1e-300,-1e-300) | 10
|
|
|
|
| (0,0) | (-10,0) | 10
|
|
|
|
| (1e-300,-1e-300) | (-10,0) | 10
|
|
|
|
| (-10,0) | (-5,-12) | 13
|
|
|
|
| (-5,-12) | (-10,0) | 13
|
|
|
|
| (-5,-12) | (0,0) | 13
|
|
|
|
| (-5,-12) | (1e-300,-1e-300) | 13
|
|
|
|
| (0,0) | (-5,-12) | 13
|
|
|
|
| (1e-300,-1e-300) | (-5,-12) | 13
|
|
|
|
| (0,0) | (10,10) | 14.142135623731
|
|
|
|
| (1e-300,-1e-300) | (10,10) | 14.142135623731
|
|
|
|
| (10,10) | (0,0) | 14.142135623731
|
|
|
|
| (10,10) | (1e-300,-1e-300) | 14.142135623731
|
|
|
|
| (-3,4) | (10,10) | 14.3178210632764
|
|
|
|
| (10,10) | (-3,4) | 14.3178210632764
|
|
|
|
| (-5,-12) | (-3,4) | 16.1245154965971
|
|
|
|
| (-3,4) | (-5,-12) | 16.1245154965971
|
|
|
|
| (-10,0) | (10,10) | 22.3606797749979
|
|
|
|
| (10,10) | (-10,0) | 22.3606797749979
|
|
|
|
| (5.1,34.5) | (10,10) | 24.9851956166046
|
|
|
|
| (10,10) | (5.1,34.5) | 24.9851956166046
|
|
|
|
| (-5,-12) | (10,10) | 26.6270539113887
|
|
|
|
| (10,10) | (-5,-12) | 26.6270539113887
|
|
|
|
| (-3,4) | (5.1,34.5) | 31.5572495632937
|
|
|
|
| (5.1,34.5) | (-3,4) | 31.5572495632937
|
|
|
|
| (0,0) | (5.1,34.5) | 34.8749193547455
|
|
|
|
| (1e-300,-1e-300) | (5.1,34.5) | 34.8749193547455
|
|
|
|
| (5.1,34.5) | (0,0) | 34.8749193547455
|
|
|
|
| (5.1,34.5) | (1e-300,-1e-300) | 34.8749193547455
|
|
|
|
| (-10,0) | (5.1,34.5) | 37.6597928831267
|
|
|
|
| (5.1,34.5) | (-10,0) | 37.6597928831267
|
|
|
|
| (-5,-12) | (5.1,34.5) | 47.5842410888311
|
|
|
|
| (5.1,34.5) | (-5,-12) | 47.5842410888311
|
|
|
|
| (-10,0) | (1e+300,Infinity) | Infinity
|
|
|
|
| (-5,-12) | (1e+300,Infinity) | Infinity
|
|
|
|
| (-3,4) | (1e+300,Infinity) | Infinity
|
|
|
|
| (0,0) | (1e+300,Infinity) | Infinity
|
|
|
|
| (1e-300,-1e-300) | (1e+300,Infinity) | Infinity
|
|
|
|
| (5.1,34.5) | (1e+300,Infinity) | Infinity
|
|
|
|
| (10,10) | (1e+300,Infinity) | Infinity
|
|
|
|
| (1e+300,Infinity) | (-10,0) | Infinity
|
|
|
|
| (1e+300,Infinity) | (-5,-12) | Infinity
|
|
|
|
| (1e+300,Infinity) | (-3,4) | Infinity
|
|
|
|
| (1e+300,Infinity) | (0,0) | Infinity
|
|
|
|
| (1e+300,Infinity) | (1e-300,-1e-300) | Infinity
|
|
|
|
| (1e+300,Infinity) | (5.1,34.5) | Infinity
|
|
|
|
| (1e+300,Infinity) | (10,10) | Infinity
|
|
|
|
| (-10,0) | (NaN,NaN) | NaN
|
|
|
|
| (-5,-12) | (NaN,NaN) | NaN
|
|
|
|
| (-3,4) | (NaN,NaN) | NaN
|
|
|
|
| (0,0) | (NaN,NaN) | NaN
|
|
|
|
| (1e-300,-1e-300) | (NaN,NaN) | NaN
|
|
|
|
| (5.1,34.5) | (NaN,NaN) | NaN
|
|
|
|
| (10,10) | (NaN,NaN) | NaN
|
|
|
|
| (1e+300,Infinity) | (1e+300,Infinity) | NaN
|
|
|
|
| (1e+300,Infinity) | (NaN,NaN) | NaN
|
|
|
|
| (NaN,NaN) | (-10,0) | NaN
|
|
|
|
| (NaN,NaN) | (-5,-12) | NaN
|
|
|
|
| (NaN,NaN) | (-3,4) | NaN
|
|
|
|
| (NaN,NaN) | (0,0) | NaN
|
|
|
|
| (NaN,NaN) | (1e-300,-1e-300) | NaN
|
|
|
|
| (NaN,NaN) | (5.1,34.5) | NaN
|
|
|
|
| (NaN,NaN) | (10,10) | NaN
|
|
|
|
| (NaN,NaN) | (1e+300,Infinity) | NaN
|
|
|
|
| (NaN,NaN) | (NaN,NaN) | NaN
|
|
|
|
(81 rows)
|
1997-04-05 13:24:54 +02:00
|
|
|
|
2000-01-05 07:06:23 +01:00
|
|
|
SELECT '' AS thirty, p1.f1 AS point1, p2.f1 AS point2
|
1997-04-05 13:24:54 +02:00
|
|
|
FROM POINT_TBL p1, POINT_TBL p2
|
1997-07-29 18:24:04 +02:00
|
|
|
WHERE (p1.f1 <-> p2.f1) > 3;
|
2018-09-26 10:45:21 +02:00
|
|
|
thirty | point1 | point2
|
|
|
|
--------+-------------------+-------------------
|
|
|
|
| (0,0) | (-10,0)
|
|
|
|
| (0,0) | (-3,4)
|
|
|
|
| (0,0) | (5.1,34.5)
|
|
|
|
| (0,0) | (-5,-12)
|
|
|
|
| (0,0) | (1e+300,Infinity)
|
|
|
|
| (0,0) | (NaN,NaN)
|
|
|
|
| (0,0) | (10,10)
|
|
|
|
| (-10,0) | (0,0)
|
|
|
|
| (-10,0) | (-3,4)
|
|
|
|
| (-10,0) | (5.1,34.5)
|
|
|
|
| (-10,0) | (-5,-12)
|
|
|
|
| (-10,0) | (1e-300,-1e-300)
|
|
|
|
| (-10,0) | (1e+300,Infinity)
|
|
|
|
| (-10,0) | (NaN,NaN)
|
|
|
|
| (-10,0) | (10,10)
|
|
|
|
| (-3,4) | (0,0)
|
|
|
|
| (-3,4) | (-10,0)
|
|
|
|
| (-3,4) | (5.1,34.5)
|
|
|
|
| (-3,4) | (-5,-12)
|
|
|
|
| (-3,4) | (1e-300,-1e-300)
|
|
|
|
| (-3,4) | (1e+300,Infinity)
|
|
|
|
| (-3,4) | (NaN,NaN)
|
|
|
|
| (-3,4) | (10,10)
|
|
|
|
| (5.1,34.5) | (0,0)
|
|
|
|
| (5.1,34.5) | (-10,0)
|
|
|
|
| (5.1,34.5) | (-3,4)
|
|
|
|
| (5.1,34.5) | (-5,-12)
|
|
|
|
| (5.1,34.5) | (1e-300,-1e-300)
|
|
|
|
| (5.1,34.5) | (1e+300,Infinity)
|
|
|
|
| (5.1,34.5) | (NaN,NaN)
|
|
|
|
| (5.1,34.5) | (10,10)
|
|
|
|
| (-5,-12) | (0,0)
|
|
|
|
| (-5,-12) | (-10,0)
|
|
|
|
| (-5,-12) | (-3,4)
|
|
|
|
| (-5,-12) | (5.1,34.5)
|
|
|
|
| (-5,-12) | (1e-300,-1e-300)
|
|
|
|
| (-5,-12) | (1e+300,Infinity)
|
|
|
|
| (-5,-12) | (NaN,NaN)
|
|
|
|
| (-5,-12) | (10,10)
|
|
|
|
| (1e-300,-1e-300) | (-10,0)
|
|
|
|
| (1e-300,-1e-300) | (-3,4)
|
|
|
|
| (1e-300,-1e-300) | (5.1,34.5)
|
|
|
|
| (1e-300,-1e-300) | (-5,-12)
|
|
|
|
| (1e-300,-1e-300) | (1e+300,Infinity)
|
|
|
|
| (1e-300,-1e-300) | (NaN,NaN)
|
|
|
|
| (1e-300,-1e-300) | (10,10)
|
|
|
|
| (1e+300,Infinity) | (0,0)
|
|
|
|
| (1e+300,Infinity) | (-10,0)
|
|
|
|
| (1e+300,Infinity) | (-3,4)
|
|
|
|
| (1e+300,Infinity) | (5.1,34.5)
|
|
|
|
| (1e+300,Infinity) | (-5,-12)
|
|
|
|
| (1e+300,Infinity) | (1e-300,-1e-300)
|
|
|
|
| (1e+300,Infinity) | (1e+300,Infinity)
|
|
|
|
| (1e+300,Infinity) | (NaN,NaN)
|
|
|
|
| (1e+300,Infinity) | (10,10)
|
|
|
|
| (NaN,NaN) | (0,0)
|
|
|
|
| (NaN,NaN) | (-10,0)
|
|
|
|
| (NaN,NaN) | (-3,4)
|
|
|
|
| (NaN,NaN) | (5.1,34.5)
|
|
|
|
| (NaN,NaN) | (-5,-12)
|
|
|
|
| (NaN,NaN) | (1e-300,-1e-300)
|
|
|
|
| (NaN,NaN) | (1e+300,Infinity)
|
|
|
|
| (NaN,NaN) | (NaN,NaN)
|
|
|
|
| (NaN,NaN) | (10,10)
|
|
|
|
| (10,10) | (0,0)
|
|
|
|
| (10,10) | (-10,0)
|
|
|
|
| (10,10) | (-3,4)
|
|
|
|
| (10,10) | (5.1,34.5)
|
|
|
|
| (10,10) | (-5,-12)
|
|
|
|
| (10,10) | (1e-300,-1e-300)
|
|
|
|
| (10,10) | (1e+300,Infinity)
|
|
|
|
| (10,10) | (NaN,NaN)
|
|
|
|
(72 rows)
|
1997-04-05 13:24:54 +02:00
|
|
|
|
2000-01-05 07:06:23 +01:00
|
|
|
-- put distance result into output to allow sorting with GEQ optimizer - tgl 97/05/10
|
|
|
|
SELECT '' AS fifteen, p1.f1 AS point1, p2.f1 AS point2, (p1.f1 <-> p2.f1) AS distance
|
1997-04-05 13:24:54 +02:00
|
|
|
FROM POINT_TBL p1, POINT_TBL p2
|
1997-07-29 18:24:04 +02:00
|
|
|
WHERE (p1.f1 <-> p2.f1) > 3 and p1.f1 << p2.f1
|
2007-01-09 03:14:16 +01:00
|
|
|
ORDER BY distance, p1.f1[0], p2.f1[0];
|
2018-09-26 10:45:21 +02:00
|
|
|
fifteen | point1 | point2 | distance
|
|
|
|
---------+------------------+-------------------+------------------
|
|
|
|
| (-3,4) | (0,0) | 5
|
|
|
|
| (-3,4) | (1e-300,-1e-300) | 5
|
|
|
|
| (-10,0) | (-3,4) | 8.06225774829855
|
|
|
|
| (-10,0) | (0,0) | 10
|
|
|
|
| (-10,0) | (1e-300,-1e-300) | 10
|
|
|
|
| (-10,0) | (-5,-12) | 13
|
|
|
|
| (-5,-12) | (0,0) | 13
|
|
|
|
| (-5,-12) | (1e-300,-1e-300) | 13
|
|
|
|
| (0,0) | (10,10) | 14.142135623731
|
|
|
|
| (1e-300,-1e-300) | (10,10) | 14.142135623731
|
|
|
|
| (-3,4) | (10,10) | 14.3178210632764
|
|
|
|
| (-5,-12) | (-3,4) | 16.1245154965971
|
|
|
|
| (-10,0) | (10,10) | 22.3606797749979
|
|
|
|
| (5.1,34.5) | (10,10) | 24.9851956166046
|
|
|
|
| (-5,-12) | (10,10) | 26.6270539113887
|
|
|
|
| (-3,4) | (5.1,34.5) | 31.5572495632937
|
|
|
|
| (0,0) | (5.1,34.5) | 34.8749193547455
|
|
|
|
| (1e-300,-1e-300) | (5.1,34.5) | 34.8749193547455
|
|
|
|
| (-10,0) | (5.1,34.5) | 37.6597928831267
|
|
|
|
| (-5,-12) | (5.1,34.5) | 47.5842410888311
|
|
|
|
| (-10,0) | (1e+300,Infinity) | Infinity
|
|
|
|
| (-5,-12) | (1e+300,Infinity) | Infinity
|
|
|
|
| (-3,4) | (1e+300,Infinity) | Infinity
|
|
|
|
| (0,0) | (1e+300,Infinity) | Infinity
|
|
|
|
| (1e-300,-1e-300) | (1e+300,Infinity) | Infinity
|
|
|
|
| (5.1,34.5) | (1e+300,Infinity) | Infinity
|
|
|
|
| (10,10) | (1e+300,Infinity) | Infinity
|
|
|
|
(27 rows)
|
1997-04-05 13:24:54 +02:00
|
|
|
|
2000-01-05 07:06:23 +01:00
|
|
|
-- put distance result into output to allow sorting with GEQ optimizer - tgl 97/05/10
|
|
|
|
SELECT '' AS three, p1.f1 AS point1, p2.f1 AS point2, (p1.f1 <-> p2.f1) AS distance
|
2010-11-23 21:27:50 +01:00
|
|
|
FROM POINT_TBL p1, POINT_TBL p2
|
1997-07-29 18:24:04 +02:00
|
|
|
WHERE (p1.f1 <-> p2.f1) > 3 and p1.f1 << p2.f1 and p1.f1 >^ p2.f1
|
1997-05-11 17:42:09 +02:00
|
|
|
ORDER BY distance;
|
2018-09-26 10:45:21 +02:00
|
|
|
three | point1 | point2 | distance
|
|
|
|
-------+------------+------------------+------------------
|
|
|
|
| (-3,4) | (0,0) | 5
|
|
|
|
| (-3,4) | (1e-300,-1e-300) | 5
|
|
|
|
| (-10,0) | (-5,-12) | 13
|
|
|
|
| (5.1,34.5) | (10,10) | 24.9851956166046
|
|
|
|
(4 rows)
|
1997-04-05 13:24:54 +02:00
|
|
|
|
Fix gist_box_same and gist_point_consistent to handle fuzziness correctly.
While there's considerable doubt that we want fuzzy behavior in the
geometric operators at all (let alone as currently implemented), nobody is
stepping forward to redesign that stuff. In the meantime it behooves us
to make sure that index searches agree with the behavior of the underlying
operators. This patch fixes two problems in this area.
First, gist_box_same was using fuzzy equality, but it really needs to use
exact equality to prevent not-quite-identical upper index keys from being
treated as identical, which for example would prevent an existing upper
key from being extended by an amount less than epsilon. This would result
in inconsistent indexes. (The next release notes will need to recommend
that users reindex GiST indexes on boxes, polygons, circles, and points,
since all four opclasses use gist_box_same.)
Second, gist_point_consistent used exact comparisons for upper-page
comparisons in ~= searches, when it needs to use fuzzy comparisons to
ensure it finds all matches; and it used fuzzy comparisons for point <@ box
searches, when it needs to use exact comparisons because that's what the
<@ operator (rather inconsistently) does.
The added regression test cases illustrate all three misbehaviors.
Back-patch to all active branches. (8.4 did not have GiST point_ops,
but it still seems prudent to apply the gist_box_same patch to it.)
Alexander Korotkov, reviewed by Noah Misch
2013-02-09 00:03:17 +01:00
|
|
|
-- Test that GiST indexes provide same behavior as sequential scan
|
|
|
|
CREATE TEMP TABLE point_gist_tbl(f1 point);
|
|
|
|
INSERT INTO point_gist_tbl SELECT '(0,0)' FROM generate_series(0,1000);
|
|
|
|
CREATE INDEX point_gist_tbl_index ON point_gist_tbl USING gist (f1);
|
|
|
|
INSERT INTO point_gist_tbl VALUES ('(0.0000009,0.0000009)');
|
|
|
|
SET enable_seqscan TO true;
|
|
|
|
SET enable_indexscan TO false;
|
|
|
|
SET enable_bitmapscan TO false;
|
|
|
|
SELECT COUNT(*) FROM point_gist_tbl WHERE f1 ~= '(0.0000009,0.0000009)'::point;
|
|
|
|
count
|
|
|
|
-------
|
|
|
|
1002
|
|
|
|
(1 row)
|
|
|
|
|
|
|
|
SELECT COUNT(*) FROM point_gist_tbl WHERE f1 <@ '(0.0000009,0.0000009),(0.0000009,0.0000009)'::box;
|
|
|
|
count
|
|
|
|
-------
|
|
|
|
1
|
|
|
|
(1 row)
|
|
|
|
|
|
|
|
SELECT COUNT(*) FROM point_gist_tbl WHERE f1 ~= '(0.0000018,0.0000018)'::point;
|
|
|
|
count
|
|
|
|
-------
|
|
|
|
1
|
|
|
|
(1 row)
|
|
|
|
|
|
|
|
SET enable_seqscan TO false;
|
|
|
|
SET enable_indexscan TO true;
|
|
|
|
SET enable_bitmapscan TO true;
|
|
|
|
SELECT COUNT(*) FROM point_gist_tbl WHERE f1 ~= '(0.0000009,0.0000009)'::point;
|
|
|
|
count
|
|
|
|
-------
|
|
|
|
1002
|
|
|
|
(1 row)
|
|
|
|
|
|
|
|
SELECT COUNT(*) FROM point_gist_tbl WHERE f1 <@ '(0.0000009,0.0000009),(0.0000009,0.0000009)'::box;
|
|
|
|
count
|
|
|
|
-------
|
|
|
|
1
|
|
|
|
(1 row)
|
|
|
|
|
|
|
|
SELECT COUNT(*) FROM point_gist_tbl WHERE f1 ~= '(0.0000018,0.0000018)'::point;
|
|
|
|
count
|
|
|
|
-------
|
|
|
|
1
|
|
|
|
(1 row)
|
|
|
|
|
|
|
|
RESET enable_seqscan;
|
|
|
|
RESET enable_indexscan;
|
|
|
|
RESET enable_bitmapscan;
|