PersistentRipsLaplacian.hpp

Go to the documentation of this file.

#ifndef PL_R_H
#define PL_R_H
 
#include "../core/PersistentLaplacian.hpp"
#include "../../include/Ripser_modified/ripser_modified.h"
 
 
namespace PersistentLaplacians{
    template <typename eigs_Algorithm, typename up_Algorithm>
    class PersistentRipsLaplacian : public PersistentLaplacian<eigs_Algorithm, up_Algorithm> {
        public:
 
            // Lower_Distance_Matrix file format
            PersistentRipsLaplacian(const char* filename, int max_dim) : PersistentLaplacians::PersistentLaplacian<eigs_Algorithm, up_Algorithm>() { // call default constructor of parent class
                std::vector<std::vector<std::tuple<int64_t,int64_t,int>>> boundaries_triples;
                std::vector<std::vector<filtration_type>> filtrations;
                float threshold = std::numeric_limits<value_t>::max();
                rips_distance_matrix_file(filename, (long long int) max_dim, threshold, boundaries_triples, filtrations);
 
                if(boundaries_triples.size() == 0){
                    
                }
 
                std::vector<SparseMatrixInt> reindexed_boundaries(boundaries_triples.size());
                reindex_boundaries_map(boundaries_triples, reindexed_boundaries, filtrations[0].size());
                this->set_boundaries_filtrations(reindexed_boundaries, filtrations);
                
            }
 
            // Point cloud format
            PersistentRipsLaplacian(const char* filename, int max_dim, float threshold) : PersistentLaplacians::PersistentLaplacian<eigs_Algorithm, up_Algorithm>() { // call default constructor of parent class
                std::vector<std::vector<std::tuple<int64_t,int64_t,int>>> boundaries_triples;
                std::vector<std::vector<filtration_type>> filtrations;
                rips_distance_matrix_file(filename, (long long int) max_dim, threshold, boundaries_triples, filtrations);
 
                std::vector<SparseMatrixInt> reindexed_boundaries(boundaries_triples.size());
 
                reindex_boundaries_map(boundaries_triples, reindexed_boundaries, filtrations[0].size());
                this->set_boundaries_filtrations(reindexed_boundaries, filtrations);
                
            }
 
            // Point cloud format
            PersistentRipsLaplacian(std::vector<std::vector<double>> points, int max_dim){
                std::vector<std::vector<std::tuple<int64_t,int64_t,int>>> boundaries_triples;
                std::vector<std::vector<filtration_type>> filtrations;
                float threshold = std::numeric_limits<value_t>::max();
                rips_point_cloud(points, max_dim, threshold, boundaries_triples, filtrations);
                std::vector<SparseMatrixInt> reindexed_boundaries(boundaries_triples.size());
 
                reindex_boundaries_map(boundaries_triples, reindexed_boundaries, filtrations[0].size());
                this->set_boundaries_filtrations(reindexed_boundaries, filtrations);
                
            }
 
            PersistentRipsLaplacian(std::vector<std::vector<double>> points, int max_dim, float threshold){
                std::vector<std::vector<std::tuple<int64_t,int64_t,int>>> boundaries_triples;
                std::vector<std::vector<filtration_type>> filtrations;
                rips_point_cloud(points, max_dim, threshold, boundaries_triples, filtrations);
                std::vector<SparseMatrixInt> reindexed_boundaries(boundaries_triples.size());
 
                reindex_boundaries_map(boundaries_triples, reindexed_boundaries, filtrations[0].size());
                this->set_boundaries_filtrations(reindexed_boundaries, filtrations);
                
            }
        private:
 
    };
}
 
#endif