Method CreateDataArray

IO/ERF/vtkHDF5Helper.cxx:109–230

----------------------------------------------------------------------------

//----------------------------------------------------------------------------
vtkAbstractArray* vtkHDF5Helper::CreateDataArray(const hid_t& fileId, const std::string& pathName)
{
  vtkSmartPointer<vtkAbstractArray> dataArray;

  if (!vtkHDF5Helper::ArrayExists(fileId, pathName.c_str()))
  {
    vtkWarningWithObjectMacro(nullptr, "Array name '" << pathName << "' isn't available.");
    return nullptr;
  }

  vtkHDF::ScopedH5DHandle arrayId = H5Dopen(fileId, pathName.c_str());
  if (arrayId < 0)
  {
    vtkWarningWithObjectMacro(nullptr, "No array named " << pathName << " available");
    return nullptr;
  }

  hsize_t length = vtkHDF5Helper::GetDataLength(arrayId);
  if (length <= 0)
  {
    return nullptr;
  }

  std::vector<hsize_t> dims = vtkHDF5Helper::GetDataDimensions(arrayId);
  if (dims.size() > 2)
  {
    vtkWarningWithObjectMacro(nullptr, "Only 1D or 2D array are supported.");
    return nullptr;
  }

  // Now the type
  hid_t rawType = H5Dget_type(arrayId);
  hid_t dataType = H5Tget_native_type(rawType, H5T_DIR_ASCEND);

  if (H5Tequal(dataType, H5T_NATIVE_INT))
  {
    dataArray = vtkIntArray::New();
    dataArray->SetNumberOfTuples(dims[0]);
    if (dims.size() == 2)
    {
      dataArray->SetNumberOfComponents(dims[1]);
      dataArray->SetNumberOfTuples(dims[0]);
    }

    int* arrayPtr = static_cast<int*>(vtkArrayDownCast<vtkIntArray>(dataArray)->GetPointer(0));
    H5Dread(arrayId, dataType, H5S_ALL, H5S_ALL, H5P_DEFAULT, arrayPtr);
  }
  else if (H5Tequal(dataType, H5T_NATIVE_FLOAT))
  {
    dataArray = vtkFloatArray::New();
    dataArray->SetNumberOfTuples(dims[0]);
    if (dims.size() == 2)
    {
      dataArray->SetNumberOfComponents(dims[1]);
      dataArray->SetNumberOfTuples(dims[0]);
    }
    else
    {