Changeset 8003

Timestamp:

Sep 4, 2010, 9:45:34 PM (14 years ago)

Author:

habili

Message:

fix for is_inside_polygon

File:

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

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

@@ -18 +18 @@
         self.verbose = verbose
 
-        # FIXME (SR) Using get_triangle_containing_point which needs to be sped up
-
         self.compute_indices()
         self.compute_area()
@@ -32 +30 @@
 
         # Check that polygon lies within the mesh.
-        for point in self.polygon + self.enquiry_pt:
+        for point in self.polygon:
                 msg = 'Point %s ' %  str(point)
                 msg += ' did not fall within the domain boundary.'
                 assert is_inside_polygon(point, bounding_polygon), msg
+                
+        point = self.enquiry_pt
+        msg = 'Enquiry Point %s ' %  str(point)
+        msg += ' did not fall within the domain boundary.'
+        assert is_inside_polygon(point, bounding_polygon), msg