<p>geos 3.5.0</p><p>i'm using geos to generate voronoi diagrams . i want to know each site i input into the function and its corresponding polygon.</p><p>here is the code according to page <a href="http://geos.osgeo.org/doxygen/classgeos_1_1triangulate_1_1quadedge_1_1QuadEdgeSubdivision.html#a3ee3bb33600640553458ba405e83f485">http://geos.osgeo.org/doxygen/classgeos_1_1triangulate_1_1quadedge_1_1QuadEdgeSubdivision.html#a3ee3bb33600640553458ba405e83f485</a></p><p> std::auto_ptr<GeometryCollection> polys = subdiv->getVoronoiDiagram(geomFact);<br> for (std::size_t i = 0; i < polys->getNumGeometries(); ++i) {<br>  Coordinate *c = (Coordinate*)((Geometry*)(polys->getGeometryN(i)))->getUserData();<br>  const Polygon* p = dynamic_cast<const Polygon*>(polys->getGeometryN(i));</p><p>  //ensure(p);<br>  std::auto_ptr<CoordinateSequence> cs(<br>   p->getExteriorRing()->getCoordinates()<br>   );<br>  size_t from = cs->size();<br>  cs->removeRepeatedPoints();<br>  //ensure_equals(from, cs->size());<br> }</p><p>coordinate of cordinate* c does not equal to any of my input sites.</p><p>so i feel in QuadEdgeSubdivision.cpp line 586, there should be some code passing a coordinate * variable to function setUserData using operator new.  </p><p><br></p><div></div>