DelaunayRipsPersistenceDiagram.cpp

Go to the documentation of this file.

#include <DelaunayRipsPersistenceDiagram.h>
 
using namespace ttk::rpd;
 
ttk::DelaunayRipsPersistenceDiagram::DelaunayRipsPersistenceDiagram() {
  // inherited from Debug: prefix will be printed at the beginning of every msg
  this->setDebugMsgPrefix("DelaunayRipsPD");
}
 
int ttk::DelaunayRipsPersistenceDiagram::execute(
  const PointCloud &points, MultidimensionalDiagram &ph) const {
 
#ifdef TTK_ENABLE_CGAL
  const unsigned dim = points[0].size();
  if(points.size() <= dim) {
    printErr("Not enough points");
    return 1;
  }
  if(dim == 2) {
    FastRipsPersistenceDiagram2 FRPD(points);
    FRPD.setDebugLevel(debugLevel_);
    FRPD.computeDelaunayRips0And1Persistence(ph);
  } else if(dim == 3) {
#ifndef CGAL_LINKED_WITH_TBB
    if(getThreadNumber() > 1) {
      printWrn("TTK was not compiled with TBB:");
      printWrn("sequential Delaunay triangulation only");
    }
#endif
    gph::runDelaunayRipsPersistenceDiagram3(points, ph, getThreadNumber());
    ph[0].emplace_back(
      FiltratedSimplex{{-1}, 0.}, FiltratedSimplex{{-1}, inf}); // infinite pair
  } else {
    if(dim > TTK_DELAUNAY_MAX_COMPILED_DIMENSION)
      printWrn("High dimension: " + std::to_string(dim) + ">"
               + std::to_string(TTK_DELAUNAY_MAX_COMPILED_DIMENSION));
    gph::tryDimensions(points, ph, getThreadNumber());
    ph[0].emplace_back(
      FiltratedSimplex{{-1}, 0.}, FiltratedSimplex{{-1}, inf}); // infinite pair
  }
  return 0;
#else
  TTK_FORCE_USE(points);
  TTK_FORCE_USE(ph);
  printErr("TTK was not compiled with CGAL:");
  printErr("this filter is not available.");
  return 1;
#endif
}
 
int ttk::DelaunayRipsPersistenceDiagram::execute(
  const PointCloud &points,
  MultidimensionalDiagram &ph,
  std::vector<Generator1> &generators1,
  std::vector<Generator2> &generators2) const {
 
#ifdef TTK_ENABLE_CGAL
  const unsigned dim = points[0].size();
  if(dim > 3) {
    printErr("Input dimension too large: " + std::to_string(dim) + ">3");
    return 1;
  }
  if(points.size() <= dim) {
    printErr("Not enough points");
    return 1;
  }
  if(dim == 2) {
    FastRipsPersistenceDiagram2 FRPD(points);
    FRPD.setDebugLevel(debugLevel_);
    FRPD.computeDelaunayRips0And1Persistence(ph);
    FRPD.exportRips1Generators(generators1);
  } else if(dim == 3) {
#ifndef CGAL_LINKED_WITH_TBB
    if(getThreadNumber() > 1) {
      printWrn("TTK was not compiled with TBB:");
      printWrn("sequential Delaunay triangulation only");
    }
#endif
    gph::runDelaunayRipsPersistenceDiagram3(
      points, ph, generators1, generators2, getThreadNumber());
    ph[0].emplace_back(
      FiltratedSimplex{{-1}, 0.}, FiltratedSimplex{{-1}, inf}); // infinite pair
  }
  return 0;
#else
  TTK_FORCE_USE(points);
  TTK_FORCE_USE(ph);
  TTK_FORCE_USE(generators1);
  TTK_FORCE_USE(generators2);
  printErr("TTK was not compiled with CGAL:");
  printErr("this filter is not available.");
  return 1;
#endif
}