1 files changed, 111 insertions, 9 deletions

diff --git a/src/solvers/h48/gendata_h48.h b/src/solvers/h48/gendata_h48.h
index 3445570..75b7abf 100644
--- a/src/solvers/h48/gendata_h48.h
+++ b/src/solvers/h48/gendata_h48.h
@@ -103,14 +103,19 @@ STATIC_INLINE bool gendata_h48k2_dfs_stop(cube_t, int8_t, h48k2_dfs_arg_t *);
 STATIC size_t gendata_h48k2_realcoord(gendata_h48_arg_t *);
 STATIC void gendata_h48k2_dfs(h48k2_dfs_arg_t *arg);
 STATIC void * gendata_h48k2_runthread(void *);
-STATIC tableinfo_t makeinfo_h48k2(gendata_h48_arg_t *, uint8_t);
+STATIC tableinfo_t makeinfo_h48k2(gendata_h48_arg_t *);
+STATIC void getdistribution_h48(const uint8_t *,
+    uint64_t [static INFO_DISTRIBUTION_LEN], uint8_t, uint8_t);
 
 STATIC uint32_t *get_cocsepdata_ptr(const void *);
 STATIC uint8_t *get_h48data_ptr(const void *);
 
 STATIC_INLINE uint8_t get_h48_pval(const uint8_t *, int64_t, uint8_t);
 STATIC_INLINE void set_h48_pval(uint8_t *, int64_t, uint8_t, uint8_t);
-STATIC_INLINE uint8_t get_h48_bound(cube_t, uint32_t, uint8_t, uint8_t, uint8_t *);
+STATIC_INLINE uint8_t get_h48_bound(
+    cube_t, uint32_t, uint8_t, uint8_t, uint8_t *);
+
+size_t gendata_h48_derive(uint8_t, const void *, void *);
 
 STATIC uint64_t
 gendata_h48short(gendata_h48short_arg_t *arg)
@@ -387,7 +392,7 @@ gendata_h48k2(gendata_h48_arg_t *arg)
                 [11] = 10
         };
 
-        uint8_t t, selectedbase, *table;
+        uint8_t t, *table;
         int64_t j;
         uint64_t i, ii, inext, count;
         h48map_t shortcubes;
@@ -414,8 +419,9 @@ gendata_h48k2(gendata_h48_arg_t *arg)
         };
         gendata_h48short(&shortarg);
 
-        selectedbase = arg->base < 20 ? arg->base : base[arg->h];
-        arg->info = makeinfo_h48k2(arg, selectedbase);
+        if (arg->base >= 20)
+                arg->base = base[arg->h];
+        arg->info = makeinfo_h48k2(arg);
 
         inext = count = 0;
         pthread_mutex_init(&shortcubes_mutex, NULL);
@@ -425,7 +431,7 @@ gendata_h48k2(gendata_h48_arg_t *arg)
                 dfsarg[i] = (h48k2_dfs_arg_t){
                         .h = arg->h,
                         .k = arg->k,
-                        .base = selectedbase,
+                        .base = arg->base,
                         .shortdepth = shortdepth,
                         .cocsepdata = arg->cocsepdata,
                         .table = table,
@@ -603,14 +609,17 @@ STATIC_INLINE bool
 gendata_h48k2_dfs_stop(cube_t cube, int8_t depth, h48k2_dfs_arg_t *arg)
 {
         uint64_t val;
-        int64_t coord;
+        int64_t coord, mutex;
         int8_t oldval;
 
         if (arg->h == 0 || arg->h == 11) {
                 /* We are in the "real coordinate" case, we can stop
                    if this coordinate has already been visited */
                 coord = coord_h48(cube, arg->cocsepdata, arg->h);
+                mutex = H48_INDEX(coord, arg->k) % CHUNKS;
+                pthread_mutex_lock(arg->table_mutex[mutex]);
                 oldval = get_h48_pval(arg->table, coord, arg->k);
+                pthread_mutex_unlock(arg->table_mutex[mutex]);
                 return oldval <= depth;
         } else {
                 /* With 0 < k < 11 we do not have a "real coordinate".
@@ -637,7 +646,7 @@ gendata_h48k2_realcoord_runthread(void *arg)
 }
 
 STATIC tableinfo_t
-makeinfo_h48k2(gendata_h48_arg_t *arg, uint8_t base)
+makeinfo_h48k2(gendata_h48_arg_t *arg)
 {
         tableinfo_t info;
 
@@ -651,7 +660,7 @@ makeinfo_h48k2(gendata_h48_arg_t *arg, uint8_t base)
                 .classes = 0,
                 .h48h = arg->h,
                 .bits = 2,
-                .base = base,
+                .base = arg->base,
                 .maxvalue = 3,
                 .next = 0,
         };
@@ -662,6 +671,25 @@ makeinfo_h48k2(gendata_h48_arg_t *arg, uint8_t base)
         return info;
 }
 
+STATIC void
+getdistribution_h48(
+        const uint8_t *table,
+        uint64_t distr[static INFO_DISTRIBUTION_LEN],
+        uint8_t h,
+        uint8_t k
+) {
+        uint8_t val;
+        int64_t i, h48max;
+
+        memset(distr, 0, INFO_DISTRIBUTION_LEN * sizeof(uint64_t));
+
+        h48max = H48_COORDMAX(h);
+        for (i = 0; i < h48max; i++) {
+                val = get_h48_pval(table, i, k);
+                distr[val]++;
+        }
+}
+
 STATIC uint32_t *
 get_cocsepdata_ptr(const void *data)
 {
@@ -695,3 +723,77 @@ get_h48_bound(cube_t cube, uint32_t cdata, uint8_t h, uint8_t k, uint8_t *table)
         coord = coord_h48_edges(cube, COCLASS(cdata), TTREP(cdata), h);
         return get_h48_pval(table, coord, k);
 }
+
+size_t
+gendata_h48_derive(uint8_t h, const void *fulltable, void *buf)
+{
+        size_t cocsepsize, h48size;
+        uint8_t val_full, val_derive, *h48full, *h48derive;
+        int64_t i, j, h48max;
+        gendata_h48_arg_t arg;
+        tableinfo_t cocsepinfo, fulltableinfo;
+
+        /* Initializing values in case of error */
+        /* TODO cleanup this */
+        fulltableinfo.bits = 2;
+        fulltableinfo.base = 8;
+
+        readtableinfo_n(fulltable, 2, &fulltableinfo);
+        arg.h = h;
+        arg.k = fulltableinfo.bits;
+        arg.maxdepth = 20;
+        arg.buf = buf;
+        arg.cocsepdata = (uint32_t *)((char *)buf + INFOSIZE);
+        arg.base = fulltableinfo.base;
+        arg.info = makeinfo_h48k2(&arg);
+
+        /* Technically this step is redundant, except that we
+           need selfsim and crep */
+        cocsepsize = gendata_cocsep(buf, arg.selfsim, arg.crep);
+        arg.h48buf = (uint8_t *)buf + cocsepsize;
+        h48size = H48_TABLESIZE(h, arg.k) + INFOSIZE;
+
+        if (buf == NULL)
+                goto gendata_h48_derive_return_size;
+
+        if (!readtableinfo(buf, &cocsepinfo)) {
+                LOG("gendata_h48: could not read info for cocsep table\n");
+                goto gendata_h48_derive_error;
+        }
+
+        cocsepinfo.next = cocsepsize;
+        if (!writetableinfo(&cocsepinfo, buf)) {
+                LOG("gendata_h48_derive: could not write info for cocsep table"
+                    " with updated 'next' value\n");
+                goto gendata_h48_derive_error;
+        }
+
+        h48full = (uint8_t *)fulltable + cocsepsize + INFOSIZE;
+        h48derive = (uint8_t *)arg.h48buf + INFOSIZE;
+        memset(h48derive, 0xFF, H48_TABLESIZE(h, arg.k));
+        memset(arg.info.distribution, 0,
+            INFO_DISTRIBUTION_LEN * sizeof(uint64_t));
+
+        h48max = H48_COORDMAX(fulltableinfo.h48h);
+        for (i = 0; i < h48max; i++) {
+                if (i % INT64_C(1000000000) == 0 && i > 0)
+                        LOG("Processing %" PRId64 "th coordinate\n", i);
+                j = i >> (int64_t)(fulltableinfo.h48h - h);
+                val_full = get_h48_pval(h48full, i, arg.k);
+                val_derive = get_h48_pval(h48derive, j, arg.k);
+                set_h48_pval(h48derive, j, arg.k, MIN(val_full, val_derive));
+        }
+
+        getdistribution_h48(h48derive, arg.info.distribution, h, arg.k);
+
+        if (!writetableinfo(&arg.info, arg.h48buf)) {
+                LOG("gendata_h48_derive: could not write info for table\n");
+                goto gendata_h48_derive_error;
+        }
+
+gendata_h48_derive_return_size:
+        return cocsepsize + h48size;
+
+gendata_h48_derive_error:
+        return 0;
+}