Method ProduceOutputFromBinaryFileDouble

IO/Geometry/vtkParticleReader.cxx:502–673

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

//------------------------------------------------------------------------------
int vtkParticleReader::ProduceOutputFromBinaryFileDouble(vtkInformationVector* outputVector)
{

  unsigned long fileLength, start, next, length, ptIdx, cellPtIdx;
  unsigned long cellLength;
  int piece, numPieces;
  double *data, *ptr;

  if (!this->FileName)
  {
    vtkErrorMacro(<< "FileName must be specified.");
    return 0;
  }

  this->OpenFile();

  // Get the size of the header from the size of the image
  this->File->seekg(0, ios::end);
  if (this->File->fail())
  {
    vtkErrorMacro("Could not seek to end of file.");
    return 0;
  }

  fileLength = (unsigned long)this->File->tellg();
  if (this->HasScalar)
  {
    this->NumberOfPoints = fileLength / (4 * sizeof(double));
  }
  else
  {
    this->NumberOfPoints = fileLength / (3 * sizeof(double));
  }

  // get the info object
  vtkInformation* outInfo = outputVector->GetInformationObject(0);

  piece = outInfo->Get(vtkStreamingDemandDrivenPipeline::UPDATE_PIECE_NUMBER());
  numPieces = outInfo->Get(vtkStreamingDemandDrivenPipeline::UPDATE_NUMBER_OF_PIECES());

  if ((unsigned long)numPieces > this->NumberOfPoints)
  {
    numPieces = (int)(this->NumberOfPoints);
  }
  if (numPieces <= 0 || piece < 0 || piece >= numPieces)
  {
    return 0;
  }

  start = static_cast<unsigned long>(piece * this->NumberOfPoints / numPieces);
  next = static_cast<unsigned long>((piece + 1) * this->NumberOfPoints / numPieces);

  length = next - start;

  if (this->HasScalar)
  {
    data = new double[length * 4];
  }