ERROR: GEOSTopologyPreserveSimplify: vector::reserve

[strk]

I got the above error trying to simplifying a pre-simplified dataset. Still trying to find the smallest dataset reproducing it. May be hard because it looks like an out of memory condition (but I'm not sure it is). I bet it's related to emptyness.

[strk]

Interesting enough the _same_ input passed to ST_AsEWKT give a parsing error: ERROR: parse error - invalid geometry

going to canonical output and back to canonical input fixes the ST_AsEWKT but not ST_SimplifyPreserveTopology

[strk]

Alright, input looks like:

 MULTIPOLYGON ( ((...)), EMPTY )

WKT parser refuses to accept it.

It looks like I can reproduce using this call:

select ST_SimplifyPreserveTopology(ST_collect('MULTIPOLYGON ( ((0 0, 10 0, 10 10, 0 10, 0 0)) )', 'POLYGON EMPTY'), 0);

[strk]

So this is a GEOS bug, see ​http://trac.osgeo.org/geos/ticket/519

Let's keep this open until the other side is fixed.

[strk]

GEOS bug fixed in version 3.3.3