Changeset 7682

Timestamp:

Apr 12, 2010, 7:00:43 PM (15 years ago)

Author:

James Hudson

Message:

[Ticket 344] - produce asc when there is a hole in the model. Bug fix: holes no longer cause an exception, and are returned as outside points.

File:

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

anuga_core/source/anuga/fit_interpolate/interpolate.py

@@ -70 +70 @@
                 use_cache=False,
                 verbose=False,
-                                output_centroids=False):
+                output_centroids=False):
     """Interpolate vertex_values to interpolation points.
 
@@ -521 +521 @@
                     centroids.append(self.mesh.centroid_coordinates[k])                                         
             else:
-                msg = 'Could not find triangle for point', x
-                raise Exception(msg)
+                 # Mesh has a hole - move this point to the outside list
+#                 num.delete(inside_poly_indices,[i], axis=0)
+                 num.append(outside_poly_indices, [i], axis=0)
+#                msg = 'Could not find triangle for point', x
+#                raise Exception(msg)
+
         return A, inside_poly_indices, outside_poly_indices, centroids