Method LoadAttribute

IO/AMR/vtkAMREnzoReaderInternal.cxx:362–573

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

//------------------------------------------------------------------------------
int vtkEnzoReaderInternal::LoadAttribute(const char* attribute, int blockIdx)
{
  // TODO: implement this
  // called by GetBlockAttribute( ... ) or GetParticlesAttribute()
  this->ReadMetaData();

  if (attribute == nullptr || blockIdx < 0 || blockIdx >= this->NumberOfBlocks)
  {
    return 0;
  }

  // this->Internal->Blocks includes a pseudo block --- the root as block #0
  blockIdx++;

  // currently only the HDF5 file format is supported
  std::string blckFile = this->Blocks[blockIdx].BlockFileName;
  hid_t fileIndx = H5Fopen(blckFile.c_str(), H5F_ACC_RDONLY, H5P_DEFAULT);
  if (fileIndx < 0)
  {
    return 0;
  }

  // retrieve the contents of the root directory to look for a group
  // corresponding to the target block, if available, open that group

  char blckName[65];
  int objIndex;
  hsize_t numbObjs;
  hid_t rootIndx = H5Gopen(fileIndx, "/");
  H5Gget_num_objs(rootIndx, &numbObjs);
  for (objIndex = 0; objIndex < static_cast<int>(numbObjs); objIndex++)
  {
    int type = H5Gget_objtype_by_idx(rootIndx, objIndex);
    if (type == H5G_GROUP)
    {
      H5Gget_objname_by_idx(rootIndx, objIndex, blckName, 64);
      auto result = vtk::scan<int>(blckName, "Grid{:d}");
      if (result && (result->value() == blockIdx || result->value() == (blockIdx + 1)))
      {
        // located the target block
        rootIndx = H5Gopen(rootIndx, blckName);
        break;
      }
    }
  } // END for all objects

  // disable error messages while looking for the attribute (if any) name
  // and enable error messages when it is done
  void* pContext = nullptr;
  H5E_auto_t erorFunc;
  H5Eget_auto(&erorFunc, &pContext);
  H5Eset_auto(nullptr, nullptr);

  hid_t attrIndx = H5Dopen(rootIndx, attribute);

  H5Eset_auto(erorFunc, pContext);
  pContext = nullptr;