#include <gudhi/Bitmap_cubical_complex.h>
#include <gudhi/Persistent_cohomology.h>
#include <iostream>
#include <string>
#include <vector>
#include <cstddef>
int main(int argc, char** argv) {
std::cout << "This program computes persistent homology, by using bitmap_cubical_complex class, of cubical " <<
"complexes provided in text files in Perseus style (the only numbered in the first line is a dimension D of a" <<
"bitmap. In the lines I between 2 and D+1 there are numbers of top dimensional cells in the direction I. Let " <<
"N denote product of the numbers in the lines between 2 and D. In the lines D+2 to D+2+N there are " <<
"filtrations of top dimensional cells. We assume that the cells are in the lexicographical order. See " <<
"CubicalOneSphere.txt or CubicalTwoSphere.txt for example.\n" << std::endl;
if (argc != 2) {
std::cerr << "Wrong number of parameters. Please provide the name of a file with a Perseus style bitmap at " <<
"the input. The program will now terminate.\n";
return 1;
}
Bitmap_cubical_complex b(argv[1]);
Persistent_cohomology pcoh(b);
int p = 2;
double min_persistence = 0;
pcoh.init_coefficients(p);
pcoh.compute_persistent_cohomology(min_persistence);
std::string output_file_name(argv[1]);
output_file_name += "_persistence";
std::size_t last_in_path = output_file_name.find_last_of("/\\");
if (last_in_path != std::string::npos) {
output_file_name = output_file_name.substr(last_in_path+1);
}
std::ofstream out(output_file_name.c_str());
pcoh.output_diagram(out);
out.close();
std::cout << "Result in file: " << output_file_name << "\n";
return 0;
}