Method WriteGeometryFile

IO/Geometry/vtkBYUWriter.cxx:148–210

}

void vtkBYUWriter::WriteGeometryFile(FILE* geomFile, int numPts)
{
  int numPolys, numEdges;
  int i;
  double* x;
  vtkIdType npts = 0;
  const vtkIdType* pts = nullptr;
  vtkPoints* inPts;
  vtkCellArray* inPolys;
  vtkPolyData* input = this->GetInput();
  //
  // Check input
  //
  inPolys = input->GetPolys();
  if ((inPts = input->GetPoints()) == nullptr || inPolys == nullptr)
  {
    vtkErrorMacro(<< "No data to write!");
    return;
  }
  //
  // Write header (not using fixed format! - potential problem in some files.)
  //
  numPolys = input->GetPolys()->GetNumberOfCells();
  for (numEdges = 0, inPolys->InitTraversal(); inPolys->GetNextCell(npts, pts);)
  {
    numEdges += npts;
  }

  vtk::print(geomFile, "{:d} {:d} {:d} {:d}\n", 1, numPts, numPolys, numEdges);
  vtk::print(geomFile, "{:d} {:d}\n", 1, numPolys);

  //
  // Write data
  //
  // write point coordinates
  for (i = 0; i < numPts; i++)
  {
    x = inPts->GetPoint(i);
    vtk::print(geomFile, "{:e} {:e} {:e} ", x[0], x[1], x[2]);
    if ((i % 2))
    {
      vtk::print(geomFile, "\n");
    }
  }
  if ((numPts % 2))
  {
    vtk::print(geomFile, "\n");
  }

  // write poly data. Remember 1-offset.
  for (inPolys->InitTraversal(); inPolys->GetNextCell(npts, pts);)
  {
    // write this polygon
    // treating vtkIdType as int
    for (i = 0; i < (npts - 1); i++)
    {
      vtk::print(geomFile, "{:d} ", static_cast<int>(pts[i] + 1));
    }