alpha shape creation update

git-svn-id: https://phat.googlecode.com/svn/trunk@50 8e3bb3c2-eed4-f18f-5264-0b6c94e6926d

1 files changed, 16 insertions, 0 deletions

diff --git a/include/phat/compute_persistence_pairs.h b/include/phat/compute_persistence_pairs.h
index e782bd7..ca46ac0 100644
--- a/include/phat/compute_persistence_pairs.h
+++ b/include/phat/compute_persistence_pairs.h
@@ -43,8 +43,24 @@ namespace phat {
     void compute_persistence_pairs_dualized( persistence_pairs& pairs, boundary_matrix< Representation >& boundary_matrix ) {
         ReductionAlgorithm reduce;
         const index nr_columns = boundary_matrix.get_num_cols();
+        
+        double start;
+        
+        start= omp_get_wtime();
+        std::clog << "Dualize ... " << std::endl;
+
         dualize( boundary_matrix );
+        
+        std::clog <<  omp_get_wtime() - start <<"s" << std::endl;
+
+        start= omp_get_wtime();
+        std::clog << "Reduce ... " << std::endl;
+        
         reduce( boundary_matrix );
+        
+        std::clog <<  omp_get_wtime() - start <<"s" << std::endl;
+
+        
         pairs.clear();
         for( index idx = 0; idx < nr_columns; idx++ ) {
             if( !boundary_matrix.is_empty( idx ) ) {