[postgis-users] Geomunion crash

[Pierre Racine]

Hi,

The following query makes geomunion(geometry set) to crash. I know there is 
two identical shapes but when I use an other parameter than "3" in buffer() 
(like 3.001), everything works fine. It also crash with "4". If I add 
"DISTINCT" in set SELECT subquery, it works fine. My understanding is that I 
should not have to...

DROP TABLE bug1;
CREATE TABLE bug1 (lat float4,long float4);
INSERT INTO bug1 (lat,long) VALUES (45.85,-1.0667);
INSERT INTO bug1 (lat,long) VALUES (47.7,-3.333);
INSERT INTO bug1 (lat,long) VALUES (47.7,-3.333);
SELECT geomunion(geomset.geom)
   FROM (SELECT buffer(geomfromtext('POINT(' || long || ' ' || lat || ')', 
4269),3) as geom FROM bug1) as geomset;

The result is:

NOTICE:  TopologyException: no outgoing dirEdge found (-3.85329,44.7512)

ERROR:  GEOS union() threw an error!

Pierre Racine