[postgis-tickets] [PostGIS] #5144: 3d intersects in postgis hangs on simple polyhedron

[PostGIS]

#5144: 3d intersects in postgis hangs on simple polyhedron
-------------------------+----------------------------
 Reporter:  lbartoletti  |      Owner:  robe
     Type:  defect       |     Status:  new
 Priority:  medium       |  Milestone:  PostGIS SFCGAL
Component:  sfcgal       |    Version:  3.2.x
 Keywords:               |
-------------------------+----------------------------
 > This postgis query hangs without returning,
 >
 >
 {{{
 SELECT ST_3DIntersects('POINT(0 0.5 0.5)', ST_MakeSolid('POLYHEDRALSURFACE
 Z (((2 0.5 0.5, 0 0 1, 0 0 0, 2 0.5 0.5)), ((2 0.5 0.5, 0 1 1, 0 0 1, 2
 0.5 0.5)), ((2 0.5 0.5, 0 1 0, 0 1 1, 2 0.5 0.5)), ((2 0.5 0.5, 0 0 0, 0 1
 0, 2 0.5 0.5)), ((0 0 0, 0 0 1, 0 1 1, 0 1 0, 0 0 0)))'))
 }}}
 >
 > The polyhedron is a pyramid with it's base being the y-z unit square and
 apex in (2, 0.5, 0.5). The problematic point (found by happenstance) seems
 to be precicely (0, 0.5, 0.5), located on the base of the rotated pyramid.
 >
 > Tested both on
 {{{
 postgis 3.2 with SFCGAL 1.3.8 (official postgis docker image)
 }}}
  as well as this environment:
 {{{
 POSTGIS="3.3.0dev 3.2.0-783-g042894569" \[EXTENSION\] PGSQL="110"
 GEOS="3.8.3-CAPI-1.13.4" SFCGAL="SFCGAL 1.4.1, CGAL 5.3, BOOST 1.78.0"
 PROJ="7.1.1" GDAL="GDAL 3.2.3, released 2021/04/27" LIBXML="2.9.9"
 LIBJSON="0.12" LIBPROTOBUF="1.2.1" WAGYU="0.5.0 (Internal)" RASTER ,
 PostgreSQL 11.2 on x86_64-pc-mingw64, compiled by gcc.exe (x86_64-posix-
 seh-rev0, Built by MinGW-W64 project) 8.1.0, 64-bit
 }}}

-- 
Ticket URL: <https://trac.osgeo.org/postgis/ticket/5144>
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.