visualize boxes in paraview

[Christian Jakob]

VTKRecorder should export boxes for post processing with Paraview.

[Christian Jakob (jakob-ifgt)]

I would try to implement, if it is not too complicated ... Surely need some help on it!

[Anton Gladky (gladky-anton)]

Hi Christian,

I would recommend you to have a look, how facets are exporting [1], [2].
It should be something similar.

Be careful, VTKRecorder is a little messy.

Anton

[1] https://github.com/yade/trunk/blob/master/pkg/dem/VTKRecorder.cpp#L138
[2] https://github.com/yade/trunk/blob/master/pkg/dem/VTKRecorder.cpp#L345

[Christian Jakob (jakob-ifgt)]

I already started ...

for facets vtkTriangle is used, but boxes are cuboids. There is no vtkCuboid for that, so what to do?
Divide the cuboid into rectangles and use vtkRectangle?

[Christian Jakob (jakob-ifgt)]

Do we need this boxesCells for vtk?
It should contain id of the body and id of the points (in this case 8 points per body). I dont know how to implement it, in interactions intrsCells is based on vtkLine (2 points), while in facets facetsCells is based on vtkTriangle (3 points). So what to do with the 8 points from a cuboid?

Please give me a hint how to proceed in lines 406,413,418 and 514 here:

https://github.com/jakob-ifgt/trunk/blob/7d204b604218ea6c3e11a91f982ea2290e46c734/pkg/dem/VTKRecorder.cpp

Thanks in advance,

Christian

[Anton Gladky (gladky-anton)]

It seems, there is something wrong with an order of setting points.

Anton

[Anton Gladky (gladky-anton)]

It should be fixed in f01916a9.

Christian, I have applied your patches, but decided to fix the issue with rectangles.
If you are not agreed, feel free to provide a patch.

Thanks.

Anton

[Christian Jakob (jakob-ifgt)]

Great job Anton!

I think it is the best way to handle boxes with vtk exporter.