Method RequiresEdgeSubdivision

Common/DataModel/vtkSmoothErrorMetric.cxx:62–101

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

//------------------------------------------------------------------------------
int vtkSmoothErrorMetric::RequiresEdgeSubdivision(
  double* leftPoint, double* midPoint, double* rightPoint, double vtkNotUsed(alpha))
{
  assert("pre: leftPoint_exists" && leftPoint != nullptr);
  assert("pre: midPoint_exists" && midPoint != nullptr);
  assert("pre: rightPoint_exists" && rightPoint != nullptr);
  //  assert( "pre: clamped_alpha" && alpha>0 && alpha < 1. ); // or else true
  if (this->GenericCell->IsGeometryLinear())
  {
    // don't need to do anything:
    return 0;
  }

  double a[3];
  double b[3];

  a[0] = leftPoint[0] - midPoint[0];
  a[1] = leftPoint[1] - midPoint[1];
  a[2] = leftPoint[2] - midPoint[2];
  b[0] = rightPoint[0] - midPoint[0];
  b[1] = rightPoint[1] - midPoint[1];
  b[2] = rightPoint[2] - midPoint[2];

  double dota = vtkMath::Dot(a, a);
  double dotb = vtkMath::Dot(b, b);
  double cosa;

  if (dota == 0 || dotb == 0)
  {
    cosa = -1.;
  }
  else
  {
    cosa = vtkMath::Dot(a, b) / sqrt(dota * dotb);
  }

  int result = (cosa > this->CosTolerance);

  return result;
}

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