Method FindCell

Testing/GenericBridge/vtkBridgeDataSet.cxx:389–432

// \post clamped_pcoords: result implies (0<=pcoords[0]<=1 && )

int vtkBridgeDataSet::FindCell(
  double x[3], vtkGenericCellIterator*& cell, double tol2, int& subId, double pcoords[3])
{
  assert("pre: not_empty" && GetNumberOfCells() > 0);
  assert("pre: cell_exists" && cell != nullptr);
  assert("pre: positive_tolerance" && tol2 > 0);

  vtkIdType cellid;
  vtkBridgeCellIterator* it = static_cast<vtkBridgeCellIterator*>(cell);

  double* ignoredWeights = new double[this->Implementation->GetMaxCellSize()];

  cellid = this->Implementation->FindCell(x, nullptr, 0, tol2, subId, pcoords, ignoredWeights);

  delete[] ignoredWeights;
  if (cellid >= 0)
  {
    it->InitWithOneCell(this, cellid); // at end
    it->Begin();
    // clamp:
    int i = 0;
    while (i < 3)
    {
      if (pcoords[i] < 0)
      {
        pcoords[i] = 0;
      }
      else if (pcoords[i] > 1)
      {
        pcoords[i] = 1;
      }
      ++i;
    }
  }

  // A=>B: !A || B
  // result => clamped pcoords
  assert("post: clamped_pcoords" &&
    ((cellid < 0) ||
      (pcoords[0] >= 0 && pcoords[0] <= 1 && pcoords[1] >= 0 && pcoords[1] <= 1 &&
        pcoords[2] >= 0 && pcoords[2] <= 1)));

  return cellid >= 0; // bool
}

//------------------------------------------------------------------------------
// Description: