Changeset 5193

Timestamp:

Apr 8, 2008, 10:03:37 AM (17 years ago)

Author:

duncan

Message:

Creating an example importing an ungenerate file

File:

• anuga_core/source/anuga/pmesh/test_mesh_interface.py (modified) (5 diffs)

anuga_core/source/anuga/pmesh/test_mesh_interface.py

@@ -43 +43 @@
                          change_points_geo_ref(inner2_polygon_absolute)
         
-        interior_regions = [(inner1_polygon, 5),(inner2_polygon, 0.2)]
+        interior_regions = [(inner1_polygon, 5),(inner2_polygon, 10)]
         m = create_mesh_from_regions(polygon,
                                      boundary_tags,
@@ -137 +137 @@
                          change_points_geo_ref(inner2_polygon_absolute)
         
-        interior_regions = [(inner1_polygon, 5),(inner2_polygon, 0.2)]
+        interior_regions = [(inner1_polygon, 5),(inner2_polygon, 10)]
 
         interior_holes = None
@@ -265 +265 @@
                          change_points_geo_ref(inner2_polygon_absolute)
         
-        interior_regions = [(inner1_polygon, 5),(inner2_polygon, 0.2)]
+        interior_regions = [(inner1_polygon, 5),(inner2_polygon, 10)]
         m = create_mesh_from_regions(polygon,
                                      boundary_tags,
@@ -321 +321 @@
                          change_points_geo_ref(inner2_polygon_absolute)
         
-        interior_regions = [(inner1_polygon, 5),(inner2_polygon, 0.2)]
+        interior_regions = [(inner1_polygon, 5),(inner2_polygon, 10)]
         m = create_mesh_from_regions(polygon,
                                      boundary_tags,
@@ -742 +742 @@
         fileName = "badmesh.tsh"
         #m.export_mesh_file(fileName)
-        
+  
+        
+    def concept_create_mesh_from_regions_with_ungenerate(self):
+        x=0
+        y=0
+        mesh_geo = geo_reference=Geo_reference(56,x,y)
+        
+        # These are the absolute values
+        polygon_absolute = [[0,0],[100,0],[100,100],[0,100]]
+        
+        x_p = -10
+        y_p = -40
+        geo_ref_poly = Geo_reference(56, x_p, y_p)
+        polygon = geo_ref_poly.change_points_geo_ref(polygon_absolute)
+
+        boundary_tags = {'walls':[0,1],'bom':[2]}
+        
+        inner1_polygon_absolute = [[10,10],[20,10],[20,20],[10,20]]
+        inner1_polygon = geo_ref_poly. \
+                         change_points_geo_ref(inner1_polygon_absolute)
+
+        inner2_polygon_absolute = [[30,30],[40,30],[40,40],[30,40]]
+        inner2_polygon = geo_ref_poly. \
+                         change_points_geo_ref(inner2_polygon_absolute)
+        
+        max_area = 10000000
+        interior_regions = [(inner1_polygon, 5),(inner2_polygon, 10)]
+        m = create_mesh_from_regions(polygon,
+                                     boundary_tags,
+                                     max_area,
+                                     interior_regions=interior_regions,
+                                     poly_geo_reference=geo_ref_poly,
+                                     mesh_geo_reference=mesh_geo)
+                    
+        m.export_mesh_file('a_test_mesh_iknterface.tsh')                 
+        
+        fileName = tempfile.mktemp(".txt")
+        file = open(fileName,"w")
+        file.write("         1       ??      ??\n\
+       90.0       90.0\n\
+       81.0       90.0\n\
+       81.0       81.0\n\
+       90.0       81.0\n\
+       90.0       90.0\n\
+END\n\
+         2      ?? ??\n\
+       10.0       80.0\n\
+       10.0       90.0\n\
+       20.0       90.0\n\
+       10.0       80.0\n\
+END\n\
+END\n")
+        file.close() 
+        
+        m.import_ungenerate_file(fileName, tag='wall')
+        os.remove(fileName)
+        m.generate_mesh(maximum_triangle_area=max_area,verbose=False)
+        m.export_mesh_file('b_test_mesh_iknterface.tsh')
 #-------------------------------------------------------------
 if __name__ == "__main__":
     suite = unittest.makeSuite(TestCase,'test')
-    #suite = unittest.makeSuite(TestCase,'test_create_mesh_from_regions4')
+    #suite = unittest.makeSuite(TestCase,'test_create_mesh_from_regions_with_ungenerate')
     runner = unittest.TextTestRunner() #verbosity=2)
     runner.run(suite)