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
<?xml version="1.0"?> <Xdmf xmlns:xi="http://www.w3.org/2001/XInclude"; Version="2.0"> <Domain> <Grid Name="FEM_Mesh_mesh"> <Time Value="13.3 megayear"/> <Topology Type="Quadrilateral" NumberOfElements="49152"> <DataItem ItemType="Function" Dimensions="49152 4" Function="JOIN($0, $1, $3, $2)"> <DataItem ItemType="HyperSlab" Dimensions="49152 1" Name="C0"> <DataItem Dimensions="3 2" Format="XML"> 0 0 1 1 49152 1 </DataItem> <DataItem Format="HDF" NumberType="Int" Dimensions="49152 1">mesh.h5:/en_map</DataItem> </DataItem> <DataItem ItemType="HyperSlab" Dimensions="49152 1" Name="C1"> <DataItem Dimensions="3 2" Format="XML"> 0 1 1 1 49152 1 </DataItem> <DataItem Format="HDF" NumberType="Int" Dimensions="49152 1">mesh.h5:/en_map</DataItem> </DataItem> <DataItem ItemType="HyperSlab" Dimensions="49152 1" Name="C2"> <DataItem Dimensions="3 2" Format="XML"> 0 2 1 1 49152 1 </DataItem> <DataItem Format="HDF" NumberType="Int" Dimensions="49152 1">mesh.h5:/en_map</DataItem> </DataItem> <DataItem ItemType="HyperSlab" Dimensions="49152 1" Name="C3"> <DataItem Dimensions="3 2" Format="XML"> 0 3 1 1 49152 1 </DataItem> <DataItem Format="HDF" NumberType="Int" Dimensions="49152 1">mesh.h5:/en_map</DataItem> </DataItem> </DataItem> </Topology> <Geometry Type="XYZ"> <DataItem ItemType="Function" Dimensions="49665 3" Function="JOIN($0, $1, 0*$1)"> <DataItem ItemType="HyperSlab" Dimensions="49665 1" Name="XCoords"> <DataItem Dimensions="3 2" Format="XML"> 0 0 1 1 49665 1 </DataItem> <DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="49665 2">mesh.h5:/vertices</DataItem> </DataItem> <DataItem ItemType="HyperSlab" Dimensions="49665 1" Name="YCoords"> <DataItem Dimensions="3 2" Format="XML"> 0 1 1 1 49665 1 </DataItem> <DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="49665 2">mesh.h5:/vertices</DataItem> </DataItem> </DataItem> </Geometry> <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> </Grid> </Domain> </Xdmf>
_______________________________________________ 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 Search the list archives at: http://markmail.org/search/?q=ParaView Follow this link to subscribe/unsubscribe: https://public.kitware.com/mailman/listinfo/paraview
