Changeset 8541 for trunk/anuga_core/source/anuga/fit_interpolate/interpolate.py

Timestamp:

Aug 30, 2012, 10:35:14 AM (13 years ago)

Author:

steve

Message:

Dealing with holes

File:

• trunk/anuga_core/source/anuga/fit_interpolate/interpolate.py (modified) (2 diffs)

trunk/anuga_core/source/anuga/fit_interpolate/interpolate.py

@@ -371 +371 @@
 
 
+    def get_outside_poly_indices(self):
+        """
+        Return index of those data points outside (and in holes)
+        the mesh
+
+        Precondition: interpolation or interpolation_block has been called
+        """
+        return self.outside_poly_indices
+
+
     def _get_point_data_z(self, f, verbose=False):
         """
@@ -413 +423 @@
         if verbose: log.critical('Getting indices inside mesh boundary')
 
-        # Quick test against boundary, but will not deal with holes in the mesh 
+        # Quick test against boundary, but will not deal with holes in the mesh,
+        # that is done below
         inside_boundary_indices, outside_poly_indices = \
             in_and_outside_polygon(point_coordinates,