[Paraview] Pb reading a Tensor6 XDMF type from HDF5 file

[Romain Beucher]

Hi All,

I have a HDF5 file that contains components of a 2D symmetric tensor in 
an array with dimensions (N, 3) 1:XX, 2:YY, 3:XY

I am trying to create an XDMF Tensor6 type to visualize the field in 
paraview.

Tensor6 has 6 components so I am creating components with zero values 
(XZ, ZZ, YZ) using a JOIN function and a series of hyperslabs...

       <Attribute Type="Tensor6" Center="Node" Name="projStressField">
         <DataItem ItemType="Function" Dimensions="49665 6" 
Function="JOIN($0, 0*$0, $2, $1, 0*$0, 0*$0)">
           <DataItem ItemType="HyperSlab" Dimensions="49665 1" Name="XX">
             <DataItem Dimensions="3 2" Format="XML"> 0 0 1 1 49665 1 
</DataItem>
             <DataItem Format="HDF" NumberType="Float" Precision="8" 
Dimensions="49665 3">projStressField-133.h5:/data</DataItem>
           </DataItem>
           <DataItem ItemType="HyperSlab" Dimensions="49665 1" Name="YY">
             <DataItem Dimensions="3 2" Format="XML"> 0 1 1 1 49665 1 
</DataItem>
             <DataItem Format="HDF" NumberType="Float" Precision="8" 
Dimensions="49665 3">projStressField-133.h5:/data</DataItem>
           </DataItem>
           <DataItem ItemType="HyperSlab" Dimensions="49665 1" Name="XY">
             <DataItem Dimensions="3 2" Format="XML"> 0 2 1 1 49665 1 
</DataItem>
             <DataItem Format="HDF" NumberType="Float" Precision="8" 
Dimensions="49665 3">projStressField-133.h5:/data</DataItem>
           </DataItem>
         </DataItem>
       </Attribute>


Paraview 5.3.0 crashes with the following error:

XDMF Error in 
/home/romain/Downloads/paraview/VTK/ThirdParty/xdmf2/vtkxdmf2/libsrc/XdmfHDF.cxx 
line 558 (Source and Target Spaces specify different sizes)
XDMF Error in 
/home/romain/Downloads/paraview/VTK/ThirdParty/xdmf2/vtkxdmf2/libsrc/XdmfHDF.cxx 
line 559 (Source = 148995 items)
XDMF Error in 
/home/romain/Downloads/paraview/VTK/ThirdParty/xdmf2/vtkxdmf2/libsrc/XdmfHDF.cxx 
line 560 (Target = 99330 items)
XDMF Error in 
/home/romain/Downloads/paraview/VTK/ThirdParty/xdmf2/vtkxdmf2/libsrc/XdmfHDF.cxx 
line 905 (Can't Read Temp Dataset)
paraview: malloc.c:2399: sysmalloc: Assertion `(old_top == initial_top 
(av) && old_size == 0) || ((unsigned long) (old_size) >= MINSIZE && 
prev_inuse (old_top) && ((unsigned long) old_end & (pagesize - 1)) == 
0)' failed.
Aborted (core dumped)

I have attached an example.

Any ideas?

Romain Beucher


-------------- next part --------------
A non-text attachment was scrubbed...
Name: XDMF.fields.00133.xmf
Type: text/xml
Size: 3434 bytes
Desc: XDMF.fields.00133.xmf
URL: <https://public.kitware.com/pipermail/paraview/attachments/20180503/218b72d3/attachment.xml>