If you can post the file, and/or a stack trace. Someone may be better able to help.
David E DeMarle Kitware, Inc. R&D Engineer 21 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-881-4909 On Thu, Jul 26, 2012 at 11:27 PM, Petr Krysl <[email protected]> wrote: > Hello, > > Paraview crashes when reading a binary VTK file. I am writing out that file > using the Matlab function below. When written out in ASCII the file is read > and interpreted correctly; when written out in binary paraview just crashes > without giving any details as to where and why it crashed. I hope someone > has an idea: after several hours of debugging I am at a loss. > > % Export mesh to a VTK 1.0 file as an unstructured grid. > % > % function vtk_export_mesh (theFile,Connectivity,Points,Cell_types,options) > % > % theFile= File name as a string, > % Connectivity= connectivity array, one row per cell, > % Points= Coordinate array, one row per point, > % Cell_types= Cell type code: L2=3, T3=5, Q4=9, T4=10, H8=12 > function vtk_export_mesh (theFile,Connectivity,Points,Cell_types,options) > binary= ~true; > scalars=[]; scalars_name= ['Data']; > if ( exist ( 'options', 'var') ) > if (isfield(options, 'binary')) > binary= options.binary; > end > if (isfield(options, 'scalars')) > scalars= options.scalars; > end > if (isfield(options, 'scalars_name')) > scalars_name= options.scalars_name; > end > end > if (~iscell(Connectivity)) > Connectivity={Connectivity}; > end > if (~iscell(Cell_types)) > Cell_types={Cell_types}; > end > [pathstr, name, ext] = fileparts(theFile); > if (~strcmp(ext,'vtk')) > theFile = [theFile '.vtk']; > end > fid=fopen(theFile,'w'); > if (fid==-1) > warning (['Could not open ' theFile]) > return > end > fprintf(fid,'# vtk DataFile Version 3.0\n'); > fprintf(fid,'Really cool data\n'); > if (~binary) > fprintf(fid,'ASCII\n'); > else > fprintf(fid,'BINARY\n'); > end > fprintf(fid,'\n'); > fprintf(fid,'DATASET UNSTRUCTURED_GRID\n'); > fprintf(fid,'POINTS %d double\n',size(Points,1));% > if (~binary) > for i= 1:size(Points, 1) > fprintf(fid,'%g %g %g\n',Points(i,:)); > end > else > fwrite(fid,cast(Points,'double'),'double'); > end > fprintf(fid,'\n'); > fprintf(fid,'\n'); > > for k=1:length(Connectivity) > f=Connectivity{k}; > ctype=zeros(size(f, 1),1)+Cell_types{k}; > fprintf(fid,'CELLS %d %d\n',size(f,1),(size(f,1)*(size(f,2)+1)));% > if (~binary) > for i= 1:size(f, 1) > fprintf(fid,'%d ',size(f,2)); > for j= 1: size(f,2) > fprintf(fid,'%d ',f(i,j)-1); > end > fprintf(fid,'\n'); > end > else > > fwrite(fid,cast([zeros(size(f,1),1)+size(f,2),f-1],'int32'),'int32'); > end > fprintf(fid,'\n'); > fprintf(fid,'\n'); > fprintf(fid,'CELL_TYPES %d\n',size(f,1));% > if (~binary) > for i= 1:size(f, 1) > fprintf(fid,'%d\n',ctype(i)); > end > else > fwrite(fid,cast(ctype,'int32'),'int32'); > end > fprintf(fid,'\n'); > fprintf(fid,'\n'); > end > > if (~isempty(scalars)) > fprintf(fid,'point_data %d\n',length(scalars));% > fprintf(fid,'SCALARS %s float\n',scalars_name);% > fprintf(fid,'LOOKUP_TABLE default\n');% > if (~binary) > for j= 1:length(scalars) > fprintf(fid,'%s\n',scalars(j)); > end > else > fwrite(fid,cast(scalars,'single'),'single'); > end > end > fprintf(fid,'\n'); > fid=fclose(fid); > end > > -- > Petr Krysl > University of California, San Diego > Skype: Petr.Krysl.UCSD.EDU > http://hogwarts.ucsd.edu/~pkrysl/ > > > _______________________________________________ > Powered by www.kitware.com > > Visit other Kitware open-source projects at > http://www.kitware.com/opensource/opensource.html > > Please keep messages on-topic and check the ParaView Wiki at: > http://paraview.org/Wiki/ParaView > > Follow this link to subscribe/unsubscribe: > http://www.paraview.org/mailman/listinfo/paraview > _______________________________________________ Powered by www.kitware.com Visit other Kitware open-source projects at http://www.kitware.com/opensource/opensource.html Please keep messages on-topic and check the ParaView Wiki at: http://paraview.org/Wiki/ParaView Follow this link to subscribe/unsubscribe: http://www.paraview.org/mailman/listinfo/paraview
