Method PopulateMeshTimeIndices

IO/Geometry/vtkOpenFOAMReader.cxx:6627–6683

------------------------------------------------------------------------------ create a Lookup Table containing the location of the points and faces files for each time steps mesh

// create a Lookup Table containing the location of the points
// and faces files for each time steps mesh
void vtkOpenFOAMReaderPrivate::PopulateMeshTimeIndices(const bool skipComputingMetaData)
{
  auto& faces = this->PolyMeshTimeIndexFaces;
  auto& points = this->PolyMeshTimeIndexPoints;

  // Ensure consistent sizing
  const vtkIdType nTimes = this->TimeValues->GetNumberOfTuples();
  if (!skipComputingMetaData)
  {
    this->PopulateMeshIndicesFileChecks->SetNumberOfComponents(3);
    this->PopulateMeshIndicesFileChecks->SetNumberOfTuples(nTimes);
    this->PopulateMeshIndicesFileChecks->SetName("PopulateMeshIndicesFileChecks");
  }

  faces.resize(nTimes, TIMEINDEX_UNVISITED);
  points.resize(nTimes, TIMEINDEX_UNVISITED);

#if VTK_OPENFOAM_TIME_PROFILING
  auto start = std::chrono::high_resolution_clock::now();
#endif
  for (vtkIdType timeIter = 0; timeIter < nTimes; ++timeIter)
  {
    // The mesh directory for this timestep
    const std::string meshDir(this->TimeRegionPath(timeIter) + "/polyMesh/");
    bool hasMeshDir;
    bool topoChanged;
    bool pointsMoved;
    if (!skipComputingMetaData)
    {
      hasMeshDir = vtksys::SystemTools::FileIsDirectory(meshDir);
      topoChanged = hasMeshDir && vtkFoamFile::IsFile(meshDir + "faces", true);
      pointsMoved = hasMeshDir && vtkFoamFile::IsFile(meshDir + "points", true);
      this->PopulateMeshIndicesFileChecks->SetValue(3 * timeIter, hasMeshDir);
      this->PopulateMeshIndicesFileChecks->SetValue(3 * timeIter + 1, topoChanged);
      this->PopulateMeshIndicesFileChecks->SetValue(3 * timeIter + 2, pointsMoved);
    }
    else
    {
      hasMeshDir = this->PopulateMeshIndicesFileChecks->GetValue(3 * timeIter);
      topoChanged = this->PopulateMeshIndicesFileChecks->GetValue(3 * timeIter + 1);
      pointsMoved = this->PopulateMeshIndicesFileChecks->GetValue(3 * timeIter + 2);
    }

    UpdateTimeInstance(faces, timeIter, topoChanged);
    UpdateTimeInstance(points, timeIter, pointsMoved);
  }
#if VTK_OPENFOAM_TIME_PROFILING
  auto end = std::chrono::high_resolution_clock::now();
  this->RequestInformationTimeInMicroseconds +=
    std::chrono::duration_cast<std::chrono::microseconds>(end - start).count();
#endif

#if VTK_FOAMFILE_DEBUG
  PrintMeshTimes("faces", faces);
  PrintMeshTimes("points", points);
#endif
}