Changeset 3193 for inundation/pmesh/test_mesh_interface.py

Timestamp:

Jun 21, 2006, 3:33:49 PM (18 years ago)

Author:

duncan

Message:

bug fixing

File:

• inundation/pmesh/test_mesh_interface.py (modified) (5 diffs)

inundation/pmesh/test_mesh_interface.py

@@ -3 +3 @@
 import tempfile
 import unittest
+import os
+from mesh import importMeshFromFile
 from mesh_interface import *
 from load_mesh.loadASCII import *
@@ -221 +223 @@
                         'FAILED!')
 
+    def test_create_mesh_from_regions4(self):
+
+        file_name = tempfile.mktemp(".tsh")
+        
+        # These are the absolute values
+        density_outer = 1000
+        min_outer = 0 
+        max_outer = 1000
+        polygon_outer = [[min_outer,min_outer],[max_outer,min_outer],
+                   [max_outer,max_outer],[min_outer,max_outer]]
+        
+        density_inner1 = 10000000
+        inner_buffer = 100
+        min_inner1 = min_outer + inner_buffer
+        max_inner1 = max_outer - inner_buffer
+        inner1_polygon = [[min_inner1,min_inner1],[max_inner1,min_inner1],
+                   [max_inner1,max_inner1],[min_inner1,max_inner1]]
+      
+        
+        boundary_tags = {'walls':[0,1],'bom':[2]}
+        
+        interior_regions = [(inner1_polygon, density_inner1)]
+        create_mesh_from_regions(polygon_outer
+                                     , boundary_tags
+                                     , density_outer
+                                     , interior_regions=interior_regions
+                                     ,filename=file_name
+                                     #,verbose=True
+                                     ,verbose=False
+                                     )
+        
+        m = importMeshFromFile(file_name)
+        
+        #print "file_name",file_name 
+        self.failUnless(len(m.meshTriangles) <= 900,
+                        'Test mesh interface failed!')
+        self.failUnless(len(m.meshTriangles) >= 200,
+                        'Test mesh interface failed!')
+        
+        create_mesh_from_regions(polygon_outer
+                                     , boundary_tags
+                                     , interior_regions=interior_regions
+                                     ,filename=file_name
+                                     #,verbose=True
+                                     ,verbose=False
+                                     )
+        
+        m = importMeshFromFile(file_name)
+        
+        #print "len(m.meshTriangles)",len(m.meshTriangles) 
+        self.failUnless(len(m.meshTriangles) <= 100,
+                        'Test mesh interface failed!')
+
+        os.remove(file_name)
+        
+    def test_create_mesh_from_regions5(self):
+
+        file_name = tempfile.mktemp(".tsh")
+        
+        # These are the absolute values
+        density_outer = 10000000 
+        min_outer = 0 
+        max_outer = 1000
+        polygon_outer = [[min_outer,min_outer],[max_outer,min_outer],
+                   [max_outer,max_outer],[min_outer,max_outer]]
+        
+        density_inner1 = 1000
+        inner_buffer = 100
+        min_inner1 = min_outer + inner_buffer
+        max_inner1 = max_outer - inner_buffer
+        inner1_polygon = [[min_inner1,min_inner1],[max_inner1,min_inner1],
+                   [max_inner1,max_inner1],[min_inner1,max_inner1]]
+      
+        
+        boundary_tags = {'walls':[0,1],'bom':[2]}
+        
+        interior_regions = [(inner1_polygon, density_inner1)]
+        create_mesh_from_regions(polygon_outer
+                                     , boundary_tags
+                                     , density_outer
+                                     , interior_regions=interior_regions
+                                     ,filename=file_name
+                                     #,verbose=True
+                                     ,verbose=False
+                                     )
+        
+        m = importMeshFromFile(file_name)
+        #print "file_name",file_name
+        #print "len(m.meshTriangles",len(m.meshTriangles) 
+        self.failUnless(len(m.meshTriangles) <= 2000, 
+                        'Test mesh interface failed!')
+ 
+        self.failUnless(len(m.meshTriangles) >= 900,
+                        'Test mesh interface failed!')
+
+        os.remove(file_name)
+        
     def test_create_mesh_from_regions_interior_regions(self):
         """Test that create_mesh_from_regions fails when an interior region is
@@ -254 +353 @@
                                          boundary_tags,
                                          10000000,
-                                         interior_regions=interior_regions)
+                                         interior_regions=interior_regions,verbose=False)
         except:
             pass
@@ -287 +386 @@
                                      boundary_tags,
                                      10000000,
-                                     interior_regions=interior_regions)
+                                     interior_regions=interior_regions,verbose=False)
 
 
@@ -337 +436 @@
                                      boundary_tags,
                                      10000000,
-                                     poly_geo_reference=geo_ref_poly)
+                                     poly_geo_reference=geo_ref_poly,verbose=False)