- Timestamp:
- Nov 24, 2010, 4:55:55 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/anuga_1_2_1/anuga_core/source/anuga/pmesh/test_mesh.py
r7873 r8083 1757 1757 self.failUnless(segs[2].tag=='bom', 1758 1758 'FAILED!') 1759 self.failUnless(segs[3].tag==' ',1759 self.failUnless(segs[3].tag=='interior', 1760 1760 'FAILED!') 1761 1761 verts = m.getUserVertices() … … 1775 1775 self.failUnless(point_x == new_point_x, ' failed') 1776 1776 self.failUnless(point_y == new_point_y, ' failed') 1777 1778 1779 1780 1781 def test_add_hole_from_polygon_none_tag(self): 1782 x=-500 1783 y=-1000 1784 m=Mesh(geo_reference=Geo_reference(56,x,y)) 1785 1786 # These are the absolute values 1787 polygon_absolute = [[0,0],[1,0],[1,1],[0,1]] 1788 1789 x_p = -10 1790 y_p = -40 1791 geo_ref_poly = Geo_reference(56, x_p, y_p) 1792 polygon = geo_ref_poly.change_points_geo_ref(polygon_absolute) 1793 1794 poly_point = m.add_hole_from_polygon(polygon, 1795 None, 1796 geo_reference=geo_ref_poly) 1797 # poly_point values are relative to the mesh geo-ref 1798 # make them absolute 1799 #print "poly_point.x+x",poly_point.x+x 1800 #print "poly_point.y+y",poly_point.y+y 1801 #print "polygon_absolute", polygon_absolute 1802 self.failUnless(is_inside_polygon([poly_point.x+x,poly_point.y+y], 1803 polygon_absolute, closed = False), 1804 'FAILED!') 1805 1806 self.failUnless(len(m.holes)==1, 1807 'FAILED!') 1808 segs = m.getUserSegments() 1809 self.failUnless(len(segs)==4, 1810 'FAILED!') 1811 self.failUnless(len(m.userVertices)==4, 1812 'FAILED!') 1813 1814 self.failUnless(segs[0].tag=='interior', 1815 'FAILED!') 1816 self.failUnless(segs[1].tag=='interior', 1817 'FAILED!') 1818 1819 self.failUnless(segs[2].tag=='interior', 1820 'FAILED!') 1821 self.failUnless(segs[3].tag=='interior', 1822 'FAILED!') 1823 verts = m.getUserVertices() 1824 #print "User verts",verts 1825 #print 'polygon',polygon 1826 #vert values are relative 1827 for point,new_point in map(None,polygon,verts): 1828 point_x = point[0] + geo_ref_poly.get_xllcorner() 1829 new_point_x = new_point.x + m.geo_reference.get_xllcorner() 1830 point_y = point[1] + geo_ref_poly.get_yllcorner() 1831 #print "new_point.y",new_point.y 1832 #print "m.geo_ref.get_yllcorner()",m.geo_reference.get_yllcorner() 1833 new_point_y = new_point.y + m.geo_reference.get_yllcorner() 1834 #print "point_y",point_y 1835 #print "new_point_y",new_point_y 1836 1837 self.failUnless(point_x == new_point_x, ' failed') 1838 self.failUnless(point_y == new_point_y, ' failed') 1777 1839 1778 1840 def test_add_circle(self):
Note: See TracChangeset
for help on using the changeset viewer.