[Paraview] slicing a large VTK_POLYHEDRON
Pierre Van Hauwaert
pierre at rtech-engineering.nl
Tue Jan 3 16:46:27 EST 2017
Hi Mathieu and Thomas,
Thank you for the suggestions.
1) I tried the tetrahedralize filter on my data but I did not get the
results I was expecting. I have a polyhedron which is not convex
(https://postimg.org/image/9s3egl8rz/) and the filter fill the gap as
you can see here (https://postimg.org/image/5u1nt6d2h/) and it is not
what I was expecting. Indeed, the initial purpose of the slice was to
visualize the hole.
2) I tried this fix:
https://gitlab.kitware.com/vtk/vtk/merge_requests/2088
<https://gitlab.kitware.com/vtk/vtk/merge_requests/2088>
only modifying the 2 files (*Common/DataModel/vtkPolyhedron.cxx*
<https://gitlab.kitware.com/vtk/vtk/merge_requests/2088/diffs#15fe983b62f8470f9429f461e8b010d03b476c0c>
*Common/Core/vtkSetGet.h) *
<https://gitlab.kitware.com/vtk/vtk/merge_requests/2088/diffs#6df1a223bec663530f25d2d214b6b11611fb9cdf>
And compile VTK alone. I used a python script (activating the python VTK
wrapping) to do the test with the vtkCutter function but I ended having
the same exact issue (segfault, no error message) with only the
polyhedron that have more than 1024 faces. I am not sure I should get
the same results if I manage to compile paraview. When I will manage I
will get back to you.
Thanks,
Pierre
On 03-01-17 16:51, Mathieu Westphal wrote:
> Hi
>
> Indeed, incorrect copy paste. thanks for pointing it out.
>
> Mathieu Westphal
>
> On Tue, Jan 3, 2017 at 4:46 PM, T.J. Corona <tj.corona at kitware.com
> <mailto:tj.corona at kitware.com>> wrote:
>
> Hi Mathieu and Pierre,
>
> Perhaps you meant to point to this “work in progress” branch?
>
> https://gitlab.kitware.com/vtk/vtk/merge_requests/2088
> <https://gitlab.kitware.com/vtk/vtk/merge_requests/2088>
>
> Sincerely,
> T.J.
>
> Thomas J. Corona, Ph.D.
> Kitware, Inc.
> R&D Engineer
> 21 Corporate Drive
> Clifton Park, NY 12065-8662
> Phone: 518-881-4443 <tel:%28518%29%20881-4443>
>
>> On Jan 3, 2017, at 10:18 AM, Mathieu Westphal
>> <mathieu.westphal at kitware.com
>> <mailto:mathieu.westphal at kitware.com>> wrote:
>>
>> Hello
>>
>> for your information there is a bug in the slicing of
>> vtkPolyHedron that can cause a segfault.
>> It looks very much like your error. It has yet to be corrected.
>> https://gitlab.kitware.com/vtk/vtk/issues/16877
>> <https://gitlab.kitware.com/vtk/vtk/issues/16877>
>>
>> You can already test this "work in progress" branch to see if
>> this fixes your issue :
>> https://gitlab.kitware.com/vtk/vtk/merge_requests/2304
>> <https://gitlab.kitware.com/vtk/vtk/merge_requests/2304>
>>
>> You could of course modify your data to tetrahedron with
>> tetrahedralize filter.
>>
>> Regards,
>>
>> Mathieu Westphal
>>
>> On Tue, Jan 3, 2017 at 4:10 PM, Pierre Van Hauwaert
>> <pierre at rtech-engineering.nl
>> <mailto:pierre at rtech-engineering.nl>> wrote:
>>
>> Hi,
>>
>> I am using paraview to visualise my data. I generate my data
>> using a combination of the 2 following types
>> (http://www.vtk.org/doc/nightly/html/vtkCellType_8h_source.html
>> <http://www.vtk.org/doc/nightly/html/vtkCellType_8h_source.html>)
>> VTK_VOXEL
>> <http://www.vtk.org/doc/nightly/html/vtkCellType_8h.html#ab1d6fd1f3177b8a2a32bb018807151f8a06a70888292dc4309e42404af594fbb4>=
>> 11
>> VTK_POLYHEDRON = 42 I can visualize that data with parview
>> without a problem. But, because I ended up having a segfault
>> without any error message when slicing (Z=0) the data. I
>> decomposed my data in separate files for each VTK_POLYHEDRON
>> in order to investigate the problem.
>>
>> I found out that only the polyhedrons with the largest size
>> were having the problem. The limit is somewhere between the 2
>> files:
>> - poly-F1000-2395.vtu : working: 466 points, 928 faces
>> - poly-F1000-2987.vtu : segfault when slicing: 546 points,
>> 1088 faces
>>
>> So my guess is that there is a limit for the number of points
>> or the number of faces that can have the polyhedron (1024 ?)
>> if I want to be able to use the slice function.
>>
>> My questions are :
>> 1) Is there actually a limitation in Paraview (VTK?)
>> regarding the number of faces or the number of the point that
>> the slice function can handle for a polyhedron ?
>> 2) If there is this limitation, is there an option to remove
>> it in Paraview ? (specific version, compilation option, etc)
>> 3) If there is this limitation, is there any work around that
>> exist so I am able to slice my data ? For example a different
>> format ?
>>
>>
>> I attached the data to be able to back up my statement:
>>
>> poly-F1000-*.vtu : files describing a polyhedron each:
>> poly-F1000.pvd : load all those files
>>
>> Biggest file for which the slicing works (size, name):
>> 34K poly-F1000-2395.vtu
>>
>> The 8 biggest files triggering the crash when doing a slice:
>> (size, name)
>> 38K poly-F1000-2987.vtu
>> 39K poly-F1000-2983.vtu
>> 39K poly-F1000-2935.vtu
>> 39K poly-F1000-2988.vtu
>> 39K poly-F1000-2931.vtu
>> 40K poly-F1000-2937.vtu
>> 40K poly-F1000-2981.vtu
>> 40K poly-F1000-2930.vtu
>> fail.pvd: load those 8 files
>>
>> the other files can be sliced without any issues
>>
>>
>> Thanks,
>>
>> Pierre
>>
>> --
>> R.Tech Engineering B.V.
>> Eekholt 42 - 1112 XH Diemen - the Netherlands
>> Tel: +31 (0) 2 04 95 02 22
>>
>> _______________________________________________
>> Powered by www.kitware.com <http://www.kitware.com/>
>>
>> Visit other Kitware open-source projects at
>> http://www.kitware.com/opensource/opensource.html
>> <http://www.kitware.com/opensource/opensource.html>
>>
>> Please keep messages on-topic and check the ParaView Wiki at:
>> http://paraview.org/Wiki/ParaView
>> <http://paraview.org/Wiki/ParaView>
>>
>> Search the list archives at:
>> http://markmail.org/search/?q=ParaView
>> <http://markmail.org/search/?q=ParaView>
>>
>> Follow this link to subscribe/unsubscribe:
>> http://public.kitware.com/mailman/listinfo/paraview
>> <http://public.kitware.com/mailman/listinfo/paraview>
>>
>>
>> _______________________________________________
>> Powered by www.kitware.com <http://www.kitware.com>
>>
>> Visit other Kitware open-source projects at
>> http://www.kitware.com/opensource/opensource.html
>> <http://www.kitware.com/opensource/opensource.html>
>>
>> Please keep messages on-topic and check the ParaView Wiki at:
>> http://paraview.org/Wiki/ParaView <http://paraview.org/Wiki/ParaView>
>>
>> Search the list archives at:
>> http://markmail.org/search/?q=ParaView
>> <http://markmail.org/search/?q=ParaView>
>>
>> Follow this link to subscribe/unsubscribe:
>> http://public.kitware.com/mailman/listinfo/paraview
>> <http://public.kitware.com/mailman/listinfo/paraview>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/paraview/attachments/20170103/08cf7d03/attachment.html>
More information about the ParaView
mailing list