Changeset 4579 for anuga_core/source/anuga/abstract_2d_finite_volumes/quantity.py

Timestamp:

Jul 3, 2007, 1:30:49 PM (17 years ago)

Author:

ole

Message:

Thoughts towards allowing 'polygon' to be specified in set_quantity.
This is far from done.

File:

• anuga_core/source/anuga/abstract_2d_finite_volumes/quantity.py (modified) (4 diffs)

anuga_core/source/anuga/abstract_2d_finite_volumes/quantity.py

@@ -179 +179 @@
                    alpha = None,
                    location = 'vertices',
+                   polygon = None,
                    indices = None,
                    verbose = False,
@@ -262 +263 @@
                   will be left undefined.
 
+
+        polygon: Restrict update of quantity to locations that fall
+                 inside polygon. Polygon works by selecting indices
+                 and calling set_values recursively.
+
+        indices: Restrict update of quantity to locations that are 
+                 identified by indices (e.g. node ids if location
+                 is 'vertices')        
+        
         verbose: True means that output to stdout is generated
 
@@ -278 +288 @@
         from types import FloatType, IntType, LongType, ListType, NoneType
         from Numeric import ArrayType
+
+
+        # Polygon situation
+        #if polygon is not None:
+        #    if indices is not None:
+        #        msg = 'Only one of polygon and indices can be specified'
+        #        raise Exception, msg
+        #
+        #    Need to get candidate points. I think we should
+        #    simplify this whole thing a bit. Do we really need location?
+        #    point_indices = inside_polygon(points, polygon)
+
+
 
         #General input checks
@@ -388 +411 @@
         """
 
+        # FIXME (Ole): Somehow indices refer to centroids
+        # rather than vertices as default. See unit test
+        # test_set_vertex_values_using_general_interface_with_subset(self):
+        
 
         if location == 'centroids':