Changeset 5570 for anuga_core/source/anuga/shallow_water/shallow_water_domain.py

Timestamp:

Jul 25, 2008, 1:35:45 PM (15 years ago)

Author:

ole

Message:

Added input tests regarding polygons and points

File:

• anuga_core/source/anuga/shallow_water/shallow_water_domain.py (modified) (4 diffs)

anuga_core/source/anuga/shallow_water/shallow_water_domain.py

@@ -108 +108 @@
 
 
-from anuga.utilities.polygon import inside_polygon, polygon_area        
+from anuga.utilities.polygon import inside_polygon, polygon_area, is_inside_polygon
 
 
@@ -1501 +1501 @@
             
                      
-        from math import pi
+        from math import pi, cos, sin
 
         self.domain = domain
@@ -1513 +1513 @@
                          # previous timestep in order to obtain rate
 
+                         
+        bounding_polygon = domain.get_boundary_polygon()
+
+
         # Update area if applicable
         self.exchange_area = None        
@@ -1521 +1525 @@
 
             self.exchange_area = radius**2*pi
+
+            # Check that circle center lies within the mesh.
+            msg = 'Center %s specified for forcing term did not' %(str(center))
+            msg += 'fall within the domain boundary.'
+            assert is_inside_polygon(center, bounding_polygon), msg
+
+            # Check that circle periphery lies within the mesh.
+            N = 100
+            periphery_points = []
+            for i in range(N):
+
+                theta = 2*pi*i/100
+                
+                x = center[0] + radius*cos(theta)
+                y = center[1] + radius*sin(theta)
+
+                periphery_points.append([x,y])
+
+
+            for point in periphery_points:
+                msg = 'Point %s on periphery for forcing term did not' %(str(point))
+                msg += ' fall within the domain boundary.'
+                assert is_inside_polygon(point, bounding_polygon), msg
+
         
         if polygon is not None:
             self.exchange_area = polygon_area(self.polygon)
+
+            # Check that polygon lies within the mesh.
+            for point in self.polygon:
+                msg = 'Point %s in polygon for forcing term did not' %(str(point))
+                msg += 'fall within the domain boundary.'
+                assert is_inside_polygon(point, bounding_polygon), msg
+        
+
+
+