Changeset 8003


Ignore:
Timestamp:
Sep 4, 2010, 9:45:34 PM (14 years ago)
Author:
habili
Message:

fix for is_inside_polygon

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/anuga_core/source/anuga/structures/inlet.py

    r7984 r8003  
    1818        self.verbose = verbose
    1919
    20         # FIXME (SR) Using get_triangle_containing_point which needs to be sped up
    21 
    2220        self.compute_indices()
    2321        self.compute_area()
     
    3230
    3331        # Check that polygon lies within the mesh.
    34         for point in self.polygon + self.enquiry_pt:
     32        for point in self.polygon:
    3533                msg = 'Point %s ' %  str(point)
    3634                msg += ' did not fall within the domain boundary.'
    3735                assert is_inside_polygon(point, bounding_polygon), msg
     36               
     37        point = self.enquiry_pt
     38        msg = 'Enquiry Point %s ' %  str(point)
     39        msg += ' did not fall within the domain boundary.'
     40        assert is_inside_polygon(point, bounding_polygon), msg
    3841
    3942
Note: See TracChangeset for help on using the changeset viewer.