[Paraview] slicing a large VTK_POLYHEDRON
Pierre Van Hauwaert
pierre at rtech-engineering.nl
Wed Jan 4 04:18:26 EST 2017
Hi,
As requested I attached with that email a list of 8 segfaulting
polyhedrons. It happens with a centre of (0,0,0) and an oZ normal.
Regards,
Pierre
On 01/04/2017 10:11 AM, Mathieu Westphal wrote:
> Hi
>
> It may be a different issue than the one i'm fixing. If you can share
> the segfaulting polyhedron, please do so i can take a look when i get
> back on this issue ( which my not be soon, feel free to find other
> solutions)
>
> Regards,
>
> Mathieu Westphal
>
> On Tue, Jan 3, 2017 at 10:46 PM, Pierre Van Hauwaert
> <pierre at rtech-engineering.nl <mailto:pierre at rtech-engineering.nl>> wrote:
>
> 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/
> <https://postimg.org/image/9s3egl8rz/>) and the filter fill the
> gap as you can see here (https://postimg.org/image/5u1nt6d2h/
> <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)
>>> 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
>>>
>>> Search the list archives at:
>>> 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
>>>
>>> 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
>>>
>>> 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/20170104/5b003498/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: segfaultingPolyhedrons.tgz
Type: application/x-gtar-compressed
Size: 74951 bytes
Desc: not available
URL: <http://public.kitware.com/pipermail/paraview/attachments/20170104/5b003498/attachment.tgz>
More information about the ParaView
mailing list