Method SimplifyBranches

Common/DataModel/vtkReebGraph.cxx:1269–1527

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

//------------------------------------------------------------------------------
int vtkReebGraph::Implementation::SimplifyBranches(
  double simplificationThreshold, vtkReebGraphSimplificationMetric* simplificationMetric)
{
  static const vtkReebLabelTag RouteOld = 100;
  static const vtkReebLabelTag RouteNew = 200;
  int nstack, mstack = 0;
  int* stack = nullptr;

  if (!simplificationThreshold)
    return 0;

  int nsimp = 0;
  int cont = 0;
  constexpr int step = 10000;
  bool redo;

  vtkDataSet* input = inputMesh;

REDO:
  nstack = 0;
  redo = false;

  for (int N = 1; N < this->MainNodeTable.Size; ++N)
  {
    if (this->GetNode(N)->ArcUpId == -2)
      continue;

    vtkReebNode* n = this->GetNode(N);

    // simplify atomic nodes
    if (!n->ArcDownId && !n->ArcUpId)
    {
      // delete the node from the graph...
      this->GetNode(N)->ArcUpId = -2;
      this->GetNode(N)->ArcDownId = this->MainNodeTable.FreeZone;
      this->MainNodeTable.FreeZone = (N);
      --(this->MainNodeTable.Number);
    }
    else if (!n->ArcDownId)
    {
      // insert into stack branches to simplify
      for (int A_ = n->ArcUpId; A_; A_ = this->GetArc(A_)->ArcDwId0)
      {
        vtkReebArc* a_ = this->GetArc(A_);
        if ((!inputMesh) || (!simplificationMetric))
        {
          if (vtkReebGraphGetArcPersistence(this, a_) < simplificationThreshold)
          {
            vtkReebGraphStackPush(A_);
          }
        }
        else
        {
          if (this->ComputeCustomMetric(simplificationMetric, a_) < simplificationThreshold)
          {
            vtkReebGraphStackPush(A_);
          }
        }