Changeset 8037 for trunk/anuga_core/source/anuga/geometry/polygon.py

Timestamp:

Oct 15, 2010, 5:16:59 PM (13 years ago)

Author:

habili

Message:

added polyline_overlap routines

File:

• trunk/anuga_core/source/anuga/geometry/polygon.py (modified) (2 diffs)

trunk/anuga_core/source/anuga/geometry/polygon.py

@@ -257 +257 @@
     return indices[count:]    
 
+def polyline_overlap(triangles, polyline, verbose=False):
+    """Determine if a polygon and triangle overlap
+
+    """
+    polyline = ensure_numeric(polyline)
+    triangles = ensure_numeric(triangles)
+    
+    M = triangles.shape[0]/3  # Number of triangles
+
+    indices = num.zeros(M, num.int)
+
+    count = _polyline_overlap(polyline, triangles, indices)
+
+    if verbose:
+        log.critical('Found %d triangles (out of %d) that polygon' % (count, M))
+
+    return indices[:count]
+    
+def not_polyline_overlap(triangles, polyline, verbose=False):
+    """Determine if a polyline and triangle overlap
+
+    """
+    polyline = ensure_numeric(polyline)
+    triangles = ensure_numeric(triangles)
+    
+    M = triangles.shape[0]/3  # Number of triangles
+
+    indices = num.zeros(M, num.int)
+
+    count = _polyline_overlap(polyline, triangles, indices)
+
+    if verbose:
+        log.critical('Found %d triangles (out of %d) that polygon' % (count, M))
+
+    return indices[count:]    
+    
+
 def is_inside_triangle(point, triangle, 
                        closed=True, 
@@ -1135 +1172 @@
     from polygon_ext import _interpolate_polyline    
     from polygon_ext import _polygon_overlap
+    from polygon_ext import _polyline_overlap
     from polygon_ext import _is_inside_triangle        
     #from polygon_ext import _intersection