From 5cec22844dedc8f27ed0c2d880ab3b0431cb1a88 Mon Sep 17 00:00:00 2001 From: pdlotko Date: Wed, 1 Nov 2017 14:24:38 +0000 Subject: Bug fixed. git-svn-id: svn+ssh://scm.gforge.inria.fr/svnroot/gudhi/branches/periodic_cubical_complex_fix@2821 636b058d-ea47-450e-bf9e-a15bfbe3eedb Former-commit-id: f5816f9d7ca51a30c5f64dc7cd58f40ffc157917 --- src/Bitmap_cubical_complex/test/Bitmap_test.cpp | 378 +++++++++++++++++++++--- 1 file changed, 337 insertions(+), 41 deletions(-) (limited to 'src/Bitmap_cubical_complex/test') diff --git a/src/Bitmap_cubical_complex/test/Bitmap_test.cpp b/src/Bitmap_cubical_complex/test/Bitmap_test.cpp index db90eb94..f815a70c 100644 --- a/src/Bitmap_cubical_complex/test/Bitmap_test.cpp +++ b/src/Bitmap_cubical_complex/test/Bitmap_test.cpp @@ -42,6 +42,8 @@ Bitmap_cubical_complex_periodic_boundary_conditions_base; typedef Gudhi::cubical_complex::Bitmap_cubical_complex Bitmap_cubical_complex_periodic_boundary_conditions; + + BOOST_AUTO_TEST_CASE(check_dimension) { std::vector< double > increasingFiltrationOfTopDimensionalCells({1, 2, 3, 4, 5, 6, 7, 8, 9}); @@ -100,24 +102,24 @@ BOOST_AUTO_TEST_CASE(compute_boundary_test_1) { std::vector boundary8; boundary8.push_back(1); boundary8.push_back(15); - boundary8.push_back(7); boundary8.push_back(9); + boundary8.push_back(7); std::vector boundary9; boundary9.push_back(2); boundary9.push_back(16); std::vector boundary10; boundary10.push_back(3); boundary10.push_back(17); - boundary10.push_back(9); boundary10.push_back(11); + boundary10.push_back(9); std::vector boundary11; boundary11.push_back(4); boundary11.push_back(18); std::vector boundary12; boundary12.push_back(5); boundary12.push_back(19); - boundary12.push_back(11); boundary12.push_back(13); + boundary12.push_back(11); std::vector boundary13; boundary13.push_back(6); boundary13.push_back(20); @@ -140,24 +142,24 @@ BOOST_AUTO_TEST_CASE(compute_boundary_test_1) { std::vector boundary22; boundary22.push_back(15); boundary22.push_back(29); - boundary22.push_back(21); boundary22.push_back(23); + boundary22.push_back(21); std::vector boundary23; boundary23.push_back(16); boundary23.push_back(30); std::vector boundary24; boundary24.push_back(17); boundary24.push_back(31); - boundary24.push_back(23); boundary24.push_back(25); + boundary24.push_back(23); std::vector boundary25; boundary25.push_back(18); boundary25.push_back(32); std::vector boundary26; boundary26.push_back(19); boundary26.push_back(33); - boundary26.push_back(25); boundary26.push_back(27); + boundary26.push_back(25); std::vector boundary27; boundary27.push_back(20); boundary27.push_back(34); @@ -180,24 +182,24 @@ BOOST_AUTO_TEST_CASE(compute_boundary_test_1) { std::vector boundary36; boundary36.push_back(29); boundary36.push_back(43); - boundary36.push_back(35); boundary36.push_back(37); + boundary36.push_back(35); std::vector boundary37; boundary37.push_back(30); boundary37.push_back(44); std::vector boundary38; boundary38.push_back(31); boundary38.push_back(45); - boundary38.push_back(37); boundary38.push_back(39); + boundary38.push_back(37); std::vector boundary39; boundary39.push_back(32); boundary39.push_back(46); std::vector boundary40; boundary40.push_back(33); boundary40.push_back(47); - boundary40.push_back(39); boundary40.push_back(41); + boundary40.push_back(39); std::vector boundary41; boundary41.push_back(34); boundary41.push_back(48); @@ -275,11 +277,13 @@ BOOST_AUTO_TEST_CASE(compute_boundary_test_1) { for (size_t i = 0; i != increasing.size(); ++i) { std::vector< size_t > bd = increasing.get_boundary_of_a_cell(i); for (size_t j = 0; j != bd.size(); ++j) { - BOOST_CHECK(boundaries[i][j] == bd[j]); - } + BOOST_CHECK(boundaries[i][j] == bd[j]); + } } } + + BOOST_AUTO_TEST_CASE(compute_boundary_test_2) { std::vector< double > increasingFiltrationOfTopDimensionalCells({1, 2, 3, 4, 5, 6, 7, 8, 9}); @@ -373,16 +377,19 @@ BOOST_AUTO_TEST_CASE(compute_boundary_test_2) { coboundaryElements.push_back(40); coboundaryElements.push_back(41); coboundaryElements.push_back(47); + size_t number = 0; for (size_t i = 0; i != increasing.size(); ++i) { std::vector< size_t > bd = increasing.get_coboundary_of_a_cell(i); for (size_t j = 0; j != bd.size(); ++j) { BOOST_CHECK(coboundaryElements[number] == bd[j]); - ++number; - } + ++number; + } } } + + BOOST_AUTO_TEST_CASE(compute_boundary_test_3) { std::vector< double > increasingFiltrationOfTopDimensionalCells({1, 2, 3, 4, 5, 6, 7, 8, 9}); @@ -578,50 +585,50 @@ BOOST_AUTO_TEST_CASE(boudary_operator_2d_bitmap_with_periodic_bcond) { std::vector boundary0; std::vector boundary1; - boundary1.push_back(0); boundary1.push_back(2); + boundary1.push_back(0); std::vector boundary2; std::vector boundary3; - boundary3.push_back(2); boundary3.push_back(0); + boundary3.push_back(2); std::vector boundary4; - boundary4.push_back(0); boundary4.push_back(8); + boundary4.push_back(0); std::vector boundary5; - boundary5.push_back(1); boundary5.push_back(9); + boundary5.push_back(1); boundary5.push_back(4); boundary5.push_back(6); std::vector boundary6; - boundary6.push_back(2); boundary6.push_back(10); + boundary6.push_back(2); std::vector boundary7; - boundary7.push_back(3); boundary7.push_back(11); + boundary7.push_back(3); boundary7.push_back(6); boundary7.push_back(4); std::vector boundary8; std::vector boundary9; - boundary9.push_back(8); boundary9.push_back(10); + boundary9.push_back(8); std::vector boundary10; std::vector boundary11; - boundary11.push_back(10); boundary11.push_back(8); + boundary11.push_back(10); std::vector boundary12; - boundary12.push_back(8); boundary12.push_back(0); + boundary12.push_back(8); std::vector boundary13; - boundary13.push_back(9); boundary13.push_back(1); + boundary13.push_back(9); boundary13.push_back(12); boundary13.push_back(14); std::vector boundary14; - boundary14.push_back(10); boundary14.push_back(2); + boundary14.push_back(10); std::vector boundary15; - boundary15.push_back(11); boundary15.push_back(3); + boundary15.push_back(11); boundary15.push_back(14); boundary15.push_back(12); @@ -647,7 +654,7 @@ BOOST_AUTO_TEST_CASE(boudary_operator_2d_bitmap_with_periodic_bcond) { std::vector< size_t > bd = cmplx.get_boundary_of_a_cell(i); for (size_t j = 0; j != bd.size(); ++j) { BOOST_CHECK(boundaries[i][j] == bd[j]); - } + } } } @@ -771,6 +778,8 @@ BOOST_AUTO_TEST_CASE(bitmap_2d_with_periodic_bcond_filtration) { BOOST_CHECK(filtration[i] == cmplx.get_cell_data(i)); } } + + BOOST_AUTO_TEST_CASE(all_cells_iterator_and_boundary_iterators_in_Bitmap_cubical_complex_base_check) { std::vector< double > expected_filtration; @@ -836,14 +845,14 @@ BOOST_AUTO_TEST_CASE(all_cells_iterator_and_boundary_iterators_in_Bitmap_cubical expected_boundary.push_back(10); expected_boundary.push_back(1); expected_boundary.push_back(11); - expected_boundary.push_back(5); expected_boundary.push_back(7); + expected_boundary.push_back(5); expected_boundary.push_back(2); expected_boundary.push_back(12); expected_boundary.push_back(3); expected_boundary.push_back(13); - expected_boundary.push_back(7); expected_boundary.push_back(9); + expected_boundary.push_back(7); expected_boundary.push_back(4); expected_boundary.push_back(14); expected_boundary.push_back(10); @@ -854,14 +863,14 @@ BOOST_AUTO_TEST_CASE(all_cells_iterator_and_boundary_iterators_in_Bitmap_cubical expected_boundary.push_back(20); expected_boundary.push_back(11); expected_boundary.push_back(21); - expected_boundary.push_back(15); expected_boundary.push_back(17); + expected_boundary.push_back(15); expected_boundary.push_back(12); expected_boundary.push_back(22); expected_boundary.push_back(13); expected_boundary.push_back(23); - expected_boundary.push_back(17); expected_boundary.push_back(19); + expected_boundary.push_back(17); expected_boundary.push_back(14); expected_boundary.push_back(24); expected_boundary.push_back(20); @@ -936,7 +945,7 @@ BOOST_AUTO_TEST_CASE(all_cells_iterator_and_boundary_iterators_in_Bitmap_cubical Bitmap_cubical_complex_base::Boundary_range bdrange = ba.boundary_range(*it); for ( Bitmap_cubical_complex_base::Boundary_iterator bd = bdrange.begin() ; bd != bdrange.end() ; ++bd ) { - BOOST_CHECK( expected_boundary[bd_it] == *bd ); + BOOST_CHECK( expected_boundary[bd_it] == *bd ); ++bd_it; } @@ -1022,14 +1031,14 @@ BOOST_AUTO_TEST_CASE(all_cells_iterator_and_boundary_iterators_in_Bitmap_cubical expected_boundary.push_back(10); expected_boundary.push_back(1); expected_boundary.push_back(11); - expected_boundary.push_back(5); expected_boundary.push_back(7); + expected_boundary.push_back(5); expected_boundary.push_back(2); expected_boundary.push_back(12); expected_boundary.push_back(3); expected_boundary.push_back(13); - expected_boundary.push_back(7); expected_boundary.push_back(9); + expected_boundary.push_back(7); expected_boundary.push_back(4); expected_boundary.push_back(14); expected_boundary.push_back(10); @@ -1040,14 +1049,14 @@ BOOST_AUTO_TEST_CASE(all_cells_iterator_and_boundary_iterators_in_Bitmap_cubical expected_boundary.push_back(20); expected_boundary.push_back(11); expected_boundary.push_back(21); - expected_boundary.push_back(15); expected_boundary.push_back(17); + expected_boundary.push_back(15); expected_boundary.push_back(12); expected_boundary.push_back(22); expected_boundary.push_back(13); expected_boundary.push_back(23); - expected_boundary.push_back(17); expected_boundary.push_back(19); + expected_boundary.push_back(17); expected_boundary.push_back(14); expected_boundary.push_back(24); expected_boundary.push_back(20); @@ -1124,7 +1133,7 @@ BOOST_AUTO_TEST_CASE(all_cells_iterator_and_boundary_iterators_in_Bitmap_cubical Bitmap_cubical_complex_base::Boundary_range bdrange = ba.boundary_range(*it); for ( Bitmap_cubical_complex_base::Boundary_iterator bd = bdrange.begin() ; bd != bdrange.end() ; ++bd ) { - BOOST_CHECK( expected_boundary[bd_it] == *bd ); + BOOST_CHECK( expected_boundary[bd_it] == *bd ); ++bd_it; } @@ -1208,14 +1217,14 @@ BOOST_AUTO_TEST_CASE(all_cells_iterator_and_boundary_iterators_in_Bitmap_cubical expected_boundary.push_back(10); expected_boundary.push_back(1); expected_boundary.push_back(11); - expected_boundary.push_back(5); expected_boundary.push_back(7); + expected_boundary.push_back(5); expected_boundary.push_back(2); expected_boundary.push_back(12); expected_boundary.push_back(3); expected_boundary.push_back(13); - expected_boundary.push_back(7); expected_boundary.push_back(9); + expected_boundary.push_back(7); expected_boundary.push_back(4); expected_boundary.push_back(14); expected_boundary.push_back(10); @@ -1226,14 +1235,14 @@ BOOST_AUTO_TEST_CASE(all_cells_iterator_and_boundary_iterators_in_Bitmap_cubical expected_boundary.push_back(20); expected_boundary.push_back(11); expected_boundary.push_back(21); - expected_boundary.push_back(15); expected_boundary.push_back(17); + expected_boundary.push_back(15); expected_boundary.push_back(12); expected_boundary.push_back(22); expected_boundary.push_back(13); expected_boundary.push_back(23); - expected_boundary.push_back(17); expected_boundary.push_back(19); + expected_boundary.push_back(17); expected_boundary.push_back(14); expected_boundary.push_back(24); expected_boundary.push_back(20); @@ -1310,7 +1319,7 @@ BOOST_AUTO_TEST_CASE(all_cells_iterator_and_boundary_iterators_in_Bitmap_cubical Bitmap_cubical_complex_base::Boundary_range bdrange = ba.boundary_range(*it); for ( Bitmap_cubical_complex_base::Boundary_iterator bd = bdrange.begin() ; bd != bdrange.end() ; ++bd ) { - BOOST_CHECK( expected_boundary[bd_it] == *bd ); + BOOST_CHECK( expected_boundary[bd_it] == *bd ); ++bd_it; } @@ -1376,3 +1385,290 @@ BOOST_AUTO_TEST_CASE(Top_dimensional_cells_iterator_range_check) } } + + + +BOOST_AUTO_TEST_CASE( check_if_boundary_of_boundary_is_zero_non_periodic_case_3_d ) +{ + std::vector< unsigned > sizes(3); + sizes[0] = 5; + sizes[1] = 5; + sizes[2] = 5; + + std::vector< double > data(125 , 0); + + int number_of_all_elements = (2*sizes[0]+1)*(2*sizes[1]+1)*(2*sizes[2]+1); + Bitmap_cubical_complex_base ba( sizes , data ); + for ( Bitmap_cubical_complex_base::All_cells_iterator it = ba.all_cells_iterator_begin() ; it != ba.all_cells_iterator_end() ; ++it ) + { + std::vector< int > elems_in_boundary( number_of_all_elements,0 ); + int i = 1; + Bitmap_cubical_complex_base::Boundary_range bdrange = ba.boundary_range(*it); + for ( Bitmap_cubical_complex_base::Boundary_iterator bd = bdrange.begin() ; bd != bdrange.end() ; ++bd ) + { + Bitmap_cubical_complex_base::Boundary_range second_bdrange = ba.boundary_range(*bd); + int j = 1; + for ( Bitmap_cubical_complex_base::Boundary_iterator bd2 = second_bdrange.begin() ; bd2 != second_bdrange.end() ; ++bd2 ) + { + elems_in_boundary[ *bd2 ] += i*j; + j *= -1; + } + i *= -1; + } + //check if there is anything nonzero in elems_in_boundary + for ( size_t i = 0 ; i != elems_in_boundary.size() ; ++i ) + { + BOOST_CHECK( elems_in_boundary[i] == 0 ); + } + } +} + + + + +BOOST_AUTO_TEST_CASE( check_if_boundary_of_boundary_is_zero_non_periodic_case_4_d ) +{ + std::vector< unsigned > sizes(4); + sizes[0] = 3; + sizes[1] = 3; + sizes[2] = 3; + sizes[3] = 3; + + std::vector< double > data(81 , 0); + + int number_of_all_elements = (2*sizes[0]+1)*(2*sizes[1]+1)*(2*sizes[2]+1)*(2*sizes[3]+1); + Bitmap_cubical_complex_base ba( sizes , data ); + for ( Bitmap_cubical_complex_base::All_cells_iterator it = ba.all_cells_iterator_begin() ; it != ba.all_cells_iterator_end() ; ++it ) + { + std::vector< int > elems_in_boundary( number_of_all_elements,0 ); + int i = 1; + Bitmap_cubical_complex_base::Boundary_range bdrange = ba.boundary_range(*it); + for ( Bitmap_cubical_complex_base::Boundary_iterator bd = bdrange.begin() ; bd != bdrange.end() ; ++bd ) + { + Bitmap_cubical_complex_base::Boundary_range second_bdrange = ba.boundary_range(*bd); + int j = 1; + for ( Bitmap_cubical_complex_base::Boundary_iterator bd2 = second_bdrange.begin() ; bd2 != second_bdrange.end() ; ++bd2 ) + { + elems_in_boundary[ *bd2 ] += i*j; + j *= -1; + } + i *= -1; + } + //check if there is anything nonzero in elems_in_boundary + for ( size_t i = 0 ; i != elems_in_boundary.size() ; ++i ) + { + BOOST_CHECK( elems_in_boundary[i] == 0 ); + } + } +} + + + + + + +BOOST_AUTO_TEST_CASE( check_if_boundary_of_boundary_is_zero_periodic_case_2d ) +{ + std::vector< unsigned > sizes(2); + sizes[0] = 10; + sizes[1] = 10; + + std::vector directions_of_periodicity(2,true); + std::vector< double > data(100 , 0); + + int number_of_all_elements = (2*sizes[0])*(2*sizes[1]);// *(2*sizes[2]); + Bitmap_cubical_complex_periodic_boundary_conditions ba( sizes , data , directions_of_periodicity ); + for ( Bitmap_cubical_complex_periodic_boundary_conditions::All_cells_iterator it = ba.all_cells_iterator_begin() ; it != ba.all_cells_iterator_end() ; ++it ) + { + std::vector< int > elems_in_boundary( number_of_all_elements,0 ); + int i = 1; + + //std::cout << "Element : " << *it << std::endl; + + Bitmap_cubical_complex_periodic_boundary_conditions_base::Boundary_range bdrange = ba.boundary_range(*it); + for ( Bitmap_cubical_complex_periodic_boundary_conditions::Boundary_iterator bd = bdrange.begin() ; bd != bdrange.end() ; ++bd ) + { + //std::cout << *bd << " "; + Bitmap_cubical_complex_periodic_boundary_conditions::Boundary_range second_bdrange = ba.boundary_range(*bd); + int j = 1; + for ( Bitmap_cubical_complex_periodic_boundary_conditions::Boundary_iterator bd2 = second_bdrange.begin() ; bd2 != second_bdrange.end() ; ++bd2 ) + { + elems_in_boundary[ *bd2 ] += i*j; + j *= -1; + } + i *= -1; + } + //getchar(); + + //check if there is anything nonzero in elems_in_boundary + for ( size_t i = 0 ; i != elems_in_boundary.size() ; ++i ) + { + BOOST_CHECK( elems_in_boundary[i] == 0 ); + } + } +} + + + +BOOST_AUTO_TEST_CASE( check_if_boundary_of_boundary_is_zero_periodic_case_3d ) +{ + std::vector< unsigned > sizes(3); + sizes[0] = 5; + sizes[1] = 5; + sizes[2] = 5; + + std::vector directions_of_periodicity(3,true); + + std::vector< double > data(125 , 0); + + int number_of_all_elements = (2*sizes[0])*(2*sizes[1])*(2*sizes[2]); + Bitmap_cubical_complex_periodic_boundary_conditions ba( sizes , data , directions_of_periodicity ); + for ( Bitmap_cubical_complex_periodic_boundary_conditions::All_cells_iterator it = ba.all_cells_iterator_begin() ; it != ba.all_cells_iterator_end() ; ++it ) + { + //std::cout << "Element : " << *it << std::endl; + + std::vector< int > elems_in_boundary( number_of_all_elements,0 ); + int i = 1; + + Bitmap_cubical_complex_periodic_boundary_conditions_base::Boundary_range bdrange = ba.boundary_range(*it); + for ( Bitmap_cubical_complex_periodic_boundary_conditions::Boundary_iterator bd = bdrange.begin() ; bd != bdrange.end() ; ++bd ) + { + Bitmap_cubical_complex_periodic_boundary_conditions::Boundary_range second_bdrange = ba.boundary_range(*bd); + //std::cout << *bd << " "; + int j = 1; + for ( Bitmap_cubical_complex_periodic_boundary_conditions::Boundary_iterator bd2 = second_bdrange.begin() ; bd2 != second_bdrange.end() ; ++bd2 ) + { + elems_in_boundary[ *bd2 ] += i*j; + j *= -1; + } + i *= -1; + } + //getchar(); + + //check if there is anything nonzero in elems_in_boundary + for ( size_t i = 0 ; i != elems_in_boundary.size() ; ++i ) + { + BOOST_CHECK( elems_in_boundary[i] == 0 ); + } + } +} + + + + + + +BOOST_AUTO_TEST_CASE( check_if_boundary_of_boundary_is_zero_periodic_case_4d ) +{ + std::vector< unsigned > sizes(4); + sizes[0] = 4; + sizes[1] = 4; + sizes[2] = 4; + sizes[3] = 4; + + std::vector directions_of_periodicity(4,true); + + std::vector< double > data(256 , 0); + + int number_of_all_elements = (2*sizes[0])*(2*sizes[1])*(2*sizes[2])*(2*sizes[3]); + Bitmap_cubical_complex_periodic_boundary_conditions ba( sizes , data , directions_of_periodicity ); + for ( Bitmap_cubical_complex_periodic_boundary_conditions::All_cells_iterator it = ba.all_cells_iterator_begin() ; it != ba.all_cells_iterator_end() ; ++it ) + { + std::vector< int > elems_in_boundary( number_of_all_elements,0 ); + int i = 1; + + Bitmap_cubical_complex_periodic_boundary_conditions_base::Boundary_range bdrange = ba.boundary_range(*it); + for ( Bitmap_cubical_complex_periodic_boundary_conditions::Boundary_iterator bd = bdrange.begin() ; bd != bdrange.end() ; ++bd ) + { + Bitmap_cubical_complex_periodic_boundary_conditions::Boundary_range second_bdrange = ba.boundary_range(*bd); + int j = 1; + for ( Bitmap_cubical_complex_periodic_boundary_conditions::Boundary_iterator bd2 = second_bdrange.begin() ; bd2 != second_bdrange.end() ; ++bd2 ) + { + elems_in_boundary[ *bd2 ] += i*j; + j *= -1; + } + i *= -1; + } + + //check if there is anything nonzero in elems_in_boundary + for ( size_t i = 0 ; i != elems_in_boundary.size() ; ++i ) + { + BOOST_CHECK( elems_in_boundary[i] == 0 ); + } + } +} + + + + + +BOOST_AUTO_TEST_CASE( compute_incidence_between_cells_test ) +{ + std::vector< unsigned > sizes(3); + sizes[0] = 5; + sizes[1] = 5; + sizes[2] = 5; + + std::vector< double > data(125 , 0); + + int number_of_all_elements = (2*sizes[0]+1)*(2*sizes[1]+1)*(2*sizes[1]+1); + Bitmap_cubical_complex_base ba( sizes , data ); + for ( Bitmap_cubical_complex_base::All_cells_iterator it = ba.all_cells_iterator_begin() ; it != ba.all_cells_iterator_end() ; ++it ) + { + std::vector< int > elems_in_boundary( number_of_all_elements,0 ); + Bitmap_cubical_complex_base::Boundary_range bdrange = ba.boundary_range(*it); + for ( Bitmap_cubical_complex_base::Boundary_iterator bd = bdrange.begin() ; bd != bdrange.end() ; ++bd ) + { + Bitmap_cubical_complex_base::Boundary_range second_bdrange = ba.boundary_range(*bd); + for ( Bitmap_cubical_complex_base::Boundary_iterator bd2 = second_bdrange.begin() ; bd2 != second_bdrange.end() ; ++bd2 ) + { + elems_in_boundary[ *bd2 ] += + ba.compute_incidence_between_cells( *it , *bd ) * + ba.compute_incidence_between_cells( *bd , *bd2 ); + } + } + //check if there is anything nonzero in elems_in_boundary + for ( size_t i = 0 ; i != elems_in_boundary.size() ; ++i ) + { + BOOST_CHECK( elems_in_boundary[i] == 0 ); + } + } +} + + + +BOOST_AUTO_TEST_CASE( compute_incidence_between_cells_test_periodic_boundary_conditions ) +{ + std::vector< unsigned > sizes(3); + sizes[0] = 5; + sizes[1] = 5; + sizes[2] = 5; + + std::vector directions_of_periodicity(2,true); + std::vector< double > data(125 , 0); + + int number_of_all_elements = (2*sizes[0])*(2*sizes[1])*(2*sizes[2]); + Bitmap_cubical_complex_periodic_boundary_conditions ba( sizes , data , directions_of_periodicity ); + for ( Bitmap_cubical_complex_periodic_boundary_conditions::All_cells_iterator it = ba.all_cells_iterator_begin() ; it != ba.all_cells_iterator_end() ; ++it ) + { + std::vector< int > elems_in_boundary( number_of_all_elements,0 ); + Bitmap_cubical_complex_periodic_boundary_conditions_base::Boundary_range bdrange = ba.boundary_range(*it); + for ( Bitmap_cubical_complex_periodic_boundary_conditions::Boundary_iterator bd = bdrange.begin() ; bd != bdrange.end() ; ++bd ) + { + //std::cout << *bd << " "; + Bitmap_cubical_complex_periodic_boundary_conditions::Boundary_range second_bdrange = ba.boundary_range(*bd); + for ( Bitmap_cubical_complex_periodic_boundary_conditions::Boundary_iterator bd2 = second_bdrange.begin() ; bd2 != second_bdrange.end() ; ++bd2 ) + { + elems_in_boundary[ *bd2 ] += + ba.compute_incidence_between_cells( *it , *bd ) * + ba.compute_incidence_between_cells( *bd , *bd2 ); + } + } + //check if there is anything nonzero in elems_in_boundary + for ( size_t i = 0 ; i != elems_in_boundary.size() ; ++i ) + { + BOOST_CHECK( elems_in_boundary[i] == 0 ); + } + } +} + -- cgit v1.2.3