From b58138cd941f8d8efd05e04eb9c6b604f24aae1d Mon Sep 17 00:00:00 2001
From: Zeno Rogue <zeno@attnam.com>
Date: Sun, 26 Jan 2020 01:18:30 +0100
Subject: [PATCH] fieldpattern:: load found solutions from shstreams

---
 fieldpattern.cpp | 53 +++++++++++++++++++++++++++++++++++++++++-------
 geom-exp.cpp     |  7 +++++++
 reg3.cpp         | 40 ++++++++++++++++++++----------------
 3 files changed, 76 insertions(+), 24 deletions(-)

diff --git a/fieldpattern.cpp b/fieldpattern.cpp
index 85a389ca..04788123 100644
--- a/fieldpattern.cpp
+++ b/fieldpattern.cpp
@@ -307,7 +307,8 @@ struct fpattern {
   vector<matrix> generate_isometries3();
   int solve3();
   bool generate_all3();
-  
+
+  void set_field(int p, int sq);
   // transmatrix full_R, full_P, full_X;
   
   #if CAP_THREAD
@@ -468,6 +469,9 @@ unsigned fpattern::compute_hash() {
   }
 
 bool fpattern::generate_all3() {
+
+  reg3::generate_fulls();
+
   matrices.clear();
   matcode.clear();
   add1(Id);
@@ -489,7 +493,6 @@ bool fpattern::generate_all3() {
   }
 
 int fpattern::solve3() {
-  reg3::generate_cellrotations();
   reg3::construct_relations();
   
   DEBB(DF_FIELD, ("generating isometries for ", Field));
@@ -551,6 +554,13 @@ int fpattern::solve3() {
   return cmb;
   }
 
+void fpattern::set_field(int p, int sq) {
+  Prime = p;
+  Field = sq ? Prime*Prime : Prime;
+  wsquare = sq;
+  for(int a=0; a<MWDIM; a++) for(int b=0; b<MWDIM; b++) Id[a][b] = a==b?1:0;
+  }
+
 int fpattern::solve() {
   
   for(int a=0; a<MWDIM; a++) for(int b=0; b<MWDIM; b++) Id[a][b] = a==b?1:0;
@@ -1063,6 +1073,9 @@ EX fpattern current_quotient_field = fpattern(0);
 EX fpattern fp_invalid = fpattern(0);
 EX bool quotient_field_changed;
 
+// these strings contain \x00
+#define STR(x) string(x, sizeof(x))
+
 EX struct fpattern& getcurrfp() {
   if(geometry == gFieldQuotient && quotient_field_changed)
     return current_quotient_field;
@@ -1075,17 +1088,21 @@ EX struct fpattern& getcurrfp() {
     // 260 cells, hash = 72414D0C (not 0C62E214)
     static fpattern fp(0);
     if(fp.Prime) return fp;
-    fp.Prime = 5; fp.force_hash = 0x72414D0C; fp.solve();
-    // fp.Prime = 11; fp.force_hash = 0x5FC4CFF0; fp.solve();
+    // fp.Prime = 5; fp.force_hash = 0x72414D0C; fp.solve();
+    // fp.Prime = 11; fp.force_hash = 0x5FC4CFF0; fp.solve();    
+    shstream ins(STR("\x05\x00\x00\x00\x02\x00\x00\x00\x03\x00\x00\x00\xfc\xff\xff\xff\x01\x00\x00\x00\x04\x00\x00\x00\xfc\xff\xff\xff\x04\x00\x00\x00\xfe\xff\xff\xff\x00\x00\x00\x00\x01\x00\x00\x00\xfe\xff\xff\xff\x04\x00\x00\x00\x02\x00\x00\x00\x01\x00\x00\x00\x00\x00\x00\x00\x03\x00\x00\x00\x04\x00\x00\x00\x01\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x02\x00\x00\x00\xff\xff\xff\xff\x00\x00\x00\x00\x00\x00\x00\x00\xfc\xff\xff\xff\x02\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x01\x00\x00\x00\x00\x00\x00\x00\xfc\xff\xff\xff\x02\x00\x00\x00\x00\x00\x00\x00\xfc\xff\xff\xff\x01\x00\x00\x00\xfd\xff\xff\xff\x00\x00\x00\x00\x02\x00\x00\x00\xfd\xff\xff\xff\x03\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x01\x00\x00\x00"));
+    hread_fpattern(ins, fp);
     return fp;
     }
   if(geometry == gSpace435) {
     // 650 cells, hash = EB201050
     static fpattern fp(0);
     if(fp.Prime) return fp;
-    fp.Prime = 5; fp.force_hash = 0x72414D0C; fp.solve();
+    // fp.Prime = 5; fp.force_hash = 0xEB201050; fp.solve();
+    // what is 0x72414D0C??
     // fp.Prime = 11; fp.force_hash = 0x65CE0C00; fp.solve();
-
+    shstream ins(STR("\x05\x00\x00\x00\x02\x00\x00\x00\x00\x00\x00\x00\xfc\xff\xff\xff\xff\xff\xff\xff\xfe\xff\xff\xff\xfc\xff\xff\xff\x04\x00\x00\x00\x02\x00\x00\x00\x04\x00\x00\x00\xff\xff\xff\xff\x02\x00\x00\x00\x03\x00\x00\x00\x03\x00\x00\x00\xfd\xff\xff\xff\x01\x00\x00\x00\x02\x00\x00\x00\x03\x00\x00\x00\x02\x00\x00\x00\x00\x00\x00\x00\xfc\xff\xff\xff\x00\x00\x00\x00\x00\x00\x00\x00\x01\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\xff\xff\xff\xff\x00\x00\x00\x00\x02\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x01\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xfd\xff\xff\xff\x00\x00\x00\x00\xfd\xff\xff\xff\x02\x00\x00\x00\x03\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\x03\x00\x00\x00\x02\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x01\x00\x00\x00"));
+    hread_fpattern(ins, fp);
     return fp;
     }
   if(geometry == gSpace336) {
@@ -1094,7 +1111,9 @@ EX struct fpattern& getcurrfp() {
     // 9408 cells in C4089F34
     static fpattern fp(0);
     if(fp.Prime) return fp;
-    fp.Prime = 7; fp.force_hash = 0xE3F6B7BCu; fp.solve();
+    // fp.Prime = 7; fp.force_hash = 0xE3F6B7BCu; fp.solve();
+    shstream ins(STR("\x07\x00\x00\x00\x03\x00\x00\x00\xfa\xff\xff\xff\x02\x00\x00\x00\x03\x00\x00\x00\x06\x00\x00\x00\x02\x00\x00\x00\xfe\xff\xff\xff\xfb\xff\xff\xff\xfc\xff\xff\xff\x03\x00\x00\x00\xfb\xff\xff\xff\xfd\xff\xff\xff\xfb\xff\xff\xff\x01\x00\x00\x00\xfd\xff\xff\xff\xfe\xff\xff\xff\xfd\xff\xff\xff\x03\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\x00\x00\x00\x00\x00\x00\x00\x00\x01\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\xfc\xff\xff\xff\x00\x00\x00\x00\x03\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x01\x00\x00\x00\x00\x00\x00\x00\xfa\xff\xff\xff\xfb\xff\xff\xff\x00\x00\x00\x00\xfa\xff\xff\xff\x02\x00\x00\x00\x06\x00\x00\x00\x00\x00\x00\x00\xfb\xff\xff\xff\x06\x00\x00\x00\x04\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x01\x00\x00\x00"));
+    hread_fpattern(ins, fp);
     return fp;
     }
   if(geometry == gSpace344) {
@@ -1145,6 +1164,8 @@ EX struct fpattern& getcurrfp() {
   return fp_invalid;
   }
 
+#undef STR
+
 // todo undefined behavior
 EX int subpathid = currfp.matcode[currfp.strtomatrix("RRRPRRRRRPRRRP")];
 EX int subpathorder = currfp.order(currfp.matrices[subpathid]);
@@ -1289,5 +1310,23 @@ EX int currfp_get_R() { return currfp.matcode[currfp.R]; }
 EX int currfp_get_P() { return currfp.matcode[currfp.P]; }
 EX int currfp_get_X() { return currfp.matcode[currfp.X]; }
 
+EX void hread_fpattern(hstream& hs, fieldpattern::fpattern& fp) {
+  hread(hs, fp.Prime);
+  hread(hs, fp.wsquare);
+  hread(hs, fp.P);
+  hread(hs, fp.R);
+  hread(hs, fp.X);
+  fp.set_field(fp.Prime, fp.wsquare);
+  fp.generate_all3();
+  }
+
+EX void hwrite_fpattern(hstream& hs, fieldpattern::fpattern& fp) {
+  hwrite(hs, fp.Prime);
+  hwrite(hs, fp.wsquare);
+  hwrite(hs, fp.P);
+  hwrite(hs, fp.R);
+  hwrite(hs, fp.X);
+  }
+
 }
 #endif
diff --git a/geom-exp.cpp b/geom-exp.cpp
index 7afbe1cb..81ca861b 100644
--- a/geom-exp.cpp
+++ b/geom-exp.cpp
@@ -988,6 +988,13 @@ int read_geom_args() {
     fieldpattern::subpathid = currfp.matcode[M];
     fieldpattern::subpathorder = currfp.order(M);
     }
+  else if(argis("-fwrite")) {
+    shstream hs;
+    hwrite_fpattern(hs, currfp);
+    string s;
+    for(char c: hs.s) s += format("\\x%02x", (unsigned char) c);
+    println(hlog, "current fieldpattern: ", s);
+    }
   else if(argis("-csp")) {
     cheat();
     currfp.findsubpath();
diff --git a/reg3.cpp b/reg3.cpp
index ecf7f24d..3b2f4daf 100644
--- a/reg3.cpp
+++ b/reg3.cpp
@@ -1084,6 +1084,27 @@ EX int matrix_order(const transmatrix A) {
   return res;
   }
 
+EX void generate_fulls() {
+  reg3::generate();
+  reg3::generate_cellrotations();
+
+  auto cons = [&] (int i0, int i1, int i2) {
+    using reg3::adjmoves;
+    transmatrix T = build_matrix(adjmoves[ 0]*C0, adjmoves[ 1]*C0, adjmoves[ 2]*C0, C0);
+    transmatrix U = build_matrix(adjmoves[i0]*C0, adjmoves[i1]*C0, adjmoves[i2]*C0, C0);
+    return U * inverse(T);
+    };
+  
+  full_P = reg3::adjmoves[0];
+  full_R = S7 == 8 ? cons(1, 7, 0) : cons(1, 2, 0);
+  full_X = S7 == 8 ? cons(1, 0, 6) : S7 == 6 ? cons(1, 0, 5) : cons(1, 0, reg3::face);
+  
+  xp_order = matrix_order(full_X * full_P);
+  r_order = matrix_order(full_R);
+  rx_order = matrix_order(full_R * full_X);
+  println(hlog, "orders = ", tie(rx_order, r_order, xp_order));
+  }
+
 EX void construct_relations() {
   if(for_cgi == &cgi) return;
   for_cgi = &cgi;
@@ -1091,6 +1112,7 @@ EX void construct_relations() {
 
   reg3::generate();
   reg3::generate_cellrotations();
+  reg3::generate_fulls();
   vector<transmatrix> all;
 
   vector<string> formulas;
@@ -1106,17 +1128,6 @@ EX void construct_relations() {
     return S;
     };
   
-  auto cons = [&] (int i0, int i1, int i2) {
-    using reg3::adjmoves;
-    transmatrix T = build_matrix(adjmoves[ 0]*C0, adjmoves[ 1]*C0, adjmoves[ 2]*C0, C0);
-    transmatrix U = build_matrix(adjmoves[i0]*C0, adjmoves[i1]*C0, adjmoves[i2]*C0, C0);
-    return U * inverse(T);
-    };
-  
-  full_P = reg3::adjmoves[0];
-  full_R = S7 == 8 ? cons(1, 7, 0) : cons(1, 2, 0);
-  full_X = S7 == 8 ? cons(1, 0, 6) : S7 == 6 ? cons(1, 0, 5) : cons(1, 0, reg3::face);
-  
   println(hlog, reg3::cellshape);
 
   println(hlog, "cellshape = ", isize(reg3::cellshape));
@@ -1153,12 +1164,7 @@ EX void construct_relations() {
     work(T * full_R, i, 'R');
     work(T * full_X, i, 'X');
     work(T * full_P, i, 'P');
-    }
-  
-  xp_order = matrix_order(full_X * full_P);
-  r_order = matrix_order(full_R);
-  rx_order = matrix_order(full_R * full_X);
-  println(hlog, "orders = ", tie(rx_order, r_order, xp_order));
+    }  
   }
 
 EX }