From 67c20567ea92564b8669a537dc3d4fe70e5ff77c Mon Sep 17 00:00:00 2001
From: Mathieu Dutour Sikiric <mathieu.dutour@gmail.com>
Date: Sat, 4 Jan 2025 13:11:46 +0100
Subject: [PATCH] Insert one test case and some other updates.

---
 CI_tests/Tspaces_IsoDelaunay/TSPACE_LinSpa | 39 ++++++++++++++++++++++
 src_group/GRP_GroupFct.h                   | 17 ++++++++++
 src_latt/IsoDelaunayDomains.h              |  3 ++
 src_latt/Makefile_tspace_isodelaunay       |  2 +-
 4 files changed, 60 insertions(+), 1 deletion(-)
 create mode 100644 CI_tests/Tspaces_IsoDelaunay/TSPACE_LinSpa

diff --git a/CI_tests/Tspaces_IsoDelaunay/TSPACE_LinSpa b/CI_tests/Tspaces_IsoDelaunay/TSPACE_LinSpa
new file mode 100644
index 00000000..e03d381a
--- /dev/null
+++ 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
index eb3aea26..797bc1f5 100644
--- 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
index 7378ba4e..ccbf9676 100644
--- 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
index 50ba77fb..39d95d58 100644
--- 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)