RipsPersistenceGenerators.h

Go to the documentation of this file.

 
#pragma once
 
// ttk common includes
#include <Debug.h>
 
#include <PairCellsWithOracle.h>
 
namespace ttk {
 
  class RipsPersistenceGenerators : virtual public Debug {
  public:
    RipsPersistenceGenerators();
    void execute(const std::vector<std::vector<double>> &points,
                 rpd::MultidimensionalDiagram &diagrams,
                 std::vector<rpd::Generator1> &generators) const {
      rpd::PairCellsWithOracle::callOracle(
        points, diagrams, SimplexMaximumDiameter, InputIsDistanceMatrix);
      rpd::PairCellsWithOracle pc(
        points, diagrams, InputIsDistanceMatrix, false);
      pc.setDebugLevel(debugLevel_);
      pc.run();
      if(!OutputCascade)
        pc.getGenerators(generators);
      else {
        rpd::EdgeSets4 criticalAndCascade;
        pc.getCascades(criticalAndCascade);
        generators.emplace_back(
          criticalAndCascade[rpd::DEATH0], std::make_pair(0., 0.));
        generators.emplace_back(
          criticalAndCascade[rpd::BIRTH1], std::make_pair(1., 1.));
        generators.emplace_back(
          criticalAndCascade[rpd::DEATH1], std::make_pair(2., 2.));
        generators.emplace_back(
          criticalAndCascade[rpd::CASC1], std::make_pair(3., 3.));
      }
    }
 
  protected:
    double SimplexMaximumDiameter{rpd::inf};
    bool InputIsDistanceMatrix{false};
    bool OutputCascade{false};
  };
 
} // namespace ttk