Add RelateNG prepared test for matches P/A

Author: dr-jts

modules/core/src/test/java/org/locationtech/jts/operation/relateng/RelateNGTest.java

@@ -11,6 +11,8 @@
  */
 package org.locationtech.jts.operation.relateng;
 
+import org.locationtech.jts.geom.IntersectionMatrix;
+
 import junit.textui.TestRunner;
 
 public class RelateNGTest extends RelateNGTestCase {
@@ -668,5 +670,17 @@ public void testPreparedAA() {
     checkPrepared(a, b);
   }
 
-  
+  public void testPreparedPA() {
+    String a = "POINT (5 5)";
+    String b = "POLYGON ((1 9, 9 9, 9 1, 1 1, 1 9))";
+    checkPrepared(a, b);
+    checkPrepared(b, a);
+    
+    //-- see https://github.com/libgeos/geos/issues/1275 (not a bug, but a good test to have)
+    String pattern = "T*****FF*";
+    String patternTrans = IntersectionMatrix.transpose(pattern);  // T*F**F***
+    checkPreparedMatches(a, b, pattern);
+    checkPreparedMatches(b, a, patternTrans); //
+  }
+
 }

modules/core/src/test/java/org/locationtech/jts/operation/relateng/RelateNGTestCase.java

@@ -113,6 +113,16 @@ void checkPrepared(String wkta, String wktb)
                             RelateNG.relate(a, b).toString());
   }
 
+  void checkPreparedMatches(String wkta, String wktb, String pattern) {
+    Geometry a = read(wkta);
+    Geometry b = read(wktb);
+    RelateNG prep_a = RelateNG.prepare(a);
+
+    assertEquals("matches " + pattern,  prep_a.evaluate(b, RelatePredicate.matches(pattern) ),
+              RelateNG.relate(a, b, RelatePredicate.matches(pattern) )
+        );
+  }
+
   TopologyPredicate trace(TopologyPredicate pred) {
     if (! isTrace)
       return pred;