Changeset 8033

Timestamp:

Oct 10, 2010, 5:15:16 PM (14 years ago)

Author:

habili

Message:

Uses polygon_overlap method to find triangles that overlap the inlet.

File:

• trunk/anuga_core/source/anuga/structures/inlet.py (modified) (4 diffs)

trunk/anuga_core/source/anuga/structures/inlet.py

@@ -1 @@
-from anuga.geometry.polygon import inside_polygon, is_inside_polygon
+from anuga.geometry.polygon import inside_polygon, is_inside_polygon, polygon_overlap
 from anuga.config import velocity_protection, g
 import math
@@ -27 +27 @@
         bounding_polygon = self.domain_bounding_polygon
         domain_centroids = self.domain.get_centroid_coordinates(absolute=True)
-
+        vertex_coordinates = self.domain.get_vertex_coordinates()
 
         # Check that polygon lies within the mesh.
@@ -40 +40 @@
         assert is_inside_polygon(point, bounding_polygon), msg
 
-
-        self.triangle_indices = inside_polygon(domain_centroids, self.polygon, verbose=self.verbose)
+        self.triangle_indices = polygon_overlap(vertex_coordinates, self.polygon)
+        #self.triangle_indices_cen = inside_polygon(domain_centroids, self.polygon, verbose=self.verbose)
 
         if len(self.triangle_indices) == 0:
@@ -49 +49 @@
 
         self.enquiry_index = self.domain.get_triangle_containing_point(self.enquiry_pt)
-