Swapping Dionysus for Dionysus2

R/RcppExports.R

@@ -0,0 +1,67 @@
+# Generated by using Rcpp::compileAttributes() -> do not edit by hand
+# Generator token: 10BE3573-1514-4C36-9D1C-5A225CD40393
+
+GridFiltration <- function(FUNvalues, gridDim, maxdimension, decomposition, printProgress) {
+    .Call(`_TDA_GridFiltration`, FUNvalues, gridDim, maxdimension, decomposition, printProgress)
+}
+
+GridDiag <- function(FUNvalues, gridDim, maxdimension, decomposition, library, location, printProgress) {
+    .Call(`_TDA_GridDiag`, FUNvalues, gridDim, maxdimension, decomposition, library, location, printProgress)
+}
+
+Bottleneck <- function(Diag1, Diag2) {
+    .Call(`_TDA_Bottleneck`, Diag1, Diag2)
+}
+
+Wasserstein <- function(Diag1, Diag2, p) {
+    .Call(`_TDA_Wasserstein`, Diag1, Diag2, p)
+}
+
+Kde <- function(X, Grid, h, kertype, weight, printProgress) {
+    .Call(`_TDA_Kde`, X, Grid, h, kertype, weight, printProgress)
+}
+
+KdeDist <- function(X, Grid, h, weight, printProgress) {
+    .Call(`_TDA_KdeDist`, X, Grid, h, weight, printProgress)
+}
+
+Dtm <- function(knnDistance, weightBound, r) {
+    .Call(`_TDA_Dtm`, knnDistance, weightBound, r)
+}
+
+DtmWeight <- function(knnDistance, weightBound, r, knnIndex, weight) {
+    .Call(`_TDA_DtmWeight`, knnDistance, weightBound, r, knnIndex, weight)
+}
+
+FiltrationDiag <- function(filtration, maxdimension, library, location, printProgress) {
+    .Call(`_TDA_FiltrationDiag`, filtration, maxdimension, library, location, printProgress)
+}
+
+FunFiltration <- function(FUNvalues, cmplx) {
+    .Call(`_TDA_FunFiltration`, FUNvalues, cmplx)
+}
+
+RipsFiltration <- function(X, maxdimension, maxscale, dist, library, printProgress) {
+    .Call(`_TDA_RipsFiltration`, X, maxdimension, maxscale, dist, library, printProgress)
+}
+
+RipsDiag <- function(X, maxdimension, maxscale, dist, libraryFiltration, libraryDiag, location, printProgress) {
+    .Call(`_TDA_RipsDiag`, X, maxdimension, maxscale, dist, libraryFiltration, libraryDiag, location, printProgress)
+}
+
+AlphaShapeFiltration <- function(X, printProgress) {
+    .Call(`_TDA_AlphaShapeFiltration`, X, printProgress)
+}
+
+AlphaShapeDiag <- function(X, maxdimension, libraryDiag, location, printProgress) {
+    .Call(`_TDA_AlphaShapeDiag`, X, maxdimension, libraryDiag, location, printProgress)
+}
+
+AlphaComplexFiltration <- function(X, printProgress) {
+    .Call(`_TDA_AlphaComplexFiltration`, X, printProgress)
+}
+
+AlphaComplexDiag <- function(X, maxdimension, libraryDiag, location, printProgress) {
+    .Call(`_TDA_AlphaComplexDiag`, X, maxdimension, libraryDiag, location, printProgress)
+}
+

R/filtrationDiag.R

@@ -15,8 +15,11 @@ filtrationDiag <- function(
   if (library == "dionysus" || library == "DIONYSUS") {
     library <- "Dionysus"
   }
-  if (library != "GUDHI" && library != "Dionysus") {
-    stop("library for computing persistence diagram should be a string: either 'GUDHI' or 'Dionysus'")
+  if (library == "D2") {
+      library <- "D2"
+  }
+  if (library != "GUDHI" && library != "Dionysus" && library != "D2") {
+    stop("library for computing persistence diagram should be a string: either 'GUDHI' or 'Dionysus' or 'Dionysus2'")
   }
   if (!is.logical(location)) {
     stop("location should be logical")
@@ -71,4 +74,4 @@ filtrationDiag <- function(
 		"deathLocation" = DeathLocation, "cycleLocation" = CycleLocation)
   }
   return (out)
-}
+}

R/gridDiag.R

@@ -127,6 +127,7 @@ function(X = NULL, FUN = NULL, lim = NULL, by = NULL, FUNvalues = NULL,
     Diag[1, 3] <- ifelse(is.null(diagLimit), max(FUNvalues), diagLimit) 
   }
   if (sublevel == FALSE) {
+    #possible bugfix with labels
     colnames(Diag) <- c("dimension", "Death", "Birth")
     Diag[, 2:3] <- -Diag[, 3:2]
   } else {

R/ripsFiltration.R

@@ -21,7 +21,10 @@ ripsFiltration <- function(
   if (library == "dionysus" || library == "DIONYSUS") {
     library <- "Dionysus"
   }
-  if (library != "GUDHI" && library != "Dionysus") {
+  if (library == "D2") {
+      library <- "D2"
+  }
+  if (library != "GUDHI" && library != "Dionysus" && library != "D2") {
     stop("library should be a string: either 'GUDHI' or 'Dionysus'")
   }
   if (!is.logical(printProgress)) {

src/RcppExports.cpp

@@ -0,0 +1,264 @@
+// Generated by using Rcpp::compileAttributes() -> do not edit by hand
+// Generator token: 10BE3573-1514-4C36-9D1C-5A225CD40393
+
+#include <RcppEigen.h>
+#include <Rcpp.h>
+
+using namespace Rcpp;
+
+// GridFiltration
+Rcpp::List GridFiltration(const Rcpp::NumericVector& FUNvalues, const Rcpp::IntegerVector& gridDim, const int maxdimension, const std::string& decomposition, const bool printProgress);
+RcppExport SEXP _TDA_GridFiltration(SEXP FUNvaluesSEXP, SEXP gridDimSEXP, SEXP maxdimensionSEXP, SEXP decompositionSEXP, SEXP printProgressSEXP) {
+BEGIN_RCPP
+    Rcpp::RObject rcpp_result_gen;
+    Rcpp::RNGScope rcpp_rngScope_gen;
+    Rcpp::traits::input_parameter< const Rcpp::NumericVector& >::type FUNvalues(FUNvaluesSEXP);
+    Rcpp::traits::input_parameter< const Rcpp::IntegerVector& >::type gridDim(gridDimSEXP);
+    Rcpp::traits::input_parameter< const int >::type maxdimension(maxdimensionSEXP);
+    Rcpp::traits::input_parameter< const std::string& >::type decomposition(decompositionSEXP);
+    Rcpp::traits::input_parameter< const bool >::type printProgress(printProgressSEXP);
+    rcpp_result_gen = Rcpp::wrap(GridFiltration(FUNvalues, gridDim, maxdimension, decomposition, printProgress));
+    return rcpp_result_gen;
+END_RCPP
+}
+// GridDiag
+Rcpp::List GridDiag(const Rcpp::NumericVector& FUNvalues, const Rcpp::IntegerVector& gridDim, const int maxdimension, const std::string& decomposition, const std::string& library, const bool location, const bool printProgress);
+RcppExport SEXP _TDA_GridDiag(SEXP FUNvaluesSEXP, SEXP gridDimSEXP, SEXP maxdimensionSEXP, SEXP decompositionSEXP, SEXP librarySEXP, SEXP locationSEXP, SEXP printProgressSEXP) {
+BEGIN_RCPP
+    Rcpp::RObject rcpp_result_gen;
+    Rcpp::RNGScope rcpp_rngScope_gen;
+    Rcpp::traits::input_parameter< const Rcpp::NumericVector& >::type FUNvalues(FUNvaluesSEXP);
+    Rcpp::traits::input_parameter< const Rcpp::IntegerVector& >::type gridDim(gridDimSEXP);
+    Rcpp::traits::input_parameter< const int >::type maxdimension(maxdimensionSEXP);
+    Rcpp::traits::input_parameter< const std::string& >::type decomposition(decompositionSEXP);
+    Rcpp::traits::input_parameter< const std::string& >::type library(librarySEXP);
+    Rcpp::traits::input_parameter< const bool >::type location(locationSEXP);
+    Rcpp::traits::input_parameter< const bool >::type printProgress(printProgressSEXP);
+    rcpp_result_gen = Rcpp::wrap(GridDiag(FUNvalues, gridDim, maxdimension, decomposition, library, location, printProgress));
+    return rcpp_result_gen;
+END_RCPP
+}
+// Bottleneck
+double Bottleneck(const Rcpp::NumericMatrix& Diag1, const Rcpp::NumericMatrix& Diag2);
+RcppExport SEXP _TDA_Bottleneck(SEXP Diag1SEXP, SEXP Diag2SEXP) {
+BEGIN_RCPP
+    Rcpp::RObject rcpp_result_gen;
+    Rcpp::RNGScope rcpp_rngScope_gen;
+    Rcpp::traits::input_parameter< const Rcpp::NumericMatrix& >::type Diag1(Diag1SEXP);
+    Rcpp::traits::input_parameter< const Rcpp::NumericMatrix& >::type Diag2(Diag2SEXP);
+    rcpp_result_gen = Rcpp::wrap(Bottleneck(Diag1, Diag2));
+    return rcpp_result_gen;
+END_RCPP
+}
+// Wasserstein
+double Wasserstein(const Rcpp::NumericMatrix& Diag1, const Rcpp::NumericMatrix& Diag2, const int p);
+RcppExport SEXP _TDA_Wasserstein(SEXP Diag1SEXP, SEXP Diag2SEXP, SEXP pSEXP) {
+BEGIN_RCPP
+    Rcpp::RObject rcpp_result_gen;
+    Rcpp::RNGScope rcpp_rngScope_gen;
+    Rcpp::traits::input_parameter< const Rcpp::NumericMatrix& >::type Diag1(Diag1SEXP);
+    Rcpp::traits::input_parameter< const Rcpp::NumericMatrix& >::type Diag2(Diag2SEXP);
+    Rcpp::traits::input_parameter< const int >::type p(pSEXP);
+    rcpp_result_gen = Rcpp::wrap(Wasserstein(Diag1, Diag2, p));
+    return rcpp_result_gen;
+END_RCPP
+}
+// Kde
+Rcpp::NumericVector Kde(const Rcpp::NumericMatrix& X, const Rcpp::NumericMatrix& Grid, const double h, const std::string& kertype, const Rcpp::NumericVector& weight, const bool printProgress);
+RcppExport SEXP _TDA_Kde(SEXP XSEXP, SEXP GridSEXP, SEXP hSEXP, SEXP kertypeSEXP, SEXP weightSEXP, SEXP printProgressSEXP) {
+BEGIN_RCPP
+    Rcpp::RObject rcpp_result_gen;
+    Rcpp::RNGScope rcpp_rngScope_gen;
+    Rcpp::traits::input_parameter< const Rcpp::NumericMatrix& >::type X(XSEXP);
+    Rcpp::traits::input_parameter< const Rcpp::NumericMatrix& >::type Grid(GridSEXP);
+    Rcpp::traits::input_parameter< const double >::type h(hSEXP);
+    Rcpp::traits::input_parameter< const std::string& >::type kertype(kertypeSEXP);
+    Rcpp::traits::input_parameter< const Rcpp::NumericVector& >::type weight(weightSEXP);
+    Rcpp::traits::input_parameter< const bool >::type printProgress(printProgressSEXP);
+    rcpp_result_gen = Rcpp::wrap(Kde(X, Grid, h, kertype, weight, printProgress));
+    return rcpp_result_gen;
+END_RCPP
+}
+// KdeDist
+Rcpp::NumericVector KdeDist(const Rcpp::NumericMatrix& X, const Rcpp::NumericMatrix& Grid, const double h, const Rcpp::NumericVector& weight, const bool printProgress);
+RcppExport SEXP _TDA_KdeDist(SEXP XSEXP, SEXP GridSEXP, SEXP hSEXP, SEXP weightSEXP, SEXP printProgressSEXP) {
+BEGIN_RCPP
+    Rcpp::RObject rcpp_result_gen;
+    Rcpp::RNGScope rcpp_rngScope_gen;
+    Rcpp::traits::input_parameter< const Rcpp::NumericMatrix& >::type X(XSEXP);
+    Rcpp::traits::input_parameter< const Rcpp::NumericMatrix& >::type Grid(GridSEXP);
+    Rcpp::traits::input_parameter< const double >::type h(hSEXP);
+    Rcpp::traits::input_parameter< const Rcpp::NumericVector& >::type weight(weightSEXP);
+    Rcpp::traits::input_parameter< const bool >::type printProgress(printProgressSEXP);
+    rcpp_result_gen = Rcpp::wrap(KdeDist(X, Grid, h, weight, printProgress));
+    return rcpp_result_gen;
+END_RCPP
+}
+// Dtm
+Rcpp::NumericVector Dtm(const Rcpp::NumericMatrix& knnDistance, const double weightBound, const double r);
+RcppExport SEXP _TDA_Dtm(SEXP knnDistanceSEXP, SEXP weightBoundSEXP, SEXP rSEXP) {
+BEGIN_RCPP
+    Rcpp::RObject rcpp_result_gen;
+    Rcpp::RNGScope rcpp_rngScope_gen;
+    Rcpp::traits::input_parameter< const Rcpp::NumericMatrix& >::type knnDistance(knnDistanceSEXP);
+    Rcpp::traits::input_parameter< const double >::type weightBound(weightBoundSEXP);
+    Rcpp::traits::input_parameter< const double >::type r(rSEXP);
+    rcpp_result_gen = Rcpp::wrap(Dtm(knnDistance, weightBound, r));
+    return rcpp_result_gen;
+END_RCPP
+}
+// DtmWeight
+Rcpp::NumericVector DtmWeight(const Rcpp::NumericMatrix& knnDistance, const double weightBound, const double r, const Rcpp::NumericMatrix& knnIndex, const Rcpp::NumericVector& weight);
+RcppExport SEXP _TDA_DtmWeight(SEXP knnDistanceSEXP, SEXP weightBoundSEXP, SEXP rSEXP, SEXP knnIndexSEXP, SEXP weightSEXP) {
+BEGIN_RCPP
+    Rcpp::RObject rcpp_result_gen;
+    Rcpp::RNGScope rcpp_rngScope_gen;
+    Rcpp::traits::input_parameter< const Rcpp::NumericMatrix& >::type knnDistance(knnDistanceSEXP);
+    Rcpp::traits::input_parameter< const double >::type weightBound(weightBoundSEXP);
+    Rcpp::traits::input_parameter< const double >::type r(rSEXP);
+    Rcpp::traits::input_parameter< const Rcpp::NumericMatrix& >::type knnIndex(knnIndexSEXP);
+    Rcpp::traits::input_parameter< const Rcpp::NumericVector& >::type weight(weightSEXP);
+    rcpp_result_gen = Rcpp::wrap(DtmWeight(knnDistance, weightBound, r, knnIndex, weight));
+    return rcpp_result_gen;
+END_RCPP
+}
+// FiltrationDiag
+Rcpp::List FiltrationDiag(const Rcpp::List& filtration, const int maxdimension, const std::string& library, const bool location, const bool printProgress);
+RcppExport SEXP _TDA_FiltrationDiag(SEXP filtrationSEXP, SEXP maxdimensionSEXP, SEXP librarySEXP, SEXP locationSEXP, SEXP printProgressSEXP) {
+BEGIN_RCPP
+    Rcpp::RObject rcpp_result_gen;
+    Rcpp::RNGScope rcpp_rngScope_gen;
+    Rcpp::traits::input_parameter< const Rcpp::List& >::type filtration(filtrationSEXP);
+    Rcpp::traits::input_parameter< const int >::type maxdimension(maxdimensionSEXP);
+    Rcpp::traits::input_parameter< const std::string& >::type library(librarySEXP);
+    Rcpp::traits::input_parameter< const bool >::type location(locationSEXP);
+    Rcpp::traits::input_parameter< const bool >::type printProgress(printProgressSEXP);
+    rcpp_result_gen = Rcpp::wrap(FiltrationDiag(filtration, maxdimension, library, location, printProgress));
+    return rcpp_result_gen;
+END_RCPP
+}
+// FunFiltration
+Rcpp::List FunFiltration(const Rcpp::NumericVector& FUNvalues, const Rcpp::List& cmplx);
+RcppExport SEXP _TDA_FunFiltration(SEXP FUNvaluesSEXP, SEXP cmplxSEXP) {
+BEGIN_RCPP
+    Rcpp::RObject rcpp_result_gen;
+    Rcpp::RNGScope rcpp_rngScope_gen;
+    Rcpp::traits::input_parameter< const Rcpp::NumericVector& >::type FUNvalues(FUNvaluesSEXP);
+    Rcpp::traits::input_parameter< const Rcpp::List& >::type cmplx(cmplxSEXP);
+    rcpp_result_gen = Rcpp::wrap(FunFiltration(FUNvalues, cmplx));
+    return rcpp_result_gen;
+END_RCPP
+}
+// RipsFiltration
+Rcpp::List RipsFiltration(const Rcpp::NumericMatrix& X, const int maxdimension, const double maxscale, const std::string& dist, const std::string& library, const bool printProgress);
+RcppExport SEXP _TDA_RipsFiltration(SEXP XSEXP, SEXP maxdimensionSEXP, SEXP maxscaleSEXP, SEXP distSEXP, SEXP librarySEXP, SEXP printProgressSEXP) {
+BEGIN_RCPP
+    Rcpp::RObject rcpp_result_gen;
+    Rcpp::RNGScope rcpp_rngScope_gen;
+    Rcpp::traits::input_parameter< const Rcpp::NumericMatrix& >::type X(XSEXP);
+    Rcpp::traits::input_parameter< const int >::type maxdimension(maxdimensionSEXP);
+    Rcpp::traits::input_parameter< const double >::type maxscale(maxscaleSEXP);
+    Rcpp::traits::input_parameter< const std::string& >::type dist(distSEXP);
+    Rcpp::traits::input_parameter< const std::string& >::type library(librarySEXP);
+    Rcpp::traits::input_parameter< const bool >::type printProgress(printProgressSEXP);
+    rcpp_result_gen = Rcpp::wrap(RipsFiltration(X, maxdimension, maxscale, dist, library, printProgress));
+    return rcpp_result_gen;
+END_RCPP
+}
+// RipsDiag
+Rcpp::List RipsDiag(const Rcpp::NumericMatrix& X, const int maxdimension, const double maxscale, const std::string& dist, const std::string& libraryFiltration, const std::string& libraryDiag, const bool location, const bool printProgress);
+RcppExport SEXP _TDA_RipsDiag(SEXP XSEXP, SEXP maxdimensionSEXP, SEXP maxscaleSEXP, SEXP distSEXP, SEXP libraryFiltrationSEXP, SEXP libraryDiagSEXP, SEXP locationSEXP, SEXP printProgressSEXP) {
+BEGIN_RCPP
+    Rcpp::RObject rcpp_result_gen;
+    Rcpp::RNGScope rcpp_rngScope_gen;
+    Rcpp::traits::input_parameter< const Rcpp::NumericMatrix& >::type X(XSEXP);
+    Rcpp::traits::input_parameter< const int >::type maxdimension(maxdimensionSEXP);
+    Rcpp::traits::input_parameter< const double >::type maxscale(maxscaleSEXP);
+    Rcpp::traits::input_parameter< const std::string& >::type dist(distSEXP);
+    Rcpp::traits::input_parameter< const std::string& >::type libraryFiltration(libraryFiltrationSEXP);
+    Rcpp::traits::input_parameter< const std::string& >::type libraryDiag(libraryDiagSEXP);
+    Rcpp::traits::input_parameter< const bool >::type location(locationSEXP);
+    Rcpp::traits::input_parameter< const bool >::type printProgress(printProgressSEXP);
+    rcpp_result_gen = Rcpp::wrap(RipsDiag(X, maxdimension, maxscale, dist, libraryFiltration, libraryDiag, location, printProgress));
+    return rcpp_result_gen;
+END_RCPP
+}
+// AlphaShapeFiltration
+Rcpp::List AlphaShapeFiltration(const Rcpp::NumericMatrix& X, const bool printProgress);
+RcppExport SEXP _TDA_AlphaShapeFiltration(SEXP XSEXP, SEXP printProgressSEXP) {
+BEGIN_RCPP
+    Rcpp::RObject rcpp_result_gen;
+    Rcpp::RNGScope rcpp_rngScope_gen;
+    Rcpp::traits::input_parameter< const Rcpp::NumericMatrix& >::type X(XSEXP);
+    Rcpp::traits::input_parameter< const bool >::type printProgress(printProgressSEXP);
+    rcpp_result_gen = Rcpp::wrap(AlphaShapeFiltration(X, printProgress));
+    return rcpp_result_gen;
+END_RCPP
+}
+// AlphaShapeDiag
+Rcpp::List AlphaShapeDiag(const Rcpp::NumericMatrix& X, const int maxdimension, const std::string& libraryDiag, const bool location, const bool printProgress);
+RcppExport SEXP _TDA_AlphaShapeDiag(SEXP XSEXP, SEXP maxdimensionSEXP, SEXP libraryDiagSEXP, SEXP locationSEXP, SEXP printProgressSEXP) {
+BEGIN_RCPP
+    Rcpp::RObject rcpp_result_gen;
+    Rcpp::RNGScope rcpp_rngScope_gen;
+    Rcpp::traits::input_parameter< const Rcpp::NumericMatrix& >::type X(XSEXP);
+    Rcpp::traits::input_parameter< const int >::type maxdimension(maxdimensionSEXP);
+    Rcpp::traits::input_parameter< const std::string& >::type libraryDiag(libraryDiagSEXP);
+    Rcpp::traits::input_parameter< const bool >::type location(locationSEXP);
+    Rcpp::traits::input_parameter< const bool >::type printProgress(printProgressSEXP);
+    rcpp_result_gen = Rcpp::wrap(AlphaShapeDiag(X, maxdimension, libraryDiag, location, printProgress));
+    return rcpp_result_gen;
+END_RCPP
+}
+// AlphaComplexFiltration
+Rcpp::List AlphaComplexFiltration(const Rcpp::NumericMatrix& X, const bool printProgress);
+RcppExport SEXP _TDA_AlphaComplexFiltration(SEXP XSEXP, SEXP printProgressSEXP) {
+BEGIN_RCPP
+    Rcpp::RObject rcpp_result_gen;
+    Rcpp::RNGScope rcpp_rngScope_gen;
+    Rcpp::traits::input_parameter< const Rcpp::NumericMatrix& >::type X(XSEXP);
+    Rcpp::traits::input_parameter< const bool >::type printProgress(printProgressSEXP);
+    rcpp_result_gen = Rcpp::wrap(AlphaComplexFiltration(X, printProgress));
+    return rcpp_result_gen;
+END_RCPP
+}
+// AlphaComplexDiag
+Rcpp::List AlphaComplexDiag(const Rcpp::NumericMatrix& X, const int maxdimension, const std::string& libraryDiag, const bool location, const bool printProgress);
+RcppExport SEXP _TDA_AlphaComplexDiag(SEXP XSEXP, SEXP maxdimensionSEXP, SEXP libraryDiagSEXP, SEXP locationSEXP, SEXP printProgressSEXP) {
+BEGIN_RCPP
+    Rcpp::RObject rcpp_result_gen;
+    Rcpp::RNGScope rcpp_rngScope_gen;
+    Rcpp::traits::input_parameter< const Rcpp::NumericMatrix& >::type X(XSEXP);
+    Rcpp::traits::input_parameter< const int >::type maxdimension(maxdimensionSEXP);
+    Rcpp::traits::input_parameter< const std::string& >::type libraryDiag(libraryDiagSEXP);
+    Rcpp::traits::input_parameter< const bool >::type location(locationSEXP);
+    Rcpp::traits::input_parameter< const bool >::type printProgress(printProgressSEXP);
+    rcpp_result_gen = Rcpp::wrap(AlphaComplexDiag(X, maxdimension, libraryDiag, location, printProgress));
+    return rcpp_result_gen;
+END_RCPP
+}
+
+static const R_CallMethodDef CallEntries[] = {
+    {"_TDA_GridFiltration", (DL_FUNC) &_TDA_GridFiltration, 5},
+    {"_TDA_GridDiag", (DL_FUNC) &_TDA_GridDiag, 7},
+    {"_TDA_Bottleneck", (DL_FUNC) &_TDA_Bottleneck, 2},
+    {"_TDA_Wasserstein", (DL_FUNC) &_TDA_Wasserstein, 3},
+    {"_TDA_Kde", (DL_FUNC) &_TDA_Kde, 6},
+    {"_TDA_KdeDist", (DL_FUNC) &_TDA_KdeDist, 5},
+    {"_TDA_Dtm", (DL_FUNC) &_TDA_Dtm, 3},
+    {"_TDA_DtmWeight", (DL_FUNC) &_TDA_DtmWeight, 5},
+    {"_TDA_FiltrationDiag", (DL_FUNC) &_TDA_FiltrationDiag, 5},
+    {"_TDA_FunFiltration", (DL_FUNC) &_TDA_FunFiltration, 2},
+    {"_TDA_RipsFiltration", (DL_FUNC) &_TDA_RipsFiltration, 6},
+    {"_TDA_RipsDiag", (DL_FUNC) &_TDA_RipsDiag, 8},
+    {"_TDA_AlphaShapeFiltration", (DL_FUNC) &_TDA_AlphaShapeFiltration, 2},
+    {"_TDA_AlphaShapeDiag", (DL_FUNC) &_TDA_AlphaShapeDiag, 5},
+    {"_TDA_AlphaComplexFiltration", (DL_FUNC) &_TDA_AlphaComplexFiltration, 2},
+    {"_TDA_AlphaComplexDiag", (DL_FUNC) &_TDA_AlphaComplexDiag, 5},
+    {NULL, NULL, 0}
+};
+
+RcppExport void R_init_TDA(DllInfo *dll) {
+    R_registerRoutines(dll, NULL, CallEntries, NULL, NULL);
+    R_useDynamicSymbols(dll, FALSE);
+}