Changeset 2266

Timestamp:

Jan 23, 2006, 1:36:15 PM (19 years ago)

Author:

jack

Message:

Heavy looping now in external C module.

Location:

inundation/pyvolution

Files:

• vtk_realtime_visualiser.py (modified) (3 diffs)
• vtk_visualiser_ext.cpp (added)

inundation/pyvolution/vtk_realtime_visualiser.py

@@ -7 +7 @@
 from Numeric import *
 from vtk.tk.vtkTkRenderWidget import vtkTkRenderWidget
+from vtk_visualiser_ext import make_vtkpoints
 
 class Visualiser(threading.Thread):
@@ -137 +138 @@
 
     def draw_quantity(self, q):
-        qty_index = zeros(self.N_vert, Float)
-        for n in range(self.N_tri):
-            for v in range(3):
-                qty_index[self.domain.triangles[n][v]] = self.domain.quantities[q].vertex_values[n][v]
-
-        self.grids[q] = vtk.vtkPoints()
-        grid = self.grids[q]
-
         if self.scale_z.has_key(q):
             scale = self.scale_z[q]
@@ -150 +143 @@
             scale = self.default_scale_z
 
-        for v in range(self.N_vert):
-            grid.InsertNextPoint(self.vert_index[v][0],
-                                 self.vert_index[v][1],
-                                 qty_index[v] * scale)
+        self.grids[q] = make_vtkpoints(self.N_tri,
+                                   self.N_vert,
+                                   scale,
+                                   self.domain.quantities[q].vertex_values,
+                                   self.vert_index,
+                                   self.domain.triangles)
+        grid = self.grids[q]
+        ###################################################
+        # Obsoleted by vtk_visualiser_ext and make_vtkpoints
+        ###################################################
+        #qty_index = zeros(self.N_vert, Float)
+        #for n in range(self.N_tri):
+        #    for v in range(3):
+        #        qty_index[self.domain.triangles[n][v]] = self.domain.quantities[q].vertex_values[n][v]
+        #self.grids[q] = vtk.vtkPoints()
+        #grid = self.grids[q]
+        #for v in range(self.N_vert):
+        #    grid.InsertNextPoint(self.vert_index[v][0],
+        #                         self.vert_index[v][1],
+        #                         qty_index[v] * scale)
+        ###################################################
 
         # Can't recycle vtkPolyData objects: Apparently they behave