[postgis-tickets] [PostGIS] #4726: Garden hang on ST_VoronoiLines

[PostGIS]

#4726: Garden hang on ST_VoronoiLines
---------------------+---------------------------
 Reporter:  robe     |      Owner:  pramsey
     Type:  defect   |     Status:  new
 Priority:  blocker  |  Milestone:  PostGIS 3.1.0
Component:  postgis  |    Version:  master
 Keywords:           |
---------------------+---------------------------
 This hangs indefinitely on my PostgreSQL 13 Geos 3.9


 {{{
 POSTGIS="3.1.0alpha2 3.1.0alpha1-168-gceca2a80d" [EXTENSION] PGSQL="130"
 GEOS="3.9.0-CAPI-1.14.0" SFCGAL="1.3.8" PROJ="6.2.1" GDAL="GDAL 3.0.2,
 released 2019/10/28" LIBXML="2.9.9" LIBJSON="0.12" LIBPROTOBUF="1.2.1"
 WAGYU="0.5.0 (Internal)" RASTER PostgreSQL 13beta2 on x86_64-w64-mingw32,
 compiled by gcc.exe (x86_64-posix-seh-rev0, Built by MinGW-W64 project)
 8.1.0, 64-bit
 }}}



 {{{
 SELECT ST_VoronoiLines(foo1.the_geom, 20.1, foo2.the_geom) As result,
 ST_AsText(foo1.the_geom) As ref1_geom, ST_AsText(foo2.the_geom) As
 ref2_geom FROM ((SELECT ST_Collect(s.the_geom) As the_geom FROM (SELECT
 ST_SetSRID(ST_Point(i,j),4326) As the_geom FROM generate_series(-10,50,15)
 As i CROSS JOIN generate_series(40,70, 15) j ) As s)) As foo1 CROSS JOIN
 ((SELECT ST_SetSRID(ST_Point(i,j),4326) As the_geom FROM (SELECT
 a*1.11111111 FROM generate_series(-10,50,2) As a) As i(i) CROSS JOIN
 generate_series(40,70, 5) j ORDER BY i,j )) As foo2 LIMIT 10;
 }}}


 So does this reduced one:


 {{{
 SELECT ST_VoronoiLines('MULTIPOINT(-10 40,5 40,20 40,35 40,50 40,-10 55,5
 55,20 55,35 55,50 55,-10 70,5 70,20 70,35 70,50 70)'::geometry, 20.1,
 'POINT(-11.1111111 40)'::geometry)
 }}}

 and it's none-cancellable.


 I tried it on


 {{{
 POSTGIS="2.5.3 r17699" [EXTENSION] PGSQL="110" GEOS="3.8.0-CAPI-1.13.1 "
 PROJ="Rel. 4.9.3, 15 August 2016" GDAL="GDAL 2.4.4, released 2020/01/08"
 LIBXML="2.7.8" LIBJSON="0.12" LIBPROTOBUF="1.2.1" RASTER
 }}}


 and doesn't hang.

  and as text rep is


 {{{
 MULTILINESTRING((-2.5 47.5,-2.5 62.5),(-2.5 62.5,5 70),(5 70,20 55),(20
 55,12.5 47.5),(12.5 47.5,-2.5 47.5),(-2.5 47.5,-70 47.5),(-70 62.5,-2.5
 62.5),(-2.5 47.5,-2.5 -20),(12.5 47.5,12.5 -20),(5 130,5 70),(20 55,27.5
 55),(27.5 55,27.5 -20),(27.5 130,27.5 55),(42.5 130,42.5 55),(42.5 55,27.5
 55),(42.5 55,42.5 -20),(110 55,42.5 55))
 }}}



 I'll try the latest GEOS 3.9 head next

-- 
Ticket URL: <https://trac.osgeo.org/postgis/ticket/4726>
PostGIS <http://trac.osgeo.org/postgis/>
The PostGIS Trac is used for bug, enhancement & task tracking, a user and developer wiki, and a view into the subversion code repository of PostGIS project.