<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal">Hi,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I have data in hdf5 files which I want to process with Paraview, so I wrote an xdmf file that describes the data. I have compiled paraview 4.0.1 with MPI on a linux cluster, and paraview with Qt on windows which I use as client.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">My data displays properly in the client, but since it is quite large I would like to read/process it in parallel.
<o:p></o:p></p>
<p class="MsoNormal">However, if I color by ProcessId (either the filter or the vtkProcessId) I see that all the data is given by a single process. I believe that everything is set up properly since if I create a sphere source and apply the same filter, I can
indeed see several colors. The structure of the xdmf file is as follows:<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><Xdmf><o:p></o:p></p>
<p class="MsoNormal"> <Domain><o:p></o:p></p>
<p class="MsoNormal"> <Grid Name="Channel" GridType="Collection" CollectionType="Temporal"><o:p></o:p></p>
<p class="MsoNormal"> <Grid Name="G@0s" GridType="Uniform" ><o:p></o:p></p>
<p class="MsoNormal"> <Topology … /><o:p></o:p></p>
<p class="MsoNormal"> <Geometry … /><o:p></o:p></p>
<p class="MsoNormal"> <Attribute Name="Pressure" Center="Cell" AttributeType="Scalar"><o:p></o:p></p>
<p class="MsoNormal"> <DataItem ItemType="HyperSlab" Dimensions="327680" Type="HyperSlab"><o:p></o:p></p>
<p class="MsoNormal"> <DataItem Dimensions="3 2" Format="XML"><o:p></o:p></p>
<p class="MsoNormal"> 0 983040<o:p></o:p></p>
<p class="MsoNormal"> 1 1<o:p></o:p></p>
<p class="MsoNormal"> 1 327680<o:p></o:p></p>
<p class="MsoNormal"> </DataItem><o:p></o:p></p>
<p class="MsoNormal"> <DataItem Format="HDF" Dimensions="71 1310720"><o:p></o:p></p>
<p class="MsoNormal"> Re350_oscillating0001.h5:/snapshots_T<o:p></o:p></p>
<p class="MsoNormal"> </DataItem><o:p></o:p></p>
<p class="MsoNormal"> </DataItem><o:p></o:p></p>
<p class="MsoNormal"> </Attribute><o:p></o:p></p>
<p class="MsoNormal"> </Grid><o:p></o:p></p>
<p class="MsoNormal"> <Grid Name="G@1s" GridType="Uniform" ><o:p></o:p></p>
<p class="MsoNormal"> <Topology … /><o:p></o:p></p>
<p class="MsoNormal"> <Geometry … /><o:p></o:p></p>
<p class="MsoNormal"> <Attribute Name="Pressure" Center="Cell" AttributeType="Scalar"><o:p></o:p></p>
<p class="MsoNormal"> <DataItem ItemType="HyperSlab" Dimensions="327680" Type="HyperSlab"><o:p></o:p></p>
<p class="MsoNormal"> <DataItem Dimensions="3 2" Format="XML"><o:p></o:p></p>
<p class="MsoNormal"> 0 983040<o:p></o:p></p>
<p class="MsoNormal"> 1 1<o:p></o:p></p>
<p class="MsoNormal"> 1 327680<o:p></o:p></p>
<p class="MsoNormal"> </DataItem><o:p></o:p></p>
<p class="MsoNormal"> <DataItem Format="HDF" Dimensions="71 1310720"><o:p></o:p></p>
<p class="MsoNormal"> Re350_oscillating0001.h5:/snapshots_T<o:p></o:p></p>
<p class="MsoNormal"> </DataItem><o:p></o:p></p>
<p class="MsoNormal"> </DataItem><o:p></o:p></p>
<p class="MsoNormal"> </Attribute><o:p></o:p></p>
<p class="MsoNormal"> </Grid><o:p></o:p></p>
<p class="MsoNormal"> <Time TimeType="List"><o:p></o:p></p>
<p class="MsoNormal"> <DataItem Format="XML" NumberType="Float" Dimensions="71"><o:p></o:p></p>
<p class="MsoNormal"> 0 1
<o:p></o:p></p>
<p class="MsoNormal"> </DataItem><o:p></o:p></p>
<p class="MsoNormal"> </Time><o:p></o:p></p>
<p class="MsoNormal"> </Grid><o:p></o:p></p>
<p class="MsoNormal"> </Domain><o:p></o:p></p>
<p class="MsoNormal"></Xdmf><o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">After much searching it seems that one has to use Spatial collections for paraview to deal with the data in parallel. Indeed with something like:<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><Xdmf><o:p></o:p></p>
<p class="MsoNormal"> <Domain><o:p></o:p></p>
<p class="MsoNormal"> <Grid Name="Channel" GridType="Collection" CollectionType="Temporal"><o:p></o:p></p>
<p class="MsoNormal"> <Grid Name="foo" GridType="Collection" CollectionType="Spatial"><o:p></o:p></p>
<p class="MsoNormal"> <Grid Name="G0@0s" GridType="Uniform" ><o:p></o:p></p>
<p class="MsoNormal"> <Topology … *<b>Contains a fraction of the cells</b>* /><o:p></o:p></p>
<p class="MsoNormal"> <Geometry … /><o:p></o:p></p>
<p class="MsoNormal"> <Attribute … /><o:p></o:p></p>
<p class="MsoNormal"> </Grid><o:p></o:p></p>
<p class="MsoNormal"> <Grid Name="G1@0s" GridType="Uniform" ><o:p></o:p></p>
<p class="MsoNormal"> <Topology … *<b>Contains another fraction of the cells</b>* /><o:p></o:p></p>
<p class="MsoNormal"> <Geometry … /><o:p></o:p></p>
<p class="MsoNormal"> <Attribute … /><o:p></o:p></p>
<p class="MsoNormal"> </Grid><o:p></o:p></p>
<p class="MsoNormal"> </Grid><o:p></o:p></p>
<p class="MsoNormal" style="margin-left:72.0pt;text-indent:36.0pt">…<o:p></o:p></p>
<p class="MsoNormal"> <Time TimeType="List"><o:p></o:p></p>
<p class="MsoNormal"> <DataItem Format="XML" NumberType="Float" Dimensions="71"><o:p></o:p></p>
<p class="MsoNormal"> 0 1
<o:p></o:p></p>
<p class="MsoNormal"> </DataItem><o:p></o:p></p>
<p class="MsoNormal"> </Time><o:p></o:p></p>
<p class="MsoNormal"> </Grid><o:p></o:p></p>
<p class="MsoNormal"> </Domain><o:p></o:p></p>
<p class="MsoNormal"></Xdmf><o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Paraview assigns 1 process per sub-mesh G0, G1, … in the spatial collection. The problem is that all these G0, G1, etc. don’t appear to be connected, and when using filters like Cell Data to Point Data, contours etc., there is a discontinuity
that is clearly visible between each part of the mesh so I don’t think that is the way to go. Or am I missing something? What is the usual way to get parallel processing of xdmf files?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Any help would be greatly appreciated,<o:p></o:p></p>
<p class="MsoNormal"><br>
Romain<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>