Method FindPointsWithinRadius

Common/DataModel/vtkPointLocator.cxx:754–801

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

//------------------------------------------------------------------------------
void vtkPointLocator::FindPointsWithinRadius(double R, const double x[3], vtkIdList* result)
{
  this->BuildLocator();
  if (this->HashTable == nullptr)
  {
    return;
  }
  int i, j;
  double dist2;
  double pt[3];
  vtkIdType ptId, nids, cno;
  vtkIdList* ptIds;
  int ijk[3], *nei;
  double R2 = R * R;
  vtkNeighborPoints buckets;

  //  Find bucket point is in.
  this->GetBucketIndices(x, ijk);

  // get all buckets within a distance
  this->GetOverlappingBuckets(&buckets, x, ijk, R, 0);
  // add the original bucket
  buckets.InsertNextPoint(ijk);

  // clear out the result
  result->Reset();

  for (i = 0; i < buckets.GetNumberOfNeighbors(); i++)
  {
    nei = buckets.GetPoint(i);
    cno = nei[0] + nei[1] * this->XD + nei[2] * this->SliceSize;

    if ((ptIds = this->HashTable[cno]) != nullptr)
    {
      nids = ptIds->GetNumberOfIds();
      for (j = 0; j < nids; j++)
      {
        ptId = ptIds->GetId(j);
        this->DataSet->GetPoint(ptId, pt);
        dist2 = vtkMath::Distance2BetweenPoints(x, pt);
        if (dist2 <= R2)
        {
          result->InsertNextId(ptId);
        }
      }
    }
  }
}

//------------------------------------------------------------------------------
void vtkPointLocator::BuildLocator()