Ignore:
Timestamp:
Jan 11, 2006, 5:31:49 PM (18 years ago)
Author:
duncan
Message:

in pmesh.mesh, can now add a region specified as a polygon. (This isn't fully implemented yet)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • inundation/utilities/test_polygon.py

    r2062 r2200  
    379379            assert not inside_polygon(point, ex_poly), '%s' %str(point)                       
    380380
    381 
     381    def test_point_in_polygon(self):
     382        polygon = [[0,0], [1,0], [1,1], [0,1]]
     383        point = point_in_polygon(polygon)
     384        assert inside_polygon(point, polygon)
     385
     386        #this get's into a vicious loop
     387        polygon = [[1e32,1e54], [1,0], [1,1], [0,1]]
     388        point = point_in_polygon(polygon)
     389        assert inside_polygon(point, polygon)
     390
     391
     392        polygon = [[0,0], [1,0], [1,1], [1e8,1e8]]
     393        point = point_in_polygon(polygon)
     394        assert inside_polygon(point, polygon)
     395       
     396        polygon = [[1e32,1e54], [-1e32,1e54],[1e32,-1e54]]
     397        point = point_in_polygon(polygon)
     398        assert inside_polygon(point, polygon)
     399       
     400        polygon = [[1e32,1e54], [1,0], [0,1]]
     401        point = point_in_polygon(polygon)
     402        assert inside_polygon(point, polygon)
    382403
    383404#-------------------------------------------------------------
Note: See TracChangeset for help on using the changeset viewer.