From svn_geos at osgeo.org  Mon Mar  2 03:22:29 2015
From: svn_geos at osgeo.org (svn_geos at osgeo.org)
Date: Mon,  2 Mar 2015 03:22:29 -0800
Subject: [geos-commits] r4048 - in trunk: src/operation/overlay tests/unit
	tests/unit/capi
Message-ID: <20150302112230.05F763903DB@trac.osgeo.org>

Author: strk
Date: 2015-03-02 03:22:29 -0800 (Mon, 02 Mar 2015)
New Revision: 4048

Added:
   trunk/tests/unit/capi/GEOSIntersectionTest.cpp
Modified:
   trunk/src/operation/overlay/OverlayOp.cpp
   trunk/tests/unit/Makefile.am
Log:
Fix memory leak in extent-limited overlay operations

The leak was introduced by recent enhancement of the class to
support extent-limited operation. It was spotted by Mick Orridge
who also provided the testcase, thanks !

Fixes #719

Modified: trunk/src/operation/overlay/OverlayOp.cpp
===================================================================
--- trunk/src/operation/overlay/OverlayOp.cpp	2015-02-03 10:11:19 UTC (rev 4047)
+++ trunk/src/operation/overlay/OverlayOp.cpp	2015-03-02 11:22:29 UTC (rev 4048)
@@ -194,7 +194,10 @@
 {
   for(size_t i=0, n=edges->size(); i<n; ++i) {
     Edge *e = (*edges)[i];
-		if ( env && ! env->intersects(e->getEnvelope()) ) continue;
+		if ( env && ! env->intersects(e->getEnvelope()) ) {
+		  dupEdges.push_back(e); // or could it be deleted directly ?
+      continue;
+    }
 #if GEOS_DEBUG
 		cerr <<" "<< e->print() << endl;
 #endif 

Modified: trunk/tests/unit/Makefile.am
===================================================================
--- trunk/tests/unit/Makefile.am	2015-02-03 10:11:19 UTC (rev 4047)
+++ trunk/tests/unit/Makefile.am	2015-03-02 11:22:29 UTC (rev 4048)
@@ -127,6 +127,7 @@
 	capi/GEOSDistanceTest.cpp \
 	capi/GEOSInterruptTest.cpp \
 	capi/GEOSIntersectsTest.cpp \
+	capi/GEOSIntersectionTest.cpp \
 	capi/GEOSNearestPointsTest.cpp \
 	capi/GEOSWithinTest.cpp \
 	capi/GEOSSimplifyTest.cpp \

Added: trunk/tests/unit/capi/GEOSIntersectionTest.cpp
===================================================================
--- trunk/tests/unit/capi/GEOSIntersectionTest.cpp	                        (rev 0)
+++ trunk/tests/unit/capi/GEOSIntersectionTest.cpp	2015-03-02 11:22:29 UTC (rev 4048)
@@ -0,0 +1,142 @@
+// 
+// Test Suite for C-API GEOSintersection
+
+#include <tut.hpp>
+// geos
+#include <geos_c.h>
+// std
+#include <cstdarg>
+#include <cstdio>
+#include <cstdlib>
+
+namespace tut
+{
+	//
+	// Test Group
+	//
+
+	// Common data used in test cases.
+	struct test_capigeosintersection_data
+	{
+		GEOSWKTWriter* wktw_;
+		GEOSGeometry* geom1_;
+		GEOSGeometry* geom2_;
+		GEOSGeometry* geom3_;
+
+		static void notice(const char *fmt, ...)
+		{
+			std::fprintf(stdout, "NOTICE: ");
+
+			va_list ap;
+			va_start(ap, fmt);
+			std::vfprintf(stdout, fmt, ap);
+			va_end(ap);
+
+			std::fprintf(stdout, "\n");
+		}
+
+		test_capigeosintersection_data()
+			: geom1_(0), geom2_(0), geom3_(0)
+		{
+			initGEOS(notice, notice);
+			wktw_ = GEOSWKTWriter_create();
+			GEOSWKTWriter_setTrim(wktw_, 1);
+			GEOSWKTWriter_setOutputDimension(wktw_, 3);
+		}
+
+		std::string toWKT(GEOSGeometry* g)
+		{
+			char* wkt = GEOSWKTWriter_write(wktw_, g);
+			std::string ret(wkt);
+			GEOSFree(wkt);
+			return ret;
+		}
+
+		~test_capigeosintersection_data()
+		{
+			GEOSWKTWriter_destroy(wktw_);
+			GEOSGeom_destroy(geom1_);
+			GEOSGeom_destroy(geom2_);
+			GEOSGeom_destroy(geom3_);
+			geom1_ = 0;
+			geom2_ = 0;
+			geom3_ = 0;
+			finishGEOS();
+		}
+
+	};
+
+	typedef test_group<test_capigeosintersection_data> group;
+	typedef group::object object;
+
+	group test_capigeosintersection_group("capi::GEOSIntersection");
+
+	//
+	// Test Cases
+	//
+
+	template<>
+	template<>
+	void object::test<1>()
+	{
+		geom1_ = GEOSGeomFromWKT("POLYGON EMPTY");
+		geom2_ = GEOSGeomFromWKT("POLYGON EMPTY");
+
+		ensure(0 != geom1_);
+		ensure(0 != geom2_);
+
+		geom3_ = GEOSIntersection(geom1_, geom2_);
+		ensure(0 != geom3_);
+		ensure_equals(toWKT(geom3_), std::string("GEOMETRYCOLLECTION EMPTY"));
+	}
+
+	template<>
+	template<>
+	void object::test<2>()
+	{
+		geom1_ = GEOSGeomFromWKT("POLYGON((1 1,1 5,5 5,5 1,1 1))");
+		geom2_ = GEOSGeomFromWKT("POINT(2 2)");
+
+		ensure(0 != geom1_);
+		ensure(0 != geom2_);
+
+		geom3_ = GEOSIntersection(geom1_, geom2_);
+		ensure(0 != geom3_);
+		ensure_equals(toWKT(geom3_), std::string("POINT (2 2)"));
+	}
+
+	template<>
+	template<>
+	void object::test<3>()
+	{
+		geom1_ = GEOSGeomFromWKT("MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)))");
+		geom2_ = GEOSGeomFromWKT("POLYGON((-1 1,-1 2,2 2,2 1,-1 1))");
+
+		ensure(0 != geom1_);
+		ensure(0 != geom2_);
+
+		geom3_ = GEOSIntersection(geom1_, geom2_);
+
+		ensure(0 != geom3_);
+		ensure_equals(toWKT(geom3_), std::string("POLYGON ((0 1, 0 2, 2 2, 2 1, 0 1))"));
+	}
+
+  /* See http://trac.osgeo.org/geos/ticket/719 */
+	template<>
+	template<>
+	void object::test<4>()
+	{
+		geom1_ = GEOSGeomFromWKT("MULTIPOLYGON(((0 0,5 10,10 0,0 0),(1 1,1 2,2 2,2 1,1 1),(100 100,100 102,102 102,102 100,100 100)))");
+		geom2_ = GEOSGeomFromWKT("POLYGON((0 1,0 2,10 2,10 1,0 1))");
+
+		ensure(0 != geom1_);
+		ensure(0 != geom2_);
+
+		geom3_ = GEOSIntersection(geom1_, geom2_);
+
+		ensure(0 != geom3_);
+		ensure_equals(toWKT(geom3_), std::string("GEOMETRYCOLLECTION (LINESTRING (1 2, 2 2), LINESTRING (2 1, 1 1), POLYGON ((0.5 1, 1 2, 1 1, 0.5 1)), POLYGON ((9 2, 9.5 1, 2 1, 2 2, 9 2)))"));
+	}
+
+} // namespace tut
+