Method ProduceOutputFromTextFileDouble

IO/Geometry/vtkParticleReader.cxx:347–420

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

//------------------------------------------------------------------------------
int vtkParticleReader::ProduceOutputFromTextFileDouble(vtkInformationVector* outputVector)
{
  // Get the size of the file.
  this->File->seekg(0, ios::end);
  if (this->File->fail())
  {
    vtkErrorMacro("Could not seek to end of file.");
    return 0;
  }

  size_t fileLength = (unsigned long)this->File->tellg();
  size_t bytesRead = 0;

  std::string s;

  vtkSmartPointer<vtkPoints> points = vtkSmartPointer<vtkPoints>::New();
  points->SetDataTypeToDouble();
  points->Reset();

  vtkSmartPointer<vtkDoubleArray> scalars = vtkSmartPointer<vtkDoubleArray>::New();
  scalars->Reset();
  scalars->SetName("Scalar");

  this->File->seekg(0, ios::beg);

  this->Alliquot = fileLength / quantum;
  this->Count = 1;
  ParseLine<double> pl;
  char buffer[256];
  while (this->File->getline(buffer, 256, '\n'))
  {
    s = buffer;
    if (!s.empty())
    {
      bytesRead += s.size();
      this->DoProgressUpdate(bytesRead, fileLength);
      double val[4];
      val[0] = val[1] = val[2] = val[3] = 0;
      if (pl(s, val))
      {
        points->InsertNextPoint(val[0], val[1], val[2]);
        if (this->HasScalar)
        {
          scalars->InsertNextValue(val[3]);
        }
      }
    }
  }

  vtkSmartPointer<vtkCellArray> vertices = vtkSmartPointer<vtkCellArray>::New();
  vertices->Reset();

  this->NumberOfPoints = points->GetNumberOfPoints();
  for (vtkIdType j = 0; j < (vtkIdType)this->NumberOfPoints; ++j)
  {
    vertices->InsertNextCell(1);
    vertices->InsertCellPoint(j);
  }