summaryrefslogtreecommitdiff
path: root/src/Persistence_representations/example/persistence_heat_maps.cpp
diff options
context:
space:
mode:
authorvrouvrea <vrouvrea@636b058d-ea47-450e-bf9e-a15bfbe3eedb>2017-05-30 15:52:00 +0000
committervrouvrea <vrouvrea@636b058d-ea47-450e-bf9e-a15bfbe3eedb>2017-05-30 15:52:00 +0000
commit9d1a526de85694b5f075bb88dbd7097a40abf10a (patch)
treebbcd0cef32610d2f5e9c0209b48c58f73fbf379a /src/Persistence_representations/example/persistence_heat_maps.cpp
parent2bcb3d7cb47ce71803f2464cc822346ed2e1b039 (diff)
clang format all sources
git-svn-id: svn+ssh://scm.gforge.inria.fr/svnroot/gudhi/branches/persistence_representation_integration@2477 636b058d-ea47-450e-bf9e-a15bfbe3eedb Former-commit-id: 326d664483d6700f82be824f79a0bf5c082b4945
Diffstat (limited to 'src/Persistence_representations/example/persistence_heat_maps.cpp')
-rw-r--r--src/Persistence_representations/example/persistence_heat_maps.cpp99
1 files changed, 44 insertions, 55 deletions
diff --git a/src/Persistence_representations/example/persistence_heat_maps.cpp b/src/Persistence_representations/example/persistence_heat_maps.cpp
index c75e2731..44227823 100644
--- a/src/Persistence_representations/example/persistence_heat_maps.cpp
+++ b/src/Persistence_representations/example/persistence_heat_maps.cpp
@@ -20,74 +20,63 @@
* along with this program. If not, see <http://www.gnu.org/licenses/>.
*/
-
-
#include <gudhi/reader_utils.h>
#include <gudhi/Persistence_heat_maps.h>
#include <iostream>
#include <vector>
-
-
using namespace Gudhi;
using namespace Gudhi::Persistence_representations;
-
double epsilon = 0.0000005;
+int main(int argc, char** argv) {
+ // create two simple vectors with birth--death pairs:
-
-int main( int argc , char** argv )
-{
- //create two simple vectors with birth--death pairs:
-
- std::vector< std::pair< double , double > > persistence1;
- std::vector< std::pair< double , double > > persistence2;
-
- persistence1.push_back( std::make_pair(1,2) );
- persistence1.push_back( std::make_pair(6,8) );
- persistence1.push_back( std::make_pair(0,4) );
- persistence1.push_back( std::make_pair(3,8) );
-
- persistence2.push_back( std::make_pair(2,9) );
- persistence2.push_back( std::make_pair(1,6) );
- persistence2.push_back( std::make_pair(3,5) );
- persistence2.push_back( std::make_pair(6,10) );
-
- //over here we define a function we sill put on a top on every birth--death pair in the persistence interval. It can be anything. Over here we will use standard Gaussian
- std::vector< std::vector<double> > filter = create_Gaussian_filter(5,1);
-
- //creating two heat maps.
- Persistence_heat_maps<constant_scaling_function> hm1( persistence1 , filter , false , 20 , 0 , 11 );
- Persistence_heat_maps<constant_scaling_function> hm2( persistence2 , filter , false , 20 , 0 , 11 );
-
- std::vector<Persistence_heat_maps<constant_scaling_function>*> vector_of_maps;
- vector_of_maps.push_back( &hm1 );
- vector_of_maps.push_back( &hm2 );
-
- //compute median/mean of a vector of heat maps:
- Persistence_heat_maps<constant_scaling_function> mean;
- mean.compute_mean( vector_of_maps );
- Persistence_heat_maps<constant_scaling_function> median;
- median.compute_median( vector_of_maps );
-
- //to compute L^1 distance between hm1 and hm2:
- std::cout << "The L^1 distance is : " << hm1.distance( hm2 , 1 ) << std::endl;
-
- //to average of hm1 and hm2:
- std::vector< Persistence_heat_maps<constant_scaling_function>* > to_average;
- to_average.push_back( &hm1 );
- to_average.push_back( &hm2 );
- Persistence_heat_maps<constant_scaling_function> av;
- av.compute_average( to_average );
-
- //to compute scalar product of hm1 and hm2:
- std::cout << "Scalar product is : " << hm1.compute_scalar_product( hm2 ) << std::endl;
-
- return 0;
-}
+ std::vector<std::pair<double, double> > persistence1;
+ std::vector<std::pair<double, double> > persistence2;
+
+ persistence1.push_back(std::make_pair(1, 2));
+ persistence1.push_back(std::make_pair(6, 8));
+ persistence1.push_back(std::make_pair(0, 4));
+ persistence1.push_back(std::make_pair(3, 8));
+ persistence2.push_back(std::make_pair(2, 9));
+ persistence2.push_back(std::make_pair(1, 6));
+ persistence2.push_back(std::make_pair(3, 5));
+ persistence2.push_back(std::make_pair(6, 10));
+ // over here we define a function we sill put on a top on every birth--death pair in the persistence interval. It can
+ // be anything. Over here we will use standard Gaussian
+ std::vector<std::vector<double> > filter = create_Gaussian_filter(5, 1);
+ // creating two heat maps.
+ Persistence_heat_maps<constant_scaling_function> hm1(persistence1, filter, false, 20, 0, 11);
+ Persistence_heat_maps<constant_scaling_function> hm2(persistence2, filter, false, 20, 0, 11);
+ std::vector<Persistence_heat_maps<constant_scaling_function>*> vector_of_maps;
+ vector_of_maps.push_back(&hm1);
+ vector_of_maps.push_back(&hm2);
+
+ // compute median/mean of a vector of heat maps:
+ Persistence_heat_maps<constant_scaling_function> mean;
+ mean.compute_mean(vector_of_maps);
+ Persistence_heat_maps<constant_scaling_function> median;
+ median.compute_median(vector_of_maps);
+
+ // to compute L^1 distance between hm1 and hm2:
+ std::cout << "The L^1 distance is : " << hm1.distance(hm2, 1) << std::endl;
+
+ // to average of hm1 and hm2:
+ std::vector<Persistence_heat_maps<constant_scaling_function>*> to_average;
+ to_average.push_back(&hm1);
+ to_average.push_back(&hm2);
+ Persistence_heat_maps<constant_scaling_function> av;
+ av.compute_average(to_average);
+
+ // to compute scalar product of hm1 and hm2:
+ std::cout << "Scalar product is : " << hm1.compute_scalar_product(hm2) << std::endl;
+
+ return 0;
+}