--- operator testing (testing is on the BOUNDING BOX (2d), not the actual geometries)

-- overlap or left of

select '77','POINT(1 1)'::GEOMETRY &< 'POINT(1 1)'::GEOMETRY as bool;
select '78','POINT(1 1)'::GEOMETRY &< 'POINT(2 1)'::GEOMETRY as bool;
select '79','POINT(2 1)'::GEOMETRY &< 'POINT(1 1)'::GEOMETRY as bool;

-- strictly left

select '80','POINT(1 1)'::GEOMETRY << 'POINT(1 1)'::GEOMETRY as bool;
select '81','POINT(1 1)'::GEOMETRY << 'POINT(2 1)'::GEOMETRY as bool;
select '82','POINT(2 1)'::GEOMETRY << 'POINT(1 1)'::GEOMETRY as bool;

-- overlap or right

select '83','POINT(1 1)'::GEOMETRY &> 'POINT(1 1)'::GEOMETRY as bool;
select '84','POINT(1 1)'::GEOMETRY &> 'POINT(2 1)'::GEOMETRY as bool;
select '85','POINT(2 1)'::GEOMETRY &> 'POINT(1 1)'::GEOMETRY as bool;

-- strictly right

select '86','POINT(1 1)'::GEOMETRY >> 'POINT(1 1)'::GEOMETRY as bool;
select '87','POINT(1 1)'::GEOMETRY >> 'POINT(2 1)'::GEOMETRY as bool;
select '88','POINT(2 1)'::GEOMETRY >> 'POINT(1 1)'::GEOMETRY as bool;

-- overlap

select '89','POINT(1 1)'::GEOMETRY && 'POINT(1 1)'::GEOMETRY as bool;
select '90','POINT(1 1)'::GEOMETRY && 'POINT(2 2)'::GEOMETRY as bool;
select '91','MULTIPOINT(0 0, 1 1)'::GEOMETRY && 'MULTIPOINT(1 1, 2 2)'::GEOMETRY as bool;
select '92','MULTIPOINT(0 0, 1 1)'::GEOMETRY && 'MULTIPOINT(1.0001 1, 2 2)'::GEOMETRY as bool;
select '93','MULTIPOINT(0 0, 1 1)'::GEOMETRY && 'MULTIPOINT(1 1.0001, 2 2)'::GEOMETRY as bool;
select '94','MULTIPOINT(0 0, 1 1)'::GEOMETRY && 'MULTIPOINT(1 0, 2 2)'::GEOMETRY as bool;
select '95','MULTIPOINT(0 0, 1 1)'::GEOMETRY && 'MULTIPOINT(1.0001 0, 2 2)'::GEOMETRY as bool;

select '96','MULTIPOINT(0 0, 1 1)'::GEOMETRY && 'MULTIPOINT(0 1, 1 2)'::GEOMETRY as bool;
select '97','MULTIPOINT(0 0, 1 1)'::GEOMETRY && 'MULTIPOINT(0 1.0001, 1 2)'::GEOMETRY as bool;

--- contains

select '98','MULTIPOINT(0 0, 10 10)'::GEOMETRY ~ 'MULTIPOINT(5 5, 7 7)'::GEOMETRY as bool;
select '99','MULTIPOINT(5 5, 7 7)'::GEOMETRY ~ 'MULTIPOINT(0 0, 10 10)'::GEOMETRY as bool;
select '100','MULTIPOINT(0 0, 7 7)'::GEOMETRY ~ 'MULTIPOINT(0 0, 10 10)'::GEOMETRY as bool;
select '101','MULTIPOINT(-0.0001 0, 7 7)'::GEOMETRY ~ 'MULTIPOINT(0 0, 10 10)'::GEOMETRY as bool;

--- contained by

select '102','MULTIPOINT(0 0, 10 10)'::GEOMETRY @ 'MULTIPOINT(5 5, 7 7)'::GEOMETRY as bool;
select '103','MULTIPOINT(5 5, 7 7)'::GEOMETRY @ 'MULTIPOINT(0 0, 10 10)'::GEOMETRY as bool;
select '104','MULTIPOINT(0 0, 7 7)'::GEOMETRY @ 'MULTIPOINT(0 0, 10 10)'::GEOMETRY as bool;
select '105','MULTIPOINT(-0.0001 0, 7 7)'::GEOMETRY @ 'MULTIPOINT(0 0, 10 10)'::GEOMETRY as bool;

-- overlap or below  &<|

select 'ovbl1',ST_MakeEnvelope(2,2,4,4) &<| ST_MakeEnvelope(2,2,4,4); --t
select 'ovbl2',ST_MakeEnvelope(2,1,4,4) &<| ST_MakeEnvelope(2,2,4,4); --t
select 'ovbl3',ST_MakeEnvelope(2,1,4,5) &<| ST_MakeEnvelope(2,2,4,4); --f
select 'ovbl4',ST_MakeEnvelope(2,0,4,1) &<| ST_MakeEnvelope(2,2,4,4); --t

-- strictly below    <<|

select 'bl1',ST_MakeEnvelope(2,1,4,4) <<| ST_MakeEnvelope(2,2,4,4); --f
select 'bl2',ST_MakeEnvelope(2,1,4,2) <<| ST_MakeEnvelope(2,2,4,4); --f
select 'bl2',ST_MakeEnvelope(2,0,4,1) <<| ST_MakeEnvelope(2,2,4,4); --t

-- overlap or above  |&>

select 'ovab1',ST_MakeEnvelope(2,2,4,4) |&> ST_MakeEnvelope(2,2,4,4); --t
select 'ovab2',ST_MakeEnvelope(2,2,4,5) |&> ST_MakeEnvelope(2,2,4,4); --t
select 'ovab3',ST_MakeEnvelope(2,1,4,5) |&> ST_MakeEnvelope(2,2,4,4); --f
select 'ovab4',ST_MakeEnvelope(2,5,4,8) |&> ST_MakeEnvelope(2,2,4,4); --t

-- strictly above    |>>

select 'ab1',ST_MakeEnvelope(2,2,4,8) |>> ST_MakeEnvelope(2,2,4,4); --f
select 'ab2',ST_MakeEnvelope(2,4,4,8) |>> ST_MakeEnvelope(2,2,4,4); --f
select 'ab3',ST_MakeEnvelope(2,5,4,8) |>> ST_MakeEnvelope(2,2,4,4); --t

-- same as           =

select 'eq1',ST_MakeEnvelope(2,2,4,4) = ST_MakeEnvelope(2,2,4,4); -- t
select 'eq2',ST_MakeEnvelope(2,4,4,8) = 'LINESTRING(2 4,4 8)'::geometry; -- f
select 'eq3',ST_MakePoint(0,0) = ST_MakePoint(1,0); -- f

-- inequality    !=  <>

select 'neq1',ST_MakeEnvelope(2,2,4,4) != ST_MakeEnvelope(2,2,4,4); -- f
select 'neq2',ST_MakeEnvelope(2,4,4,8) <> 'LINESTRING(2 4,4 8)'::geometry; -- t

-- box centroid distance  <->

select 'cd1', 'LINESTRING(0 0,0 10)'::geometry <->
              'LINESTRING(4 0,4 10)'::geometry; -- 4

-- box distance           <#>

select 'bd1', 'LINESTRING(0 0,0 10,10 10)'::geometry <#>
              'LINESTRING(6 2,6 8)'::geometry; -- 0
select 'bd2', 'LINESTRING(0 0,0 10,10 10)'::geometry <#>
              'LINESTRING(11 0,19 10)'::geometry; -- 1

-- nd overlap             &&&

select 'ndov1', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry &&&
                'POINT(3 3 3 5)'::geometry; -- f
select 'ndov2', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry &&&
                'POINT(3 3 5 3)'::geometry; -- f
select 'ndov3', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry &&&
                'POINT(3 5 3 3)'::geometry; -- f
select 'ndov4', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry &&&
                'POINT(5 3 3 3)'::geometry; -- f
select 'ndov5', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry &&&
                'POINT(3 3 3 3)'::geometry; -- t
select 'ndov6', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry &&&
                'POINT(2 4 2 4)'::geometry; -- t
select 'ndov7', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry &&&
                'POINT(4 2 4 2)'::geometry; -- t

-- &&& with mixed dimensions

WITH v(i,g) AS ( VALUES
 (1,'POINT(0 0)'::geometry),			-- true for {1,2,3,4,5,6,7,8,9}
 (2,'POINTZ(0 0 1)'),					-- true for {1,2,4,6,7}
 (3,'POINTZ(0 0 0)'),					-- true for {1,3,4,5,7}
 (4,'POINTM(0 0 1)'),					-- true for {1,2,3,4,5,8}
 (5,'POINTZM(0 0 0 1)'),				-- true for {1,3,4,5}
 (6,'POINTZM(0 0 1 0)'),					-- true for {1,2,6}
 (7,'LINESTRINGM(-1 0 2,1 0 3)'),		-- true for {1,2,3,7,8,9}
 (8,'LINESTRINGZ(-1 0 2,1 0 3)'),		-- true for {1,4,7,8,9}
 (9,'LINESTRINGZM(-1 0 2 2,1 0 3 3)')	-- true for {1,7,8,9}
 )
SELECT 'ndovm', v1.i, array_agg(v2.i) FROM v v1, v v2 WHERE v1.g &&& v2.g
group by v1.i
ORDER BY 2;

-- nd contains             ~~

select 'ndcont1', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry ~~
                'POINT(3 3 3 5)'::geometry; -- f
select 'ndcont2', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry ~~
                'POINT(3 3 5 3)'::geometry; -- f
select 'ndcont3', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry ~~
                'POINT(3 5 3 3)'::geometry; -- f
select 'ndcont4', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry ~~
                'POINT(5 3 3 3)'::geometry; -- f
select 'ndcont5', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry ~~
                'POINT(3 3 3 3)'::geometry; -- t
select 'ndcont6', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry ~~
                'POINT(2 4 2 4)'::geometry; -- t
select 'ndcont7', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry ~~
                'POINT(4 2 4 2)'::geometry; -- t

-- ~~ with mixed dimensions

WITH v(i,g) AS ( VALUES
 (1,'POINT(0 0)'::geometry),
 (2,'POINTZ(0 0 1)'),
 (3,'POINTZ(0 0 0)'),
 (4,'POINTM(0 0 1)'),
 (5,'POINTZM(0 0 0 1)'),
 (6,'POINTZM(0 0 1 0)'),
 (7,'LINESTRINGM(-1 0 2,1 0 3)'),
 (8,'LINESTRINGZ(-1 0 2,1 0 3)'),
 (9,'LINESTRINGZM(-1 0 2 2,1 0 3 3)')
 )
SELECT 'ndcontm', v1.i, array_agg(v2.i) FROM v v1, v v2 WHERE v1.g ~~ v2.g
group by v1.i
ORDER BY 2;

-- nd within             @@

select 'ndwithin1', 'POINT(3 3 3 5)'::geometry @@
                'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry; -- f
select 'ndwithin2', 'POINT(3 3 5 3)'::geometry @@
                'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry; -- f
select 'ndwithin3', 'POINT(3 5 3 3)'::geometry @@
                'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry; -- f
select 'ndwithin4', 'POINT(5 3 3 3)'::geometry @@
                'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry; -- f
select 'ndwithin5', 'POINT(3 3 3 3)'::geometry @@
                'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry; -- t
select 'ndwithin6', 'POINT(2 4 2 4)'::geometry @@
                'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry; -- t
select 'ndwithin7', 'POINT(4 2 4 2)'::geometry @@
                'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry; -- t

-- @@ with mixed dimensions

WITH v(i,g) AS ( VALUES
 (1,'POINT(0 0)'::geometry),
 (2,'POINTZ(0 0 1)'),
 (3,'POINTZ(0 0 0)'),
 (4,'POINTM(0 0 1)'),
 (5,'POINTZM(0 0 0 1)'),
 (6,'POINTZM(0 0 1 0)'),
 (7,'LINESTRINGM(-1 0 2,1 0 3)'),
 (8,'LINESTRINGZ(-1 0 2,1 0 3)'),
 (9,'LINESTRINGZM(-1 0 2 2,1 0 3 3)')
 )
SELECT 'ndwithinm', v1.i, array_agg(v2.i) FROM v v1, v v2 WHERE v1.g @@ v2.g
group by v1.i
ORDER BY 2;

-- nd same             ~~=

select 'ndsame1', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry ~~=
                'POINT(3 3 3 3)'::geometry; -- f
select 'ndsame2', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry ~~=
                'LINESTRING(2 2 2 2, 4 4 4 5)'::geometry; -- f
select 'ndsame3', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry ~~=
                'LINESTRING(2 2 2 2, 4 4 5 4)'::geometry; -- f
select 'ndsame4', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry ~~=
                'LINESTRING(2 2 2 2, 4 5 4 4)'::geometry; -- f
select 'ndsame5', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry ~~=
                'LINESTRING(2 2 2 2, 5 4 4 4)'::geometry; -- f
select 'ndsame6', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry ~~=
                'LINESTRING(4 4 4 4, 2 2 2 2)'::geometry; -- t
select 'ndsame7', 'LINESTRING(2 2 2 2, 4 4 4 4)'::geometry ~~=
                'LINESTRING(2 2 2 4, 2 2 4 2, 2 4 2 2, 4 2 2 2)'::geometry; -- t

-- ~~= with mixed dimensions

WITH v(i,g) AS ( VALUES
 (1,'POINT(0 0)'::geometry),
 (2,'POINTZ(0 0 1)'),
 (3,'POINTZ(0 0 0)'),
 (4,'POINTM(0 0 1)'),
 (5,'POINTZM(0 0 0 1)'),
 (6,'POINTZM(0 0 1 0)'),
 (7,'LINESTRINGM(-1 0 2,1 0 3)'),
 (8,'LINESTRINGZ(-1 0 2,1 0 3)'),
 (9,'LINESTRINGZM(-1 0 2 2,1 0 3 3)')
 )
SELECT 'ndsamem', v1.i, array_agg(v2.i) FROM v v1, v v2 WHERE v1.g ~~= v2.g
group by v1.i
ORDER BY 2;

-- GROUP BY on empty
SELECT '#3777', ST_AsText(geom), count(*)
FROM (VALUES
    ('POINT(0 0)'::geometry),
    ('POINT(0 0)'::geometry),
    ('POINT(0 0)'::geometry),
    ('POINT(0 1)'::geometry),
    ('LINESTRING(0 0,0 1)'::geometry),
    ('GEOMETRYCOLLECTION EMPTY'::geometry),
    ('POINT EMPTY'::geometry)
) AS f(geom)
GROUP BY geom ORDER BY 2;

SELECT '#3777.1', ST_AsText(geom), count(*)
FROM (VALUES
    ('POINT(0 0)'::geometry),
    ('POINT(0 0)'::geometry),
    ('POINT EMPTY'::geometry),
    ('POINT(0 0)'::geometry),
    ('POINT(0 1)'::geometry),
    ('LINESTRING(0 0,0 1)'::geometry),
    ('GEOMETRYCOLLECTION EMPTY'::geometry)
) AS f(geom)
GROUP BY geom ORDER BY 2;
