[vtkusers] Re: What's Wrong With This?!?

[Andrea Gavana]

Thank you anyway Goodwin to have looked at the scripts. It seems strange to
me that no one else has had the same problems with VTK + Python. I have
somewhat to assume that either there are no Windows+VTK+Python users, or
none of them is using VTK to build structured grids.
I haven't ever tried to compile a debug version of Python, and I suppose it
would be quite a nightmare on Windows. And, even if a miracle happens and I
am able to find the bug, I don't think it would make any difference: the bug
will stay there until VTK/Python developers find a way to correct it, which
may take weeks or months (years?).
That's a pity. I have left VTK one year ago because I was frustrated. Now I
can't even think to restart.

Thank you very much for your suggestions.

Andrea


On 3/22/06, Goodwin Lawlor <goodwin.lawlor at ucd.ie> wrote:
>
> Hi Andrea,
>
> I've run out of ideas except to compile a debug version of python and
> vtk and step through the script.
>
> I've run Examples\DataManipulation\Tcl\Arrays.tcl and its ok even when I
> change the number of tuples to 200- its seems to be just a python problem.
>
> I can't see anything wrong with your python scripts below... sorry!
>
> Goodwin
>
> Andrea Gavana wrote:
> > Hello Goodwin,
> >
> >     thank you very much for your suggestions. This is what I have tried:
> >
> > 1) Call SetNumberOfComponents() before SetNumberOfTuples() => VTK
> Crashes
> >
> > 2) Take vtkExtractGrid out of the pipeline => VTK Crashes
> >
> > 3) Try using vtkPoints::Allocate() and vtkPoints::SetTuple3(): VTK
> > crashes even before showing the 3D window: the problem is here, as you
> > suggested in your point. I am unable to finish the call to SetTuple3 in
> > a loop if the number of points is greater than 100 or something. I can
> > finish the loop if I choose, i.e., 50 points. I attach the example with
> > SetTuple3 as "prova3.py" (using real simple grid data)
> >
> > 4) I get the same runtime error from Windows if I simply use SetData()
> > instead of SetTuple3() and Allocate(), if I choose more than 100 points.
> > If I use 50 points, for example, I get the 3D window with a nice cube. I
> > attach another example as " prova4.py" with SetData() instead of
> > SetTuple3() and Allocate().
> >
> > Probably I have made some mistake somewhere for the SetTuple3/Allocate
> > example, but I am quite sure about the SetData() example. It is
> > fantastically difficult to find examples of use of VTK in Python on the
> > net...
> >
> > For the samples I used Numpy instead of Numeric because Numpy is much
> > more friendly in creating a 3D grid... I get the error message when I
> > change the values for nx, ny, nz to bigger values.
> >
> > Does anyone have some suggestion for me? It's really complicated to find
> > out what the problem is when I get only a windows error message... at
> > least a Python/VTK error messages would be much more useful...
> >
> > Thank you very much for your help and suggestions.
> >
> > Andrea.
> >
> >
> >
> >
> >
> >
> > ------------------------------------------------------------------------
> >
> > import numpy as N
> > import vtk
> > import wx
> > from vtk.wx.wxVTKRenderWindowInteractor import
> wxVTKRenderWindowInteractor
> >
> > app = wx.App(0)
> >
> > # create regular grid
> > xmin, ymin, zmin = -N.pi, -N.pi, 0
> > xmax, ymax, zmax =  N.pi, N.pi, 2*N.pi
> >
> > nx, ny, nz = 5, 5, 2
> > dx, dy, dz = (xmax-xmin)/(nx-1), (ymax-ymin)/(ny-1), (zmax-zmin)/(nz-1)
> >
> > ZZ, YY, XX = N.mgrid[zmin:zmax+dz:dz, ymin:ymax+dy:dy, xmin:xmax+dx:dx]
> > XX = XX.astype(N.Float32)
> > YY = YY.astype(N.Float32)
> > ZZ = ZZ.astype(N.Float32)
> >
> > # set scalar function f(x,y,z)
> > RR = N.sqrt(XX**2 + YY**2 + ZZ**2)
> > WW = N.cos(2*RR)/(RR + 0.5)
> >
> > xyz = N.zeros( (nx*ny*nz,3), N.Float32 )
> >
> > xyz[:,0] = XX.flatten()
> > xyz[:,1] = YY.flatten()
> > xyz[:,2] = ZZ.flatten()
> >
> > # 3-Dcoordinates
> > vtk_xyz = vtk.vtkFloatArray()
> > vtk_xyz.SetNumberOfTuples(nx*ny*nz)
> > vtk_xyz.SetNumberOfComponents(3)
> >
> > vtk_pts = vtk.vtkPoints()
> > vtk_pts.Allocate(nx*ny*nz, nx*ny*nz)
> > vtk_pts.SetDataTypeToFloat()
> >
> > count = 0
> >
> > for ii in xrange(nx):
> >     for jj in xrange(ny):
> >         for kk in xrange(nz):
> >             vtk_xyz.SetTuple3(count, xyz[count, 0], xyz[count, 1],
> xyz[count, 2])
> >             count = count + 1
> >
> > vtk_pts.SetData(vtk_xyz)
> >
> > # create and fill vtk array
> > vtk_zz = vtk.vtkFloatArray()
> > vtk_zz.SetNumberOfTuples(nx*ny*nz)
> > vtk_zz.SetNumberOfComponents(1)
> > vtk_zz.SetVoidArray(WW.flatten(), nx*ny*nz, 1)
> >
> > # create vtk data
> > grid = vtk.vtkStructuredGrid()
> > grid.SetDimensions(nx, ny, nz)
> > grid.SetPoints(vtk_pts)
> > grid.GetPointData().SetScalars(vtk_zz)
> >
> > # Create Actor
> > surfaceMapper = vtk.vtkDataSetMapper()
> > surfaceMapper.SetInput(grid)
> >
> > surfaceActor = vtk.vtkActor()
> > surfaceActor.SetMapper(surfaceMapper)
> >
> > # show
> > frame = wx.Frame(None, -1, "BlaBla")
> >
> > widget = wxVTKRenderWindowInteractor(frame, -1)
> > hiren = vtk.vtkRenderer()
> > widget.GetRenderWindow().AddRenderer(hiren)
> >
> > hiren.AddActor(surfaceActor)
> >
> > widget.AddObserver("ExitEvent", lambda o,e,f=frame: f.Close())
> > widget.Enable(1)
> >
> > sizer = wx.BoxSizer(wx.VERTICAL)
> > sizer.Add(widget, 1, wx.EXPAND)
> > frame.SetSizer(sizer)
> > sizer.Layout()
> >
> > frame.SetSize((800, 600))
> > frame.Show()
> >
> > app.MainLoop()
> >
> >
> >
> > ------------------------------------------------------------------------
> >
> > import numpy as N
> > import vtk
> > import wx
> > from vtk.wx.wxVTKRenderWindowInteractor import
> wxVTKRenderWindowInteractor
> >
> > app = wx.App(0)
> >
> > # create regular grid
> > xmin, ymin, zmin = -N.pi, -N.pi, 0
> > xmax, ymax, zmax =  N.pi, N.pi, 2*N.pi
> >
> > nx, ny, nz = 5, 5, 2
> > dx, dy, dz = (xmax-xmin)/(nx-1), (ymax-ymin)/(ny-1), (zmax-zmin)/(nz-1)
> >
> > ZZ, YY, XX = N.mgrid[zmin:zmax+dz:dz, ymin:ymax+dy:dy, xmin:xmax+dx:dx]
> > XX = XX.astype(N.Float32)
> > YY = YY.astype(N.Float32)
> > ZZ = ZZ.astype(N.Float32)
> >
> > # set scalar function f(x,y,z)
> > RR = N.sqrt(XX**2 + YY**2 + ZZ**2)
> > WW = N.cos(2*RR)/(RR + 0.5)
> >
> > xyz = N.zeros( (nx*ny*nz,3), N.Float32 )
> >
> > xyz[:,0] = XX.flatten()
> > xyz[:,1] = YY.flatten()
> > xyz[:,2] = ZZ.flatten()
> >
> > # 3-Dcoordinates
> > vtk_xyz = vtk.vtkFloatArray()
> > vtk_xyz.SetNumberOfTuples(nx*ny*nz)
> > vtk_xyz.SetNumberOfComponents(3)
> > vtk_xyz.SetVoidArray(xyz.flatten(), 3*nx*ny*nz, 1)
> >
> > vtk_pts = vtk.vtkPoints()
> > vtk_pts.Allocate(nx*ny*nz, nx*ny*nz)
> > vtk_pts.SetDataTypeToFloat()
> > vtk_pts.SetData(vtk_xyz)
> >
> > # create and fill vtk array
> > vtk_zz = vtk.vtkFloatArray()
> > vtk_zz.SetNumberOfTuples(nx*ny*nz)
> > vtk_zz.SetNumberOfComponents(1)
> > vtk_zz.SetVoidArray(WW.flatten(), nx*ny*nz, 1)
> >
> > # create vtk data
> > grid = vtk.vtkStructuredGrid()
> > grid.SetDimensions(nx, ny, nz)
> > grid.SetPoints(vtk_pts)
> > grid.GetPointData().SetScalars(vtk_zz)
> >
> > # Create Actor
> > surfaceMapper = vtk.vtkDataSetMapper()
> > surfaceMapper.SetInput(grid)
> >
> > surfaceActor = vtk.vtkActor()
> > surfaceActor.SetMapper(surfaceMapper)
> >
> > # show
> > frame = wx.Frame(None, -1, "BlaBla")
> >
> > widget = wxVTKRenderWindowInteractor(frame, -1)
> > hiren = vtk.vtkRenderer()
> > widget.GetRenderWindow().AddRenderer(hiren)
> >
> > hiren.AddActor(surfaceActor)
> >
> > widget.AddObserver("ExitEvent", lambda o,e,f=frame: f.Close())
> > widget.Enable(1)
> >
> > sizer = wx.BoxSizer(wx.VERTICAL)
> > sizer.Add(widget, 1, wx.EXPAND)
> > frame.SetSizer(sizer)
> > sizer.Layout()
> >
> > frame.SetSize((800, 600))
> > frame.Show()
> >
> > app.MainLoop()
> >
> >
> >
> > ------------------------------------------------------------------------
> >
> > _______________________________________________
> > This is the private VTK discussion list.
> > Please keep messages on-topic. Check the FAQ at:
> http://www.vtk.org/Wiki/VTK_FAQ
> > Follow this link to subscribe/unsubscribe:
> > http://www.vtk.org/mailman/listinfo/vtkusers
>
> _______________________________________________
> This is the private VTK discussion list.
> Please keep messages on-topic. Check the FAQ at:
> http://www.vtk.org/Wiki/VTK_FAQ
> Follow this link to subscribe/unsubscribe:
> http://www.vtk.org/mailman/listinfo/vtkusers
>



--
"Imagination Is The Only Weapon In The War Against Reality."

http://xoomer.virgilio.it/infinity77/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.vtk.org/pipermail/vtkusers/attachments/20060323/631334b0/attachment.htm>