[geos-commits] r3341 - in trunk/tests/unit: . capi

[svn_geos at osgeo.org]

Author: strk
Date: 2011-05-10 04:11:54 -0700 (Tue, 10 May 2011)
New Revision: 3341

Added:
   trunk/tests/unit/capi/GEOSIntersectsTest.cpp
Modified:
   trunk/tests/unit/Makefile.am
Log:
Add test for 357 (doesn't succeed at failing)

Modified: trunk/tests/unit/Makefile.am
===================================================================
--- trunk/tests/unit/Makefile.am	2011-05-10 09:51:31 UTC (rev 3340)
+++ trunk/tests/unit/Makefile.am	2011-05-10 11:11:54 UTC (rev 3341)
@@ -98,6 +98,7 @@
 	capi/GEOSGeomFromWKBTest.cpp \
 	capi/GEOSGeomToWKTTest.cpp \
 	capi/GEOSContainsTest.cpp \
+	capi/GEOSIntersectsTest.cpp \
 	capi/GEOSWithinTest.cpp \
 	capi/GEOSSimplifyTest.cpp \
 	capi/GEOSPreparedGeometryTest.cpp \

Added: trunk/tests/unit/capi/GEOSIntersectsTest.cpp
===================================================================
--- trunk/tests/unit/capi/GEOSIntersectsTest.cpp	                        (rev 0)
+++ trunk/tests/unit/capi/GEOSIntersectsTest.cpp	2011-05-10 11:11:54 UTC (rev 3341)
@@ -0,0 +1,146 @@
+// $Id$
+// 
+// Test Suite for C-API GEOSIntersects
+
+#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_capigeosintersects_data
+    {
+        GEOSGeometry* geom1_;
+        GEOSGeometry* geom2_;
+
+        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_capigeosintersects_data()
+            : geom1_(0), geom2_(0)
+        {
+            initGEOS(notice, notice);
+        }       
+
+        ~test_capigeosintersects_data()
+        {
+            GEOSGeom_destroy(geom1_);
+            GEOSGeom_destroy(geom2_);
+            geom1_ = 0;
+            geom2_ = 0;
+            finishGEOS();
+        }
+
+    };
+
+    typedef test_group<test_capigeosintersects_data> group;
+    typedef group::object object;
+
+    group test_capigeosintersects_group("capi::GEOSIntersects");
+
+    //
+    // Test Cases
+    //
+
+    template<>
+    template<>
+    void object::test<1>()
+    {
+        geom1_ = GEOSGeomFromWKT("POLYGON EMPTY");
+        geom2_ = GEOSGeomFromWKT("POLYGON EMPTY");
+
+        ensure( 0 != geom1_ );
+        ensure( 0 != geom2_ );
+
+        char const r1 = GEOSIntersects(geom1_, geom2_);
+
+        ensure_equals(r1, 0);
+
+        char const r2 = GEOSIntersects(geom2_, geom1_);
+
+        ensure_equals(r2, 0);
+    }
+
+    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_ );
+
+        char const r1 = GEOSIntersects(geom1_, geom2_);
+
+        ensure_equals(int(r1), 1);
+
+        char const r2 = GEOSIntersects(geom2_, geom1_);
+
+        ensure_equals(int(r2), 1);
+    }
+    
+    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_ );
+
+        char const r1 = GEOSIntersects(geom1_, geom2_);
+
+        ensure_equals(int(r1), 1);
+        
+        char const r2 = GEOSIntersects(geom2_, geom1_);
+
+        ensure_equals(int(r2), 1);
+    }
+
+    // This is a tef for bug #357
+    template<>
+    template<>
+    void object::test<4>()
+    {
+        GEOSCoordSequence* cs = GEOSCoordSeq_create(5, 2);
+
+        double nan = -0.0/0;
+        GEOSCoordSeq_setX(cs, 0, 1); GEOSCoordSeq_setY(cs, 0, 1);
+        for (unsigned int i=1; i<4; ++i) {
+            GEOSCoordSeq_setX(cs, i, nan);
+            GEOSCoordSeq_setY(cs, i, nan);
+        }
+        GEOSCoordSeq_setX(cs, 4, 1); GEOSCoordSeq_setY(cs, 4, 1);
+
+	geom1_ = GEOSGeom_createPolygon(
+		GEOSGeom_createLinearRing(cs),
+		NULL, 0);
+
+        char const r1 = GEOSIntersects(geom1_, geom1_);
+
+        ensure_equals(int(r1), 1);
+        
+    }
+ 
+} // namespace tut
+