<div dir="ltr">Hmmm. These look right at first glance. Can you send me the files? Off the list since the list has a size limitation.<div><br></div><div>-berk</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On Tue, Jun 10, 2014 at 12:45 PM, André Malcher <span dir="ltr"><<a href="mailto:andre.malcher@mytum.de" target="_blank">andre.malcher@mytum.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hello,<br>
<br>
I'm looking for resources/solutions to set up an parallel structured grid file (*.pvts).<br>
I have an parallel Lattice-Boltzmann solver written in MPI in which each rank computes one part of a cubic domain.<br>
I was able to output the *.vts files for each sub-domain, yet I am unable to compile those outputs to be shown in Paraview as a composition.<br>
The problem seems to be related to the *.pvts file, but unfortunately it is very hard to find resources on the VTK xml formats in the internet.<br>
Also, the "file-formats.pdf" provided at <a href="http://www.vtk.org" target="_blank">www.vtk.org</a> lacks information about some details.<br>
<br>
To describe my problem and for simplicity, I set up a cube with side length 8 (no. of elements) that is computed by 8 ranks,<br>
i.e. I split the domain into eight cubes (of length 4 then) that are processed by each of the ranks.<br>
My *.vts files have the following format (e.g. for rank 0):<br>
<br>
<?xml version="1.0"?><br>
<VTKFile type="StructuredGrid" version="0.1" byte_order="LittleEndian"><br>
<StructuredGrid WholeExtent="0 3 0 3 0 3" ><br>
<Piece Extent="0 3 0 3 0 3" ><br>
<PointData Scalars="density" Vectors="velocity"><br>
<DataArray Name="velocity" type="Float32" NumberOfComponents="3" format="ascii"><br>
</DataArray><br>
<DataArray Name="density" type="Float32"><br>
<!-- Some values here --><br>
</DataArray><br>
</PointData><br>
<CellData></CellData><br>
<Points><br>
<DataArray type="Float32" Name="coordinate" NumberOfComponents="3" format="ascii"><br>
<!-- Some values here --><br>
</DataArray><br>
</Points><br>
</Piece><br>
</StructuredGrid><br>
</VTKFile><br>
<br>
These work fine and show the expected result when loaded into ParaView.<br>
However, combining these files to one output fails. This is my *.pvts file that I create:<br>
<br>
<?xml version="1.0"?><br>
<VTKFile type="PStructuredGrid" version="0.1" byte_order="LittleEndian"><br>
<PStructuredGrid WholeExtent="0 7 0 7 0 7 " GhostLevel="#"><br>
<PPointData Scalars="density" Vectors="velocity"><br>
<DataArray type="Float32" Name="velocity" NumberOfComponents="3" format="ascii" /><br>
<DataArray type="Float32" Name="density" /><br>
</PPointData><br>
<PCellData></PCellData><br>
<PPoints><br>
<DataArray type="Float32" Name="coordinate" NumberOfComponents="3" format="ascii" /><br>
</PPoints><br>
<Piece Extent="0 3 0 3 0 3" Source="lbm_out_0.0.vts"><br>
</Piece><br>
<Piece Extent="3 7 0 3 0 3" Source="lbm_out_1.0.vts"><br>
</Piece><br>
<Piece Extent="0 3 3 7 0 3" Source="lbm_out_2.0.vts"><br>
</Piece><br>
<Piece Extent="3 7 3 7 0 3" Source="lbm_out_3.0.vts"><br>
</Piece><br>
<Piece Extent="0 3 0 3 3 7" Source="lbm_out_4.0.vts"><br>
</Piece><br>
<Piece Extent="3 7 0 3 3 7" Source="lbm_out_5.0.vts"><br>
</Piece><br>
<Piece Extent="0 3 3 7 3 7" Source="lbm_out_6.0.vts"><br>
</Piece><br>
<Piece Extent="3 7 3 7 3 7" Source="lbm_out_7.0.vts"><br>
</Piece><br>
</PStructuredGrid><br>
</VTKFile><br>
<br>
When I load this into ParaView and click apply it shows the right extents (0-7) but the density/velocity values as well as the bounds/ranges are totally wrong.<br>
I attached a screenshot of the ParaView information window to give a better idea what I am talking about. By the way, since the data is useless,<br>
ParaView also doesn't show any outline after clicking on "Apply". I would appreciate if you could provide me with some details about what I may be doing wrong.<br>
<br>
Besides, I would like to know how the overlapping extent areas work. When they don't overlap, I get error messages and I have read about the fact that they are required to overlap.<br>
But how does the extent include four elements when in (0..3) in the special case described above but five elements in (3..7)?<br>
<br>
And last but not least, is it possible to load legacy *.vtk files as extent source?<br>
<br>
Thank you very much in advance.<br>
André<br>
_______________________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the VTK FAQ at: <a href="http://www.vtk.org/Wiki/VTK_FAQ" target="_blank">http://www.vtk.org/Wiki/VTK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.vtk.org/mailman/listinfo/vtkusers" target="_blank">http://www.vtk.org/mailman/listinfo/vtkusers</a><br>
</blockquote></div><br></div>