Method ReadCelFile

IO/Geometry/vtkProStarReader.cxx:244–446

read in the cells from the .cel file ---------------------------------------------------------------------------*\ Line 1: PROSTAR_CELL [newline] Line 2: 0 0 0 0 0 0 0 [newline] Body: [newline] .. .. with shapeId: * 1 = point * 2 = line * 3 = shell * 11 = hexa * 12 = prism * 13 = tetra *

\*---------------------------------------------------------------------------*/
bool vtkProStarReader::ReadCelFile(vtkUnstructuredGrid* output, const idMapping& mapPointId)
{
  FILE* in = this->OpenFile(".cel");
  if (in == nullptr)
  {
    return false;
  }

  constexpr int MAX_LINE = 1024;
  char rawLine[MAX_LINE];

  int errorCount = 0;
  vtk::scan_result_type<std::string_view, int> resultLabel;
  if (fgets(rawLine, MAX_LINE, in) != nullptr && strncmp(rawLine, "PROSTAR_CELL", 12) == 0 &&
    fgets(rawLine, MAX_LINE, in) != nullptr &&
    ((resultLabel = vtk::scan_int<int>(std::string_view(rawLine)))) && resultLabel->value() >= 4000)
  {
    vtkDebugMacro(<< "Got PROSTAR_CELL header");
  }
  else
  {
    vtkErrorMacro(<< "Error reading header for PROSTAR_CELL file");
    ++errorCount;
  }
  int lineLabel = resultLabel->value();

  // don't know the number of cells a priori -- just pick some number
  output->Allocate(10000, 20000);

  // add a cellTableId array
  vtkIntArray* cellTableId = vtkIntArray::New();
  cellTableId->Allocate(10000, 20000);
  cellTableId->SetName("cellTableId");

  int shapeId, nLabels, tableId, typeId;
  std::vector<vtkIdType> starLabels;
  starLabels.reserve(256);

  // face-stream for a polyhedral cell
  // [numFace0Pts, id1, id2, id3, numFace1Pts, id1, id2, id3, ...]
  std::vector<vtkIdType> faceStream;
  faceStream.reserve(256);

  // use string buffer for easier parsing
  std::istringstream strbuf;

  int lineNr = 2;
  while (!errorCount && fgets(rawLine, MAX_LINE, in) != nullptr)
  {
    ++lineNr;
    auto result =
      vtk::scan<int, int, int, int, int>(std::string_view(rawLine), "{:d} {:d} {:d} {:d} {:d}");
    if (result)
    {
      std::tie(lineLabel, shapeId, nLabels, tableId, typeId) = result->values();
      starLabels.clear();
      starLabels.reserve(nLabels);