Shuenhoy/vtk-m
1
0
Fork 0
mirror of https://gitlab.kitware.com/vtk/vtk-m synced 2024-09-08 13:23:51 +00:00
Code Issues 2 Actions Packages Projects Releases Wiki Activity

Fix HDF5 float data read bug

Browse Source
This commit is contained in:
oruebel 2022-03-14 15:09:43 -07:00 committed by Gunther H. Weber
parent a30d115789
commit 5aed80d7d3
1 changed files with 21 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
examples/contour_tree_distributed/ContourTreeAppDataIO.h
View File

@ -274,6 +274,22 @@ bool read3DHDF5File(const int& mpi_rank,
}
}
}
else if (H5Tequal(H5Dget_type(dataset), H5T_NATIVE_FLOAT))
{
float data_out[count[0]][count[1]][count[2]]; // output buffer
status = H5Dread(dataset, H5T_NATIVE_FLOAT, memspace, dataspace, H5P_DEFAULT, data_out);
// Copy data to 1D array of the expected ValueType
for (vtkm::Id k = 0; k < count[0]; k++)
{
for (vtkm::Id j = 0; j < count[1]; j++)
{
for (vtkm::Id i = 0; i < count[2]; i++)
{
values[to1DIndex(vtkm::Id3(k, j, i), blockSize)] = ValueType(data_out[k][j][i]);
}
}
}
}
else if (H5Tequal(H5Dget_type(dataset), H5T_NATIVE_INT))
{
int data_out[count[0]][count[1]][count[2]]; // output buffer
@ -306,6 +322,11 @@ bool read3DHDF5File(const int& mpi_rank,
}
}
}
else
{
VTKM_LOG_S(vtkm::cont::LogLevel::Error, "Data type not supported by the example HDF5 reader");
throw "Data type not supported by the example HDF5 reader";
}
}
// Release HDF5 resources
H5Sclose(dataspace);