Changeset 7682


Ignore:
Timestamp:
Apr 12, 2010, 7:00:43 PM (10 years ago)
Author:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/fit_interpolate/interpolate.py

    r7679 r7682  
    7070                use_cache=False,
    7171                verbose=False,
    72                                 output_centroids=False):
     72                output_centroids=False):
    7373    """Interpolate vertex_values to interpolation points.
    7474
     
    521521                    centroids.append(self.mesh.centroid_coordinates[k])                                         
    522522            else:
    523                 msg = 'Could not find triangle for point', x
    524                 raise Exception(msg)
     523                 # Mesh has a hole - move this point to the outside list
     524#                 num.delete(inside_poly_indices,[i], axis=0)
     525                 num.append(outside_poly_indices, [i], axis=0)
     526#                msg = 'Could not find triangle for point', x
     527#                raise Exception(msg)
     528
    525529        return A, inside_poly_indices, outside_poly_indices, centroids
    526530
Note: See TracChangeset for help on using the changeset viewer.