Method BuildTube

Filters/General/vtkHyperStreamline.cxx:657–835

}

int vtkHyperStreamline::BuildTube(vtkDataSet* input, vtkPolyData* output)
{
  vtkHyperPoint *sPrev, *sPtr;
  vtkPoints* newPts;
  vtkFloatArray* newVectors;
  vtkFloatArray* newNormals;
  vtkFloatArray* newScalars = nullptr;
  vtkCellArray* newStrips;
  vtkIdType i, npts, ptOffset = 0;
  int ptId, j, id, k, i1, i2;
  double dOffset, x[3], v[3], s, r, r1[3], r2[3], stepLength;
  double xT[3], sFactor, normal[3], w[3];
  double theta = 2.0 * vtkMath::Pi() / this->NumberOfSides;
  vtkPointData* outPD;
  int iv, ix, iy;
  vtkIdType numIntPts;
  //
  // Initialize
  //
  vtkDebugMacro(<< "Creating hyperstreamline tube");
  if (this->NumberOfStreamers <= 0)
  {
    return 0;
  }

  stepLength = input->GetLength() * this->StepLength;
  outPD = output->GetPointData();

  iv = this->IntegrationEigenvector;
  ix = (iv + 1) % 3;
  iy = (iv + 2) % 3;
  //
  // Allocate
  //
  newPts = vtkPoints::New();
  newPts->Allocate(2500);
  if (input->GetPointData()->GetScalars())
  {
    newScalars = vtkFloatArray::New();
    newScalars->Allocate(2500);
  }
  newVectors = vtkFloatArray::New();
  newVectors->SetNumberOfComponents(3);
  newVectors->Allocate(7500);
  newNormals = vtkFloatArray::New();
  newNormals->SetNumberOfComponents(3);
  newNormals->Allocate(7500);
  newStrips = vtkCellArray::New();
  newStrips->AllocateEstimate(3 * this->NumberOfStreamers, VTK_CELL_SIZE);
  //
  // Loop over all hyperstreamlines generating points
  //
  for (ptId = 0; ptId < this->NumberOfStreamers; ptId++)
  {
    if ((numIntPts = this->Streamers[ptId].GetNumberOfPoints()) < 2)
    {
      continue;
    }