Method Clip

Common/DataModel/vtkGenericAdaptorCell.cxx:452–650

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

//------------------------------------------------------------------------------
void vtkGenericAdaptorCell::Clip(double value, vtkImplicitFunction* f,
  vtkGenericAttributeCollection* attributes, vtkGenericCellTessellator* tess, int insideOut,
  vtkIncrementalPointLocator* locator, vtkCellArray* connectivity, vtkPointData* outPd,
  vtkCellData* outCd, vtkPointData* internalPd, vtkPointData* secondaryPd, vtkCellData* secondaryCd)
{
  assert("pre: attributes_exist" && attributes != nullptr);
  assert("pre: tessellator_exists" && tess != nullptr);
  assert("pre: locator_exists" && locator != nullptr);
  assert("pre: connectivity_exist" && connectivity != nullptr);
  assert("pre: internalPd_exists" && internalPd != nullptr);
  assert("pre: secondaryPd_exists" && secondaryPd != nullptr);
  assert("pre: secondaryCd_exists" && secondaryCd != nullptr);

  int i;
  int j;
  int c;
  double contVal = -1000;
  double* values;

  vtkCell* linearCell = nullptr;
  vtkIdType ptsCount = 0;

  this->Reset();

  // for each cell-centered attribute: copy the value in the secondary
  // cell data.
  secondaryCd->Reset();
  int attrib = 0;
  while (attrib < attributes->GetNumberOfAttributes())
  {
    if (attributes->GetAttribute(attrib)->GetCentering() == vtkCellCentered)
    {
      vtkDataArray* array = secondaryCd->GetArray(attributes->GetAttribute(attrib)->GetName());
      values = attributes->GetAttribute(attrib)->GetTuple(this);
      array->InsertNextTuple(values);
    }
    attrib++;
  }

  int attribute = this->GetHighestOrderAttribute(attributes);
  if (this->IsGeometryLinear() &&
    (attribute == -1 || this->IsAttributeLinear(attributes->GetAttribute(attribute))))
  {
    // linear case
    switch (this->GetType())
    {
      case VTK_HIGHER_ORDER_TRIANGLE:
        linearCell = this->Triangle;
        ptsCount = 3;
        break;
      case VTK_HIGHER_ORDER_QUAD:
        linearCell = this->Quad;
        ptsCount = 4;
        break;
      case VTK_HIGHER_ORDER_TETRAHEDRON:
        linearCell = this->Tetra;
        ptsCount = 4;
        break;