[vtkusers] paraview crashes in VTK with > 2^31 points

[Marco Ehlert]

Hello,

does paraview (resp. VTK) support data sets of
more than 2^31 points on 64bit x86 ?

I am using a raw file with dimension
1630x1250x1520 of 8bit-points.

During file reading paraview crashes with SIGFAULT.
The crash seems based on a wrong 'int' calculation in

    vtkImageData::AllocateScalars().

After first attempt of fixing it roughly

/*
     scalars->SetNumberOfTuples((extent[1] - extent[0] + 1)*
                                (extent[3] - extent[2] + 1)*
                                (extent[5] - extent[4] + 1));
*/
     scalars->SetNumberOfTuples((long)(extent[1] - extent[0] + 1)*
                                (long)(extent[3] - extent[2] + 1)*
                                (long)(extent[5] - extent[4] + 1));


the file is read fine now, but paraview does not display
any points. On the other hand arrays of 400x400x400
are displayed perfectly. Maybe there are more such bugs
in paraview components. At least the above bug is still
present in current CVS VTK tree.

Any idea how to debug this situation and fix it completely
is appreciated.

Regards,
Marco Ehlert

IT-Support                        <Marco.Ehlert at mpikg.mpg.de>
Max-Planck-Institute of Colloids  Tel:  +49-(0)331-567-9650
and Interfaces
Am Muehlenberg 1, D-14476, Golm, Germany