PersistentSheafLaplacian.hpp

Go to the documentation of this file.

#ifndef PSL_H
#define PSL_H
 
#include "../core/PersistentLaplacian.hpp"
#include "sheaf_simplex_tree.hpp"
 
//
//      CAUTION:
//      sheaf simplex tree uses coboundaries
//      package assumes boundaries
//      probably need to either transpose or rewrite some function for coboundaries
//
//
//
//
//
//
//
namespace PersistentLaplacians{
    template <typename eigs_Algorithm = PersistentLaplacians::selfadjoint, typename up_Algorithm = PersistentLaplacians::schur>
    class PersistentSheafLaplacian : public PersistentLaplacian<eigs_Algorithm,up_Algorithm> {
        public:
            PersistentSheafLaplacian(sheaf_simplex_tree& sst) : PersistentLaplacians::PersistentLaplacian<eigs_Algorithm>() { // call default constructor of parent class
                std::vector<FilteredBoundaryMatrix<float>> filtered_boundaries_temp = sst.apply_restriction_function();
                // this->filtered_boundaries = sst.apply_restriction_function();
                for (auto fbm : filtered_boundaries_temp){
                    this->filtered_boundaries.push_back(fbm);
                }
                this->top_dim = this->filtered_boundaries.size();
                std::cout << "Created filtered boundary matrices for cellular sheaf:" << std::endl;
                for (int i = 0; i <  (int) this->filtered_boundaries.size(); i++){
                    std::cout << "\nDim " << i << " matrix:" << std::endl;
                    this->filtered_boundaries[i].print();
                    std::cout << "Dim " << i << " domain filtration values:";
                    this->filtered_boundaries[i].print_domain_filtration();
                    std::cout << "Dim " << i << " range filtration values:";
                    this->filtered_boundaries[i].print_range_filtration();
                }
                std::cout << "got filtered coboundaries, now transpose" << std::endl;
                
                
                for (int i = 0; i < (int) this->filtered_boundaries.size(); i++){
                    this->filtered_boundaries[i] = this->filtered_boundaries[i].transpose();
                }
                std::cout << "done transposing" << std::endl;
            }
        private:
 
 
    };
}
 
#endif