Insert one test case and some other updates.

MathieuDutSik · Jan 4, 2025 · 67c2056 · 67c2056
1 parent 2fd5d45
commit 67c2056
Show file tree

Hide file tree

Showing 4 changed files with 60 additions and 1 deletion.
diff --git a/CI_tests/Tspaces_IsoDelaunay/TSPACE_LinSpa b/CI_tests/Tspaces_IsoDelaunay/TSPACE_LinSpa
@@ -0,0 +1,39 @@
+4 4
+ 2 1 0 0
+ 1 2 1 1
+ 0 1 2 0
+ 0 1 0 2
+3
+4 4
+ 0 0 0 1
+ 0 0 1 0
+ 0 1 2 0
+ 1 0 0 0
+4 4
+ -2 -1 0 0
+ -1 0 1 0
+ 0 1 0 1
+ 0 0 1 0
+4 4
+ 2 1 0 0
+ 1 2 1 1
+ 0 1 2 0
+ 0 1 0 2
+0
+0
+3
+4 4
+ -1 0 0 0
+ 0 -1 0 0
+ 0 0 -1 0
+ 0 0 0 -1
+4 4
+ 1 0 0 0
+ 1 -1 0 0
+ 0 0 -1 0
+ 1 -2 1 1
+4 4
+ -1 1 -1 0
+ 0 -1 0 1
+ 1 -1 0 1
+ 0 -1 0 0
diff --git a/src_group/GRP_GroupFct.h b/src_group/GRP_GroupFct.h
@@ -151,6 +151,23 @@ void WriteGroup(std::ostream &os, Tgroup const &TheGRP) {
   }
 }
 
+template<typename Tgroup>
+void PrintRepresentativeAction_OnSets_GRP_f1_f2(Tgroup const& GRP, Face const& f1, Face const& f2) {
+  std::string prefix = "RepresentativeAction_OnSets_GRP_f1_f2_idx";
+  std::string FileOut = FindAvailableFileFromPrefix(prefix);
+  std::ofstream os(FileOut);
+  WriteGroup(os, GRP);
+  auto f_print=[&](Face const& f) -> void {
+    for (size_t i=0; i<f.size(); i++) {
+      os << " " << f[i];
+    }
+    os << "\n";
+  };
+  f_print(f1);
+  f_print(f2);
+}
+
+
 template <typename Tgroup> std::string StringGroup(Tgroup const &TheGRP) {
   using Telt = typename Tgroup::Telt;
   using Tidx = typename Telt::Tidx;

diff --git a/src_latt/IsoDelaunayDomains.h b/src_latt/IsoDelaunayDomains.h
@@ -1336,6 +1336,9 @@ FlippingLtype(DelaunayTesselation<Tvert, Tgroup> const &ListOrbitDelaunay,
     CheckFacetInequality(EXT, eInc, "get_matching_old_tessel EXT eInc");
 #endif
     for (auto &eAdj : ListOrbitDelaunay.l_dels[iDelaunayOld].ListAdj) {
+#ifdef PRINT_ISO_DELAUNAY_DOMAIN_REPRESENTATIVE_ACTION
+      PrintRepresentativeAction_OnSets_GRP_f1_f2(ListOrbitDelaunay.l_dels[iDelaunayOld].GRP, eAdj.eInc, eInc);
+#endif
       std::optional<Telt> opt =
           ListOrbitDelaunay.l_dels[iDelaunayOld]
               .GRP.RepresentativeAction_OnSets(eAdj.eInc, eInc);

diff --git a/src_latt/Makefile_tspace_isodelaunay b/src_latt/Makefile_tspace_isodelaunay
@@ -11,7 +11,7 @@ OUTPUTFLAG 	= -o
 CHASHMAP        = -I../basic_common_cpp/sparse-map/include/tsl/ -I../basic_common_cpp/robin-map/include/tsl/ -I../basic_common_cpp/hopscotch-map/include/tsl/
 CLOCAL          = -I../src_basic -I../src_number -I../src_matrix -I../src_lorentzian -I../src_graph -I../permutalib/src -I../src_comb -I../src_poly -I../src_group -I../src_dualdesc -I../src_isotropy -I../src_short -I.
 #COPTIONS        = -DDEBUG -DSANITY_CHECK
-COPTIONS        =
+COPTIONS        = -DPRINT_ISO_DELAUNAY_DOMAIN_REPRESENTATIVE_ACTION
 
 #CFLAGS 		= -I$(GMP_INCDIR) -I$(BOOST_INCDIR) -I$(EIGEN_PATH)
 CFLAGS 		= -I$(GMP_INCDIR) -I$(BOOST_INCDIR) -I$(EIGEN_PATH) $(NAUTY_INCLUDE) $(MPI_INCLUDE_CPP)