Method GetPointFieldAtTimeStep

IO/Geometry/vtkOpenFOAMReader.cxx:9608–9854

------------------------------------------------------------------------------ Read point field at a timestep

//------------------------------------------------------------------------------
// Read point field at a timestep
void vtkOpenFOAMReaderPrivate::GetPointFieldAtTimeStep(const std::string& varName)
{
  // Where to map data
  vtkUnstructuredGrid* internalMesh = this->InternalMesh;
  vtkMultiBlockDataSet* boundaryMesh = this->BoundaryMesh;

  // Boundary information
  const auto& patches = this->BoundaryDict;

  vtkFoamIOobject io(this->CasePath, this->Parent);
  vtkFoamDict dict;
  if (!this->ReadFieldFile(io, dict, varName, this->Parent->PointDataArraySelection))
  {
#if VTK_OPENFOAM_TIME_PROFILING
    this->RequestDataTimeInMicroseconds += io.TimeInMicroseconds;
    this->RequestDataBytes += io.Bytes;
#endif
    return;
  }

  if (io.GetClassName().compare(0, 5, "point") != 0)
  {
    vtkErrorMacro(<< io.GetFileName() << " is not a pointField");
#if VTK_OPENFOAM_TIME_PROFILING
    this->RequestDataTimeInMicroseconds += io.TimeInMicroseconds;
    this->RequestDataBytes += io.Bytes;
#endif
    return;
  }

  // Eg, from "pointScalarField" -> SCALAR_TYPE
  const auto fieldDataType(vtkFoamTypes::FieldToEnum(io.GetClassName(), 5));

  // -------------------------
  // Handle dictionary lookups first

  // The "dimensions" entry - stringify
  const std::string dimString(this->ConstructDimensions(dict));

  // The "internalField" entry
  vtkFoamEntry* ifieldEntry = nullptr;
  {
    ifieldEntry = dict.Lookup("internalField");
    if (ifieldEntry == nullptr)
    {
      vtkErrorMacro(<< "internalField not found in " << io.GetFileName());
#if VTK_OPENFOAM_TIME_PROFILING
      this->RequestDataTimeInMicroseconds += io.TimeInMicroseconds;
      this->RequestDataBytes += io.Bytes;
#endif
      return;
    }
    else if (ifieldEntry->FirstValue().GetType() == vtkFoamToken::EMPTYLIST)
    {
      if (this->NumPoints)
      {
        vtkErrorMacro(<< "internalField of " << io.GetFileName() << " is empty");
#if VTK_OPENFOAM_TIME_PROFILING