Changeset 8574

Timestamp:

Sep 13, 2012, 9:02:22 PM (13 years ago)

Author:

steve

Message:

Fixing a problem with node_index used in erosion operator

Location:

trunk/anuga_core/source/anuga

Files:

• abstract_2d_finite_volumes/generic_domain.py (modified) (1 diff)
• operators/erosion_operators.py (modified) (1 diff)

trunk/anuga_core/source/anuga/abstract_2d_finite_volumes/generic_domain.py

@@ -134 +134 @@
         self.number_of_triangles_per_node = \
                                     self.mesh.number_of_triangles_per_node
-
+        self.node_index = self.mesh.node_index
         self.vertex_value_indices = self.mesh.vertex_value_indices
         self.number_of_triangles = self.mesh.number_of_triangles

trunk/anuga_core/source/anuga/operators/erosion_operators.py

@@ -380 +380 @@
         for nid in self.node_ids:
             #print nid,self.domain.number_of_triangles_per_node[nid]
-            for vid in range(node_index[nid], node_index[nid+1]):
+            for vid in range(self.node_index[nid], self.node_index[nid+1]):
                 vidd = self.domain.vertex_value_indices[vid]
                 vertex_ids.append(vidd)