[Paraview] Particles-in-Cells: Two XDMF files?

[Felipe Bordeu]

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

I think the way to doit is to :

>>  by combining them on the <Domain> level in two separate
>> "Temporal Collections" causes ParaView to refuse volume-rendering the
>> fields (not available in the drop down menu).

then use the "Extract Block" to extract the 3DCoRectMesh. The problem is
that the output still a Multi-block Dataset, so you have to use a "Merge
Block" to eliminate the multi-block structure.

Maybe a good idea is to modified the "Extract Block" (or add an option)
to eliminate the multi-block structure when only one mesh is selected.

You can test this creating a "Wavelet" and a "Disk" then group them
using "GroupDataset" and apply the procedure.


Felipe

Le 28/02/2014 10:02, Huebl, Axel a écrit :
> Hey,
> 
> 
> did anyone else try combining a mesh and a polydata data set in *one*
> xdmf file with paraview before?
> 
> Best,
> Axel
> 
> On 24.02.2014 15:45, Huebl, Axel wrote:
>> Hi ParaView mailing list,
>>
>>
>> we are adding paraview/xdmf support to our particle-in-cell code
>> PIConGPU [1] right now. To start with the general topology:
>>
>> The particle in cell algorithm stores field data (scalars and vectors)
>> on a grid, in the simplest case a 3DCoRectMesh. In between the mesh
>> points live particles, described with a Polyvertex topology.
>>
>> The problem that now arises: we would like to describe both grid and
>> polyvertex data in the *same xdmf file*.
>>
>> Doing so, e.g. by combining them on the <Domain> level in two separate
>> "Temporal Collections" causes ParaView to refuse volume-rendering the
>> fields (not available in the drop down menu).
>> Adding two Domains to one file causes ParaView to omit the second data
>> domain during load.
>> Adding both grid types in the same "Temporal Collection" causes a crash.
>>
>> Do we really have to write two xdmf files or did we miss a point?
>>
>> A snippet of our data structures, stored in a hdf5 file per time step:
>>
>> <?xml version="1.0" ?>
>> <Xdmf>
>>   <Domain>
>>     <Grid CollectionType="Temporal" GridType="Collection" Name="Grids">
>>       <Grid GridType="Uniform" Name="Grid_0_0">
>>         <Topology Dimensions="256 768 256" TopologyType="3DCoRectMesh"/>
>>         <Geometry Type="ORIGIN_DXDYDZ">
>>           <DataItem Dimensions="3" Format="XML">0.0 0.0 0.0</DataItem>
>>           <DataItem Dimensions="3" Format="XML">1.0 1.0 1.0</DataItem>
>>         </Geometry>
>>         <Attribute Name="fields/Density_e">
>>           <DataItem Dimensions="256 768 256" Format="HDF"
>> NumberType="Float" Precision="4">h5_0.h5:Density_e</DataItem>
>>         </Attribute>
>>       </Grid>
>>       <!-- .... -->
>>     </Grid>
>>   </Domain>
>> </Xdmf>
>>
>> and the particles:
>>
>> <Xdmf>
>>   <Domain>
>>     <Grid CollectionType="Temporal" GridType="Collection" Name="Polys">
>>       <Grid GridType="Uniform" Name="Poly_0_0">
>>         <Topology NodesPerElement="94109696" TopologyType="Polyvertex"/>
>>         <Attribute Name="particles/e/weighting">
>>           <DataItem Dimensions="94109696" Format="HDF"
>> NumberType="Float" Precision="4">h5_0.h5:weighting</DataItem>
>>         </Attribute>
>>         <Time TimeType="Single" Value="0"/>
>>         <Attribute AttributeType="Vector" Name="particles/e/momentum">
>>           <DataItem Dimensions="94109696 3" Function="JOIN($0,$1,$2)"
>> ItemType="Function">
>>             <DataItem Dimensions="94109696" Format="HDF"
>> NumberType="Float" Precision="4">h5_0.h5:momentum/x</DataItem>
>>             <DataItem Dimensions="94109696" Format="HDF"
>> NumberType="Float" Precision="4">h5_0.h5:momentum/y</DataItem>
>>             <DataItem Dimensions="94109696" Format="HDF"
>> NumberType="Float" Precision="4">h5_0.h5:momentum/z</DataItem>
>>           </DataItem>
>>         </Attribute>
>>       </Grid>
>>       <!-- .... -->
>>     </Grid>
>>   </Domain>
>> </Xdmf>
>>
>>
>> Best regards and thanks for your Feedback,
>> Axel Huebl
>>
>> [1] http://picongpu.hzdr.de
>>
> 
> 
> 
> _______________________________________________
> Powered by 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
> 
> Follow this link to subscribe/unsubscribe:
> http://www.paraview.org/mailman/listinfo/paraview
> 


- -- 
Felipe Bordeu Weldt
Ingénieur de Recherche
- -------------------------------------
Tél. : 33 (0)2 40 37 16 57
Fax. : 33 (0)2 40 74 74 06
Felipe.Bordeu at ec-nantes.fr
Institut GeM - UMR CNRS 6183
École Centrale Nantes
1 Rue de La Noë, 44321 Nantes, FRANCE
- -------------------------------------
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.11 (GNU/Linux)
Comment: Using GnuPG with Thunderbird - http://www.enigmail.net/

iQEcBAEBAgAGBQJTEFWcAAoJEE/fMfNgU9/D+d8IAMR9G5rf92L8yGiwXCzvgT2u
V7II0RiY0x5+nOgZXnHvkeSNOUqPfwkyVlZhrzPDOms0R9oLVRV5hSD/QoLDhacy
GSBz+AveYs471coEwo7ywFLhRnY/J9Oz2YP2lrAcLnsWaruuuQNKB6qRACB6MbE/
7ByHq8cA81KTfx2M1+PS5cZjcviPf5Yg1ZRqBshZabN2Ea2a1BvBMZdcCnOzvOTu
H7FGTylb8fVmCCuPIN3m2hWPfLHrR1DaNKHc/k9ofdqLi8DKjP7za0+ldOh6oM3D
nhyelKcKct4LHcDV4jd3EDH+2hicqqfyDI1gf/dAG0KYwKkZxYg9EBIyNjok6hU=
=fhw9
-----END PGP SIGNATURE-----