<div>Hello Charl,</div>
<div> </div>
<div> thank you for your answer. Ok, I am using Numpy to generate the grid, and this is not good. Well, how am I supposed to actually represent the data? My real data comes from a grid file built externally by another application: I simply read it in a matrix and try to pass it to VTK as a matrix. What a VTK user should do instead? Should I always pass through the construction of a vtk file (using pyvtk) and then re-read it into memory? This would be quite slow...
</div>
<div> </div>
<div>However, I found a couple of things:</div>
<div> </div>
<div>1) As Charl suggested, using the debug build of VTK, I found that the error comes from line 328 of vtkdataarraytemplate.txx, in which I get:</div>
<div> </div>
<div>"Unhandled exception at 0x019e30db (vtkCommon.dll) in Python.exe"</div>
<div> </div>
<div>This line says:</div>
<div> </div>
<div> for(int j=0; j < this->NumberOfComponents; ++j)<br> {<br> this->Tuple[j] = static_cast<double>(t[j]);<br> }<br> return this->Tuple;</div>
<div> <br> </div>
<div>2) Doing the 2-steps method:</div>
<div> </div>
<div>a) vtks = pyvtk.VtkData(pyvtk.StructuredGrid([nx, ny, nz], MyMatrix))</div>
<div> vtks.tofile("MyFile.vtk", "ascii")</div>
<div> </div>
<div>b) grid = vtk.vtkStructuredGridReader()</div>
<div> grid.SetFileName("MyFile.vtk")</div>
<div> surfaceMapper = vtk.vtkDataSetMapper()</div>
<div> surfaceMapper.SetInput(grid.GetOutput())</div>
<div> </div>
<div>works: VTK doesn't crash. But this method may become very slow if I have, say, 200,000 cells or more (and I have them)...</div>
<div> </div>
<div>Is there other things I could try? </div>
<div> </div>
<div>Thank you very much for your suggestions and help.</div>
<div> </div>
<div>Andrea.</div>