marching-cubes2 (empty) → 0.1.0.0

raw patch · 12 files changed

+2789/−0 lines, 12 filesdep +basedep +containersdep +extrasetup-changed

Dependencies added: base, containers, extra, linear, matrix, split, vector

Files

• C/marchingCubes.c +566/−0
• C/marchingCubes.h +128/−0
• C/tables.c +1334/−0
• C/utils.c +407/−0
• CHANGELOG.md +8/−0
• LICENSE +30/−0
• README.md +12/−0
• Setup.hs +2/−0
• marching-cubes2.cabal +50/−0
• src/MarchingCubes.hs +12/−0
• src/MarchingCubes/MarchingCubes.hs +117/−0
• src/MarchingCubes/Mesh.hs +123/−0

@@ -0,0 +1,566 @@+#ifndef __HEADER__+#define __HEADER__+#include "marchingCubes.h"+#endif++unsigned* FacesNo7(int* faces,+                   size_t* p1,+                   double* values,+                   size_t l,+                   unsigned j) {+  unsigned* index = malloc(l * sizeof(unsigned));+  for(size_t i = 0; i < l; i++) {+    unsigned f = abs(faces[i]) - 1;+    unsigned e1 = FacePoints[f][1];+    unsigned e2 = FacePoints[f][2];+    unsigned e3 = FacePoints[f][3];+    unsigned e4 = FacePoints[f][4];+    size_t p = p1[i] - 2;+    double A = values[p + e1];+    double B = values[p + e2];+    double C = values[p + e3];+    double D = values[p + e4];+    int temp = faces[i] > 0 ? 1 : -1;+    temp *= (A * B - C * D > 0 ? 1 : -1);+    index[i] = temp == 1 ? upow(2, j - 1) : 0;+  }+  return (index);+}++unsigned* Faces7(int* faces, size_t* p1, double* values, size_t l, unsigned j) {+  unsigned* index = malloc(l * sizeof(unsigned));+  for(size_t i = 0; i < l; i++) {+    size_t p = p1[i] - 1;+    double A0 = values[p];+    double B0 = values[p + 3];+    double C0 = values[p + 2];+    double D0 = values[p + 1];+    double A1 = values[p + 4];+    double B1 = values[p + 7];+    double C1 = values[p + 6];+    double D1 = values[p + 5];+    double a = (A1 - A0) * (C1 - C0) - (B1 - B0) * (D1 - D0);+    double b =+        C0 * (A1 - A0) + A0 * (C1 - C0) - D0 * (B1 - B0) - B0 * (D1 - D0);+    double c = A0 * C0 - B0 * D0;+    double tmax = -b / (2 * a);+    double maximum = a * tmax * tmax + b * tmax + c;+    maximum = isnan(maximum) ? -1 : maximum;+    // printf("maximum: %f\n", maximum);+    unsigned cond1 = a < 0 ? 1 : 0;+    unsigned cond2 = tmax > 0 ? 1 : 0;+    unsigned cond3 = tmax < 1 ? 1 : 0;+    unsigned cond4 = maximum > 0 ? 1 : 0;+    unsigned totalcond = cond1 * cond2 * cond3 * cond4;+    int temp = faces[i] > 0 ? 1 : -1;+    temp *= (totalcond == 1 ? 1 : -1);+    index[i] = temp == 1 ? upow(2, j - 1) : 0;+  }+  return index;+}++size_t** faceType(double** M,+                  unsigned m,+                  unsigned n,+                  double level,+                  double max) {+  size_t** L = levelMatrix(M, m, n, level, level < max);+  size_t** minorMat = minorMatrix(L, m, n, m - 1, n - 1);+  size_t** sminorMat2 = scaleMinorMatrix(2, L, m, n, 0, n - 1);+  size_t** sminorMat4 = scaleMinorMatrix(4, L, m, n, 0, 0);+  size_t** sminorMat8 = scaleMinorMatrix(8, L, m, n, m - 1, 0);+  size_t** sum1 = matricialSum(minorMat, sminorMat2, m - 1, n - 1);+  size_t** sum2 = matricialSum(sminorMat4, sminorMat8, m - 1, n - 1);+  size_t** out = matricialSum(sum1, sum2, m - 1, n - 1);+  freeMatrix_s(L, m);+  freeMatrix_s(sminorMat2, m - 1);+  freeMatrix_s(sminorMat4, m - 1);+  freeMatrix_s(sminorMat8, m - 1);+  freeMatrix_s(sum1, m - 1);+  freeMatrix_s(sum2, m - 1);+  return out;+}++unsigned** levCells(double*** A,+                    unsigned nx,+                    unsigned ny,+                    unsigned nz,+                    double level,+                    double max,+                    size_t* outnrow) {+  unsigned** cells = malloc((nz - 1) * sizeof(unsigned*));+  unsigned** types = malloc((nz - 1) * sizeof(unsigned*));+  size_t** bottomTypes = faceType(toMatrix(A, nx, ny, 0), nx, ny, level, max);+  size_t totallength = 0;+  unsigned lengths[nz - 1];+  for(unsigned k = 0; k < nz - 1; k++) {+    size_t** topTypes =+        faceType(toMatrix(A, nx, ny, k + 1), nx, ny, level, max);+    size_t** cellTypes = matricialSum(+        bottomTypes, scaleMatrix(16, topTypes, nx - 1, ny - 1), nx - 1, ny - 1);+    unsigned length;+    unsigned** goodcells01 =+        whichIndicesAndItems(cellTypes, nx - 1, ny - 1, &length);+    cells[k] = malloc(length * sizeof(unsigned));+    types[k] = malloc(length * sizeof(unsigned));+    for(unsigned l = 0; l < length; l++) {+      cells[k][l] = goodcells01[0][l] + (nx - 1) * (ny - 1) * k + 1;+      types[k][l] = goodcells01[1][l];+    }+    bottomTypes = topTypes;+    lengths[k] = length;+    totallength += (size_t)length;+    freeMatrix_s(cellTypes, nx - 1);+    freeMatrix_u(goodcells01, 2);+  }+  unsigned** out = malloc(4 * sizeof(unsigned*));+  out[0] = malloc(totallength * sizeof(unsigned));+  out[1] = malloc(totallength * sizeof(unsigned));+  out[2] = malloc(totallength * sizeof(unsigned));+  out[3] = malloc(totallength * sizeof(unsigned));+  size_t count = 0;+  for(unsigned k = 0; k < nz - 1; k++) {+    for(unsigned l = 0; l < lengths[k]; l++) {+      unsigned c = cells[k][l] - 1;+      out[0][count] = c % (nx - 1) + 1;+      out[1][count] = (c / (nx - 1)) % (ny - 1) + 1;+      out[2][count] = c / ((nx - 1) * (ny - 1)) + 1;+      out[3][count] = types[k][l];+      count++;+    }+  }+  freeMatrix_u(cells, nz - 1);+  freeMatrix_u(types, nz - 1);+  *outnrow = totallength;+  return out;+}++size_t** GetBasic1(unsigned* R, size_t nR, size_t** vivjvk) {+  size_t** indexPtr = malloc(8 * sizeof(size_t*));+  for(unsigned i = 0; i < 8; i++) {+    indexPtr[i] = malloc(3 * sizeof(size_t));+    for(unsigned j = 0; j < 3; j++) {+      indexPtr[i][j] = indexArray[i][j];+    }+  }+  size_t** cube1 = malloc(nR * sizeof(size_t*));+  for(size_t i = 0; i < nR; i++) {+    cube1[i] = malloc(3 * sizeof(size_t));+    size_t* vivjvkRi = vivjvk[R[i]];+    for(unsigned j = 0; j < 3; j++) {+      cube1[i][j] = vivjvkRi[j];+    }+  }+  size_t** k1 = kro1(indexPtr, 8, 3, nR);+  size_t** k2 = kro2(cube1, nR, 3, 8);+  size_t** cubeco = matricialSum(k1, k2, 8 * nR + 1, 3);+  freeMatrix_s(indexPtr, 8);+  freeMatrix_s(cube1, nR);+  freeMatrix_s(k1, 8 * nR + 1);+  freeMatrix_s(k2, 8 * nR + 1);+  return cubeco;+}++double* GetBasic2prime(double*** A, double level, size_t** cubeco, size_t nR) {+  double* values = malloc((8 * nR + 1) * sizeof(double));+  for(size_t i = 0; i < 8 * nR; i++) {+    values[i] = A[cubeco[i][0] - 1][cubeco[i][1] - 1][cubeco[i][2] - 1] - level;+  }+  values[8 * nR] = 0;+  return values;+}++unsigned* get_tcase(unsigned* types, size_t nrow) {+  unsigned CRF[256] = {+      1,  2,  2,  3,  2,  4,  3,  6,  2,  3,  4,  6,  3,  6,  6,  9,  2,  3,+      4,  6,  5,  7,  7,  12, 4,  6,  8,  10, 7,  15, 13, 6,  2,  4,  3,  6,+      4,  8,  6,  10, 5,  7,  7,  15, 7,  13, 12, 6,  3,  6,  6,  9,  7,  13,+      15, 6,  7,  12, 13, 6,  11, 7,  7,  3,  2,  5,  4,  7,  3,  7,  6,  15,+      4,  7,  8,  13, 6,  12, 10, 6,  4,  7,  8,  13, 7,  11, 13, 7,  8,  13,+      14, 8,  13, 7,  8,  4,  3,  7,  6,  12, 6,  13, 9,  6,  7,  11, 13, 7,+      15, 7,  6,  3,  6,  15, 10, 6,  12, 7,  6,  3,  13, 7,  8,  4,  7,  5,+      4,  2,  2,  4,  5,  7,  4,  8,  7,  13, 3,  6,  7,  12, 6,  10, 15, 6,+      3,  6,  7,  15, 7,  13, 11, 7,  6,  9,  13, 6,  12, 6,  7,  3,  4,  8,+      7,  13, 8,  14, 13, 8,  7,  13, 11, 7,  13, 8,  7,  4,  6,  10, 12, 6,+      13, 8,  7,  4,  15, 6,  7,  3,  7,  4,  5,  2,  3,  7,  7,  11, 6,  13,+      12, 7,  6,  15, 13, 7,  9,  6,  6,  3,  6,  12, 13, 7,  15, 7,  7,  5,+      10, 6,  8,  4,  6,  3,  4,  2,  6,  13, 15, 7,  10, 8,  6,  4,  12, 7,+      7,  5,  6,  4,  3,  2,  9,  6,  6,  3,  6,  4,  3,  2,  6,  3,  4,  2,+      3,  2,  2,  1};+  unsigned* out = malloc(nrow * sizeof(unsigned));+  for(size_t i = 0; i < nrow; i++) {+    out[i] = CRF[types[i]] - 1;+  }+  return out;+}++unsigned* getR(unsigned* tcase, size_t nrow, size_t* nR) {+  unsigned* out = malloc(nrow * sizeof(unsigned));+  size_t count = 0;+  for(size_t i = 0; i < nrow; i++) {+    unsigned tc = tcase[i];+    if(tc == 1 || tc == 2 || tc == 5 || tc == 8 || tc == 9 || tc == 11 ||+       tc == 14) {+      out[count] = (unsigned)i;+      count++;+    }+  }+  *nR = count;+  return (out);+}++double** LambdaMu(double* x1, size_t n) {+  double** lambdamu = malloc(2 * sizeof(double*));+  lambdamu[0] = malloc(n * sizeof(double));+  lambdamu[1] = malloc(n * sizeof(double));+  for(size_t i = 0; i < n; i++) {+    lambdamu[0][i] = floor(x1[i] / 9);+    lambdamu[1][i] = 1 - lambdamu[0][i];+  }+  return lambdamu;+}++double* average(double** lambdamu, double* w1, double* w2, size_t n) {+  double* lambda = lambdamu[0];+  double* mu = lambdamu[1];+  double* out = malloc(n * sizeof(double));+  for(size_t i = 0; i < n; i++) {+    out[i] = mu[i] * w1[i] + lambda[i] * w2[i];+  }+  return out;+}++double* average7(double** lambdamu, double* w, size_t n) {+  double* lambda = lambdamu[0];+  double* mu = lambdamu[1];+  double* out = malloc(n * sizeof(double));+  for(size_t i = 0; i < n; i++) {+    out[i] = mu[i] * w[i] + lambda[i];+  }+  return out;+}++double* average8(double** lambdamu, double* w, size_t n) {+  double* lambda = lambdamu[0];+  double* mu = lambdamu[1];+  double* out = malloc(n * sizeof(double));+  for(size_t i = 0; i < n; i++) {+    out[i] = mu[i] * w[i] - lambda[i];+  }+  return out;+}++double** GetPoints(size_t** cubeco,+                   double* values,+                   size_t* p1,+                   unsigned* x1,+                   unsigned* x2,+                   size_t n) {+  unsigned* p1x1 = malloc(n * sizeof(unsigned));+  unsigned* p1x2 = malloc(n * sizeof(unsigned));+  double* xx1 = malloc(n * sizeof(double));+  for(size_t i = 0; i < n; i++) {+    p1x1[i] = p1[i] + x1[i];+    p1x2[i] = p1[i] + x2[i];+    xx1[i] = (double)x1[i];+  }+  double** lambdamu = LambdaMu(xx1, n);+  double* v1 = malloc(n * sizeof(double));+  double* v2 = malloc(n * sizeof(double));+  double* v3 = malloc(n * sizeof(double));+  double* v4 = malloc(n * sizeof(double));+  double* v5 = malloc(n * sizeof(double));+  double* v6 = malloc(n * sizeof(double));+  double* v7 = malloc(n * sizeof(double));+  double* v8 = malloc(n * sizeof(double));+  double* w1 = malloc(n * sizeof(double));+  double* w2 = malloc(n * sizeof(double));+  double* w3 = malloc(n * sizeof(double));+  double* w4 = malloc(n * sizeof(double));+  double* w5 = malloc(n * sizeof(double));+  double* w6 = malloc(n * sizeof(double));+  for(size_t i = 0; i < n; i++) {+    v1[i] = (double)cubeco[p1x1[i] - 2][0];+    w1[i] = (double)cubeco[p1[i] - 1][0];+    v2[i] = (double)cubeco[p1x2[i] - 2][0];+    w2[i] = (double)cubeco[p1[i]][0];+    v3[i] = (double)cubeco[p1x1[i] - 2][1];+    w3[i] = (double)cubeco[p1[i]][1];+    v4[i] = (double)cubeco[p1x2[i] - 2][1];+    w4[i] = (double)cubeco[p1[i] + 1][1];+    v5[i] = (double)cubeco[p1x1[i] - 2][2];+    w5[i] = (double)cubeco[p1[i]][2];+    v6[i] = (double)cubeco[p1x2[i] - 2][2];+    w6[i] = (double)cubeco[p1[i] + 4][2];+    v7[i] = values[p1x1[i] - 2];+    v8[i] = values[p1x2[i] - 2];+  }+  free(p1x1);+  free(p1x2);+  double** out = malloc(8 * sizeof(double*));+  out[0] = average(lambdamu, v1, w1, n);+  free(v1);+  free(w1);+  out[1] = average(lambdamu, v2, w2, n);+  free(v2);+  free(w2);+  out[2] = average(lambdamu, v3, w3, n);+  free(v3);+  free(w3);+  out[3] = average(lambdamu, v4, w4, n);+  free(v4);+  free(w4);+  out[4] = average(lambdamu, v5, w5, n);+  free(v5);+  free(w5);+  out[5] = average(lambdamu, v6, w6, n);+  free(v6);+  free(w6);+  out[6] = average7(lambdamu, v7, n);+  free(v7);+  out[7] = average8(lambdamu, v8, n);+  free(v8);+  freeMatrix_d(lambdamu, 2);+  free(xx1);+  return (out);+}++double** CalPoints(double** points, size_t n) {+  double* x1 = points[0];+  double* x2 = points[1];+  double* y1 = points[2];+  double* y2 = points[3];+  double* z1 = points[4];+  double* z2 = points[5];+  double* v1 = points[6];+  double* v2 = points[7];+  double* x = malloc(n * sizeof(double));+  double* y = malloc(n * sizeof(double));+  double* z = malloc(n * sizeof(double));+  for(size_t i = 0; i < n; i++) {+    double s = v1[i] / (v1[i] - v2[i]);+    x[i] = x1[i] + s * (x2[i] - x1[i]);+    y[i] = y1[i] + s * (y2[i] - y1[i]);+    z[i] = z1[i] + s * (z2[i] - z1[i]);+  }+  double** out = malloc(3 * sizeof(double*));+  out[0] = x;+  out[1] = y;+  out[2] = z;+  double** tout = transpose(out, 3, n);+  freeMatrix_d(out, 3);+  // free(x);free(y);free(z); fait planter+  return (tout);+}++double** computeContour3d(double* voxelAsVector,+                          unsigned nx,+                          unsigned ny,+                          unsigned nz,+                          double max,+                          double level,+                          size_t* ntriangles) {+  printf("Start marching cubes algorithm\n");+  double*** voxel = vector2array(voxelAsVector, nx, ny, nz);+  size_t nrow;+  unsigned** ijkt = levCells(voxel, nx, ny, nz, level, max, &nrow);+  unsigned* tcase = get_tcase(ijkt[3], nrow);+  size_t nR;+  unsigned* R = getR(tcase, nrow, &nR);+  if(nR == 0) {+    *ntriangles = 0;+    free(R);+    free(tcase);+    freeMatrix_u(ijkt, 4);+    return 0;+  } else {+    size_t** vivjvk = malloc(nrow * sizeof(size_t*));+    for(size_t i = 0; i < nrow; i++) {+      vivjvk[i] = malloc(3 * sizeof(size_t));+      for(unsigned j = 0; j < 3; j++) {+        vivjvk[i][j] = (size_t)ijkt[j][i];+      }+    }+    size_t** cubeco = GetBasic1(R, nR, vivjvk);+    double* values = GetBasic2prime(voxel, level, cubeco, nR);+    size_t* p1 = malloc(nR * sizeof(size_t));+    for(size_t i = 0; i < nR; i++) {+      p1[i] = i * 8 + 1;+    }+    unsigned* vt = ijkt[3];+    unsigned* cases = malloc(nR * sizeof(unsigned));+    for(size_t i = 0; i < nR; i++) {+      cases[i] = vt[R[i]] - 1;+    }+    free(R);+    unsigned* edgeslengths = malloc(nR * sizeof(unsigned));+    size_t totalLength = 0;+    for(size_t i = 0; i < nR; i++) {+      edgeslengths[i] = edgesLengths[cases[i]];+      totalLength += (size_t)edgeslengths[i];+    }+    *ntriangles = totalLength;+    size_t* p1rep = replicate(p1, edgeslengths, nR);+    size_t* edges = malloc(totalLength * sizeof(size_t));+    size_t edgeiter = 0;+    for(size_t i = 0; i < nR; i++) {+      for(unsigned j = 0; j < edgeslengths[i]; j++) {+        edges[edgeiter] = (size_t)(*Edges[cases[i]])[j];+        edgeiter++;+      }+    }+    free(edgeslengths);+    unsigned* x1 = malloc(totalLength * sizeof(unsigned));+    unsigned* x2 = malloc(totalLength * sizeof(unsigned));+    for(size_t i = 0; i < totalLength; i++) {+      unsigned* EPi = EdgePoints[edges[i] - 1];+      x1[i] = EPi[1];+      x2[i] = EPi[2];+    }+    free(edges);+    double** points = GetPoints(cubeco, values, p1rep, x1, x2, totalLength);+    double** triangles = CalPoints(points, totalLength);++    unsigned special = 0;+    for(size_t i = 0; i < nrow; i++) {+      unsigned tc = tcase[i];+      if(tc == 3 || tc == 4 || tc == 6 || tc == 7 || tc == 10 || tc == 12 ||+         tc == 13) {+        special = tc;+        break;+      }+    }+    // if(special){+    //   printf("there are special cases\n");+    // }++    for(short c = 0; c < 7; c++) {+      // printf("c: %hi\n", c);+      unsigned nR3;+      unsigned* R3 = whichEqual(tcase, special_name[c], nrow, &nR3);+      // printf("nR3: %u\n", nR3);+      if(nR3 > 0) {  // crash si nR3=2+        size_t** cubeco3 = GetBasic1(R3, nR3, vivjvk);+        double* values3 = GetBasic2prime(voxel, level, cubeco3, nR3);+        size_t* p13 = malloc(nR3 * sizeof(size_t));+        for(size_t i = 0; i < nR3; i++) {+          p13[i] = i * 8 + 1;+        }+        unsigned* cases3 = malloc(nR3 * sizeof(unsigned));+        for(size_t i = 0; i < nR3; i++) {+          cases3[i] = vt[R3[i]] - 1;+        }+        free(R3);+        unsigned nedge = special_nedge[c];+        unsigned outlength3;  // = nR3 ?+        unsigned* index3;+        int* faces3 = unlist(Faces, FacesSizes, cases3, nR3, &outlength3);+        // printf("outlength3: %u\n", outlength3); // yes, outlength3 = nR3 ! ?+        // check+        if(c == 0) {+          index3 = FacesNo7(faces3, p13, values3, nR3, 1);+        } else if(c == 1) {+          index3 = Faces7(faces3, p13, values3, nR3, 1);+        } else {+          unsigned nface = special_nface[c];+          int* facelast = jthColumn(faces3, outlength3, nface, nface - 1);+          index3 = Faces7(facelast, p13, values3, nR3, nface);+          free(facelast);+          for(unsigned j = 0; j < nface - 1; j++) {+            int* facej = jthColumn(faces3, outlength3, nface, j);+            unsigned* temp = FacesNo7(facej, p13, values3, nR3, j + 1);+            index3 = vectorialSum(index3, temp, nR3);+            free(facej);+            free(temp);+          }+        }+        free(faces3);+        unsigned* unlisted_edges3 =+            unlist_u(Edges2, Edges2Sizes, cases3, nR3, &outlength3);+        free(cases3);+        // printf("outlength3 bis: %u\n", outlength3); // outlength3 = nedge *+        // nR3 printf("nedge: %u\n", nedge);+        unsigned** edges3 = vector2matrix(unlisted_edges3, outlength3, nedge);+        free(unlisted_edges3);+        unsigned** edgesp1index = cbind(edges3, p13, index3, nR3, nedge);+        freeMatrix_u(edges3, nR3);+        free(p13);+        unsigned ind3Size = special_indSizes[c];+        unsigned ind3[ind3Size];+        for(unsigned i = 0; i < ind3Size; i++) {+          ind3[i] = (*special_ind[c])[i];+        }+        for(unsigned j = 0; j < ind3Size; j++) {+          unsigned lwrows;+          unsigned* wrows = whichEqual(index3, ind3[j], nR3, &lwrows);+          if(lwrows > 0) {+            unsigned lwcols = *(*(special_posSize)[c])[j] + 1;+            unsigned* wcols = malloc(lwcols * sizeof(unsigned));+            wcols[0] = nedge;+            for(unsigned k = 1; k < lwcols; k++) {+              wcols[k] = (*(*special_pos[c])[j])[k - 1] - 1;+            }+            unsigned** ed =+                subsetMatrix(edgesp1index, wrows, wcols, lwrows, lwcols);+            free(wrows);+            free(wcols);+            size_t* col0ed = malloc(lwrows * sizeof(size_t));+            for(unsigned i = 0; i < lwrows; i++) {+              col0ed[i] = (size_t)ed[i][0];+            }+            size_t* col0edrep = repeach(col0ed, lwcols - 1, lwrows);+            free(col0ed);+            unsigned* edge1 = matrix2vectorMinusFirstColumn(ed, lwrows, lwcols);+            freeMatrix_u(ed, lwrows);+            unsigned totalLength3 = lwrows * (lwcols - 1);+            unsigned* xx1 = malloc(totalLength3 * sizeof(unsigned));+            unsigned* xx2 = malloc(totalLength3 * sizeof(unsigned));+            for(unsigned i = 0; i < totalLength3; i++) {+              unsigned* EPi = EdgePoints[edge1[i] - 1];+              xx1[i] = EPi[1];+              xx2[i] = EPi[2];+            }+            free(edge1);+            double** points3 = GetPoints(cubeco3, values3, col0edrep, xx1, xx2,+                                         (size_t)totalLength3);+            double** triangles3 = CalPoints(points3, totalLength3);+            free(xx1);+            free(xx2);+            free(col0edrep);+            free(points3);+            triangles = realloc(+                triangles, (*ntriangles + totalLength3) * sizeof(*triangles));+            if(triangles == NULL) {+              printf("Error reallocating memory!");+              freeMatrix_d(triangles, totalLength3);+              // free les free de la fin ?+              exit(1);+            }+            for(size_t i = *ntriangles; i < *ntriangles + totalLength3; i++) {+              triangles[i] = malloc(3 * sizeof(double));+              for(short j = 0; j < 3; j++) {+                triangles[i][j] = triangles3[i - *ntriangles][j];+              }+              //*(triangles[i]) = *(triangles3[i - *ntriangles]);+            }+            *ntriangles += totalLength3;+            freeMatrix_d(triangles3, totalLength3);+          }+        } /* end loop for(unsigned j=0; j<ind3Size; j++) */+        free(values3);+      } /* end if(nR3>0) */+    }   /* end for(short i=0; i<7; i++) (loop over special_name) */++    freeMatrix_s(vivjvk, 3);+    freeMatrix_u(ijkt, 4);    // ""+    freeMatrix_d(points, 8);  // peut mettre avant la boucle sur special_name -+                              // non, Pretzel plante+    free(x1);+    free(x2);                 // ""+    freeMatrix_s(cubeco, 3);  // ""+    free(values);             // ""+    free(tcase);+    freeArray(voxel, nx, ny);+    return triangles;+  }+}

@@ -0,0 +1,128 @@+#include <stdlib.h> // for malloc+#include <stdio.h> // for printf+#include <math.h>   // for floor and isnan - compile with flag -lm+++void displayMatrix_u(unsigned**, unsigned, unsigned);++void displayMatrix(double**, unsigned, unsigned);++extern int CaseRotationFlip[256][3];++extern size_t indexArray[8][3];++extern unsigned edgesLengths[254];++extern unsigned EdgePoints[13][3];++extern unsigned (*Edges[])[];++void freeMatrix_u(unsigned**, unsigned);++void freeMatrix_s(size_t**, size_t);++void freeMatrix_d(double**, unsigned);++unsigned** copyMatrix(unsigned**, unsigned, unsigned);++size_t** minorMatrix(size_t**, unsigned, unsigned, unsigned, unsigned);++size_t** scaleMatrix(size_t, size_t**, unsigned, unsigned);++size_t** scaleMinorMatrix(size_t, size_t**, unsigned, unsigned, unsigned, unsigned);++size_t** matricialSum(size_t**, size_t**, size_t, unsigned);++size_t** levelMatrix(double**, unsigned, unsigned, double, unsigned);++double** toMatrix(double***, unsigned, unsigned, unsigned);++unsigned** whichIndicesAndItems(size_t**, unsigned, unsigned, unsigned*);++void freeArray(double***, unsigned, unsigned);++size_t** kro1(size_t**, unsigned, unsigned, size_t);++size_t** kro2(size_t**, size_t, unsigned, unsigned);++size_t* replicate(size_t*, unsigned*, size_t);++double*** vector2array(double*, unsigned, unsigned, unsigned);++size_t** faceType(double**, unsigned, unsigned, double, double);++unsigned** levCells(double***, unsigned, unsigned, unsigned, double, double, size_t*);++size_t** GetBasic1(unsigned*, size_t, size_t**);++double* GetBasic2(double***, double, unsigned*, size_t, size_t**);++unsigned* get_tcase(unsigned*, size_t);++unsigned* getR(unsigned*, size_t, size_t*);++double** LambdaMu(double*, size_t);++double* average(double**, double*, double*, size_t);++double* average7(double**, double*, size_t);++double* average8(double**, double*, size_t);++double** GetPoints(size_t**, double*, size_t*, unsigned*, unsigned*, size_t);++double** CalPoints(double**, size_t);++double** computeContour3d(double*, unsigned, unsigned, unsigned, double, double, size_t*);+++/* Special cases */++extern int (*Faces[])[];++extern unsigned FacesSizes[254];++extern unsigned FacePoints[6][6];++extern unsigned (*Edges2[])[];++extern unsigned Edges2Sizes[254];++extern unsigned special_name[7];+extern unsigned special_nface[7];+extern unsigned special_nedge[7];+extern unsigned (*special_ind[])[];+extern unsigned special_indSizes[7];+extern unsigned (*(*special_pos[])[])[];+extern unsigned* (*(special_posSize)[])[];++int* unlist(int (*jagged[])[], unsigned*, unsigned*, unsigned, unsigned*);+unsigned* unlist_u(unsigned (*jagged[])[], unsigned*, unsigned*, unsigned, unsigned*);++unsigned upow(unsigned, unsigned);++unsigned* FacesNo7(int*, size_t*, double*, size_t, unsigned);++unsigned* Faces7(int*, size_t*, double*, size_t, unsigned);++unsigned** vector2matrix(unsigned*, unsigned, unsigned);++unsigned** cbind(unsigned**, size_t*, unsigned*, unsigned, unsigned);++unsigned** subsetMatrix(unsigned**, unsigned*, unsigned*, unsigned, unsigned);++unsigned* whichEqual(unsigned*, unsigned, unsigned, unsigned*);++//unsigned* matrix2vector(unsigned**, unsigned, unsigned);++unsigned* matrix2vectorMinusFirstColumn(unsigned**, unsigned, unsigned);++size_t* repeach(size_t*, unsigned, size_t);++unsigned* replicatex(unsigned, unsigned);++int* jthColumn(int*, unsigned, unsigned, unsigned);++double** transpose(double**, size_t, size_t);++unsigned* vectorialSum(unsigned*, unsigned*, unsigned);

@@ -0,0 +1,1334 @@+#ifndef __HEADER__+#define __HEADER__+#include "marchingCubes.h"+#endif++int CaseRotationFlip[256][3] =+ {{1, 24, 1},+  {2, 19, 1},+  {2, 17, 1},+  {3, 17, 1},+  {2, 24, 1},+  {4, 24, 1},+  {3, 24, 1},+  {6, 10, 1},+  {2, 15, 1},+  {3, 19, 1},+  {4, 17, 1},+  {6, 9, 1},+  {3, 9, 1},+  {6, 8, 1},+  {6, 1, 1},+  {9, 23, -1},+  {2, 20, 1},+  {3, 18, 1},+  {4, 7, 1},+  {6, 16, 1},+  {5, 24, 1},+  {7, 5, 1},+  {7, 24, 1},+  {12, 9, 1},+  {4, 20, 1},+  {6, 22, 1},+  {8, 24, 1},+  {10, 24, -1},+  {7, 9, 1},+  {15, 24, -1},+  {13, 20, -1},+  {6, 20, -1},+  {2, 21, 1},+  {4, 6, 1},+  {3, 16, 1},+  {6, 4, 1},+  {4, 16, 1},+  {8, 23, 1},+  {6, 14, 1},+  {10, 23, -1},+  {5, 21, 1},+  {7, 10, 1},+  {7, 16, 1},+  {15, 9, 1},+  {7, 2, 1},+  {13, 8, 1},+  {12, 23, -1},+  {6, 6, -1},+  {3, 6, 1},+  {6, 17, 1},+  {6, 18, 1},+  {9, 18, 1},+  {7, 4, 1},+  {13, 17, 1},+  {15, 18, 1},+  {6, 13, -1},+  {7, 6, 1},+  {12, 16, 1},+  {13, 18, 1},+  {6, 2, -1},+  {11, 24, 1},+  {7, 3, -1},+  {7, 12, -1},+  {3, 12, -1},+  {2, 23, 1},+  {5, 23, 1},+  {4, 23, 1},+  {7, 1, 1},+  {3, 14, 1},+  {7, 14, 1},+  {6, 21, 1},+  {15, 23, -1},+  {4, 15, 1},+  {7, 19, 1},+  {8, 19, 1},+  {13, 23, -1},+  {6, 11, 1},+  {12, 17, -1},+  {10, 19, -1},+  {6, 23, -1},+  {4, 12, 1},+  {7, 18, 1},+  {8, 22, 1},+  {13, 16, 1},+  {7, 13, 1},+  {11, 23, 1},+  {13, 21, 1},+  {7, 15, -1},+  {8, 21, 1},+  {13, 22, 1},+  {14, 24, -1},+  {8, 15, -1},+  {13, 11, 1},+  {7, 7, -1},+  {8, 12, -1},+  {4, 22, -1},+  {3, 23, 1},+  {7, 23, 1},+  {6, 24, 1},+  {12, 18, 1},+  {6, 7, 1},+  {13, 19, -1},+  {9, 24, 1},+  {6, 19, -1},+  {7, 21, 1},+  {11, 18, 1},+  {13, 24, 1},+  {7, 20, -1},+  {15, 16, -1},+  {7, 22, -1},+  {6, 15, -1},+  {3, 22, -1},+  {6, 3, 1},+  {15, 17, 1},+  {10, 22, -1},+  {6, 12, -1},+  {12, 24, 1},+  {7, 11, -1},+  {6, 5, -1},+  {3, 15, -1},+  {13, 10, -1},+  {7, 8, -1},+  {8, 20, -1},+  {4, 9, -1},+  {7, 17, -1},+  {5, 22, -1},+  {4, 18, -1},+  {2, 22, -1},+  {2, 22, 1},+  {4, 18, 1},+  {5, 22, 1},+  {7, 17, 1},+  {4, 9, 1},+  {8, 20, 1},+  {7, 8, 1},+  {13, 10, 1},+  {3, 15, 1},+  {6, 5, 1},+  {7, 11, 1},+  {12, 24, -1},+  {6, 12, 1},+  {10, 22, 1},+  {15, 17, -1},+  {6, 3, -1},+  {3, 22, 1},+  {6, 15, 1},+  {7, 22, 1},+  {15, 16, 1},+  {7, 20, 1},+  {13, 24, -1},+  {11, 18, -1},+  {7, 21, -1},+  {6, 19, 1},+  {9, 24, -1},+  {13, 19, 1},+  {6, 7, -1},+  {12, 18, -1},+  {6, 24, -1},+  {7, 23, -1},+  {3, 23, -1},+  {4, 22, 1},+  {8, 12, 1},+  {7, 7, 1},+  {13, 11, -1},+  {8, 15, 1},+  {14, 24, 1},+  {13, 22, -1},+  {8, 21, -1},+  {7, 15, 1},+  {13, 21, -1},+  {11, 23, -1},+  {7, 13, -1},+  {13, 16, -1},+  {8, 22, -1},+  {7, 18, -1},+  {4, 12, -1},+  {6, 23, 1},+  {10, 19, 1},+  {12, 17, 1},+  {6, 11, -1},+  {13, 23, 1},+  {8, 19, -1},+  {7, 19, -1},+  {4, 15, -1},+  {15, 23, 1},+  {6, 21, -1},+  {7, 14, -1},+  {3, 14, -1},+  {7, 1, -1},+  {4, 23, -1},+  {5, 23, -1},+  {2, 23, -1},+  {3, 12, 1},+  {7, 12, 1},+  {7, 3, 1},+  {11, 24, -1},+  {6, 2, 1},+  {13, 18, -1},+  {12, 16, -1},+  {7, 6, -1},+  {6, 13, 1},+  {15, 18, -1},+  {13, 17, -1},+  {7, 4, -1},+  {9, 18, -1},+  {6, 18, -1},+  {6, 17, -1},+  {3, 6, -1},+  {6, 6, 1},+  {12, 23, 1},+  {13, 8, -1},+  {7, 2, -1},+  {15, 9, -1},+  {7, 16, -1},+  {7, 10, -1},+  {5, 21, -1},+  {10, 23, 1},+  {6, 14, -1},+  {8, 23, -1},+  {4, 16, -1},+  {6, 4, -1},+  {3, 16, -1},+  {4, 6, -1},+  {2, 21, -1},+  {6, 20, 1},+  {13, 20, 1},+  {15, 24, 1},+  {7, 9, -1},+  {10, 24, 1},+  {8, 24, -1},+  {6, 22, -1},+  {4, 20, -1},+  {12, 9, -1},+  {7, 24, -1},+  {7, 5, -1},+  {5, 24, -1},+  {6, 16, -1},+  {4, 7, -1},+  {3, 18, -1},+  {2, 20, -1},+  {9, 23, 1},+  {6, 1, -1},+  {6, 8, -1},+  {3, 9, -1},+  {6, 9, -1},+  {4, 17, -1},+  {3, 19, -1},+  {2, 15, -1},+  {6, 10, -1},+  {3, 24, -1},+  {4, 24, -1},+  {2, 24, -1},+  {3, 17, -1},+  {2, 17, -1},+  {2, 19, -1},+  {1, 24, -1}};++size_t indexArray[8][3] =+{{0, 0, 0},+ {1, 0, 0},+ {1, 1, 0},+ {0, 1, 0},+ {0, 0, 1},+ {1, 0, 1},+ {1, 1, 1},+ {0, 1, 1}};++unsigned edgesLengths[254] =+  {3, 3, 6, 3, 0, 6, 9, 3, 6, 0, 9, 6, 9, 9, 6, 3, 6, 0, 9, 0,+   0, 0, 12, 0, 9, 0, 12, 0, 12, 0, 9, 3, 0, 6, 9, 0, 0, 9, 12,+   0, 0, 0, 12, 0, 0, 12, 9, 6, 9, 9, 6, 0, 0, 12, 9, 0, 12, 0,+   9, 0, 0, 0, 6, 3, 0, 0, 0, 6, 0, 9, 12, 0, 0, 0, 0, 9, 12, 12,+   9, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 6, 0, 9, 12,+   9, 0, 6, 9, 0, 0, 0, 0, 12, 0, 9, 6, 9, 12, 12, 9, 12, 0, 9,+   6, 0, 0, 0, 0, 0, 0, 0, 3, 3, 0, 0, 0, 0, 0, 0, 0, 6, 9, 0, 12,+   9, 12, 12, 9, 6, 9, 0, 12, 0, 0, 0, 0, 9, 6, 0, 9, 12, 9, 0,+   6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 9, 12, 12,+   9, 0, 0, 0, 0, 12, 9, 0, 6, 0, 0, 0, 3, 6, 0, 0, 0, 9, 0, 12,+   0, 9, 12, 0, 0, 6, 9, 9, 6, 9, 12, 0, 0, 12, 0, 0, 0, 12, 9,+   0, 0, 9, 6, 0, 3, 9, 0, 12, 0, 12, 0, 9, 0, 12, 0, 0, 0, 9, 0,+   6, 3, 6, 9, 9, 6, 9, 0, 6, 3, 9, 6, 0, 3, 6, 3, 3};++unsigned EdgePoints[13][3] =+ {{1, 1, 2},+  {2, 2, 3},+  {3, 3, 4},+  {4, 4, 1},+  {5, 5, 6},+  {6, 6, 7},+  {7, 7, 8},+  {8, 8, 5},+  {9, 1, 5},+  {10, 2, 6},+  {11, 3, 7},+  {12, 4, 8},+  {13, 9, 9}};+++// https://stackoverflow.com/questions/1083658/do-jagged-arrays-exist-in-c-c+unsigned edgesRow0[] = {4, 9, 1};+unsigned edgesRow1[] = {1, 10, 2};+unsigned edgesRow2[] = {9, 10, 4, 10, 2, 4};+unsigned edgesRow3[] = {2, 11, 3};+unsigned edgesRow4[] = {0};+unsigned edgesRow5[] = {10, 11, 1, 11, 3, 1};+unsigned edgesRow6[] = {4, 9, 3, 3, 9, 11, 9, 10, 11};+unsigned edgesRow7[] = {12, 4, 3};+unsigned edgesRow8[] = {12, 9, 3, 9, 1, 3};+unsigned edgesRow9[] = {0};+unsigned edgesRow10[] = {3, 12, 2, 2, 12, 10, 12, 9, 10};+unsigned edgesRow11[] = {4, 2, 12, 2, 11, 12};+unsigned edgesRow12[] = {2, 11, 1, 1, 11, 9, 11, 12, 9};+unsigned edgesRow13[] = {1, 10, 4, 4, 10, 12, 10, 11, 12};+unsigned edgesRow14[] = {11, 9, 10, 11, 12, 9};+unsigned edgesRow15[] = {8, 5, 9};+unsigned edgesRow16[] = {5, 1, 8, 1, 4, 8};+unsigned edgesRow17[] = {0};+unsigned edgesRow18[] = {10, 2, 5, 5, 2, 8, 2, 4, 8};+unsigned edgesRow19[] = {0};+unsigned edgesRow20[] = {0};+unsigned edgesRow21[] = {0};+unsigned edgesRow22[] = {3, 4, 8, 3, 8, 10, 3, 10, 11, 8, 5, 10};+unsigned edgesRow23[] = {0};+unsigned edgesRow24[] = {8, 5, 12, 12, 5, 3, 5, 1, 3};+unsigned edgesRow25[] = {0};+unsigned edgesRow26[] = {2, 5, 10, 2, 3, 8, 2, 8, 5, 8, 3, 12};+unsigned edgesRow27[] = {0};+unsigned edgesRow28[] = {2, 11, 12, 2, 12, 5, 2, 5, 1, 8, 12, 5};+unsigned edgesRow29[] = {0};+unsigned edgesRow30[] = {8, 5, 12, 5, 10, 12, 12, 10, 11};+unsigned edgesRow31[] = {6, 10, 5};+unsigned edgesRow32[] = {0};+unsigned edgesRow33[] = {1, 5, 2, 5, 6, 2};+unsigned edgesRow34[] = {5, 6, 9, 9, 6, 4, 6, 2, 4};+unsigned edgesRow35[] = {0};+unsigned edgesRow36[] = {0};+unsigned edgesRow37[] = {11, 3, 6, 6, 3, 5, 3, 1, 5};+unsigned edgesRow38[] = {6, 9, 5, 6, 11, 4, 6, 4, 9, 4, 11, 3};+unsigned edgesRow39[] = {0};+unsigned edgesRow40[] = {0};+unsigned edgesRow41[] = {0};+unsigned edgesRow42[] = {3, 6, 2, 3, 9, 6, 3, 12, 9, 5, 9, 6};+unsigned edgesRow43[] = {0};+unsigned edgesRow44[] = {0};+unsigned edgesRow45[] = {6, 1, 5, 6, 12, 1, 6, 11, 12, 1, 12, 4};+unsigned edgesRow46[] = {5, 6, 9, 6, 11, 9, 9, 11, 12};+unsigned edgesRow47[] = {8, 6, 9, 6, 10, 9};+unsigned edgesRow48[] = {1, 4, 10, 10, 4, 6, 4, 8, 6};+unsigned edgesRow49[] = {9, 8, 1, 1, 8, 2, 8, 6, 2};+unsigned edgesRow50[] = {4, 8, 6, 4, 6, 2};+unsigned edgesRow51[] = {0};+unsigned edgesRow52[] = {0};+unsigned edgesRow53[] = {9, 3, 1, 9, 6, 3, 9, 8, 6, 11, 6, 3};+unsigned edgesRow54[] = {11, 3, 6, 3, 4, 6, 6, 4, 8};+unsigned edgesRow55[] = {0};+unsigned edgesRow56[] = {10, 1, 3, 10, 3, 8, 10, 8, 6, 3, 12, 8};+unsigned edgesRow57[] = {0};+unsigned edgesRow58[] = {3, 12, 2, 12, 8, 2, 2, 8, 6};+unsigned edgesRow59[] = {0};+unsigned edgesRow60[] = {0};+unsigned edgesRow61[] = {0};+unsigned edgesRow62[] = {12, 8, 11, 11, 8, 6};+unsigned edgesRow63[] = {6, 7, 11};+unsigned edgesRow64[] = {0};+unsigned edgesRow65[] = {0};+unsigned edgesRow66[] = {0};+unsigned edgesRow67[] = {2, 6, 3, 6, 7, 3};+unsigned edgesRow68[] = {0};+unsigned edgesRow69[] = {6, 7, 10, 10, 7, 1, 7, 3, 1};+unsigned edgesRow70[] = {6, 7, 3, 6, 3, 9, 6, 9, 10, 4, 3, 9};+unsigned edgesRow71[] = {0};+unsigned edgesRow72[] = {0};+unsigned edgesRow73[] = {0};+unsigned edgesRow74[] = {0};+unsigned edgesRow75[] = {12, 4, 7, 7, 4, 6, 4, 2, 6};+unsigned edgesRow76[] = {1, 12, 9, 1, 6, 12, 1, 2, 6, 12, 6, 7};+unsigned edgesRow77[] = {4, 7, 12, 4, 1, 6, 4, 6, 7, 6, 1, 10};+unsigned edgesRow78[] = {6, 7, 10, 7, 12, 10, 10, 12, 9};+unsigned edgesRow79[] = {0};+unsigned edgesRow80[] = {0};+unsigned edgesRow81[] = {0};+unsigned edgesRow82[] = {0};+unsigned edgesRow83[] = {0};+unsigned edgesRow84[] = {0};+unsigned edgesRow85[] = {0};+unsigned edgesRow86[] = {0};+unsigned edgesRow87[] = {0};+unsigned edgesRow88[] = {0};+unsigned edgesRow89[] = {0};+unsigned edgesRow90[] = {0};+unsigned edgesRow91[] = {0};+unsigned edgesRow92[] = {0};+unsigned edgesRow93[] = {0};+unsigned edgesRow94[] = {0};+unsigned edgesRow95[] = {5, 7, 10, 7, 11, 10};+unsigned edgesRow96[] = {0};+unsigned edgesRow97[] = {2, 1, 11, 11, 1, 7, 1, 5, 7};+unsigned edgesRow98[] = {9, 5, 7, 9, 7, 2, 9, 2, 4, 7, 11, 2};+unsigned edgesRow99[] = {10, 5, 2, 2, 5, 3, 5, 7, 3};+unsigned edgesRow100[] = {0};+unsigned edgesRow101[] = {3, 1, 5, 3, 5, 7};+unsigned edgesRow102[] = {4, 9, 3, 9, 5, 3, 3, 5, 7};+unsigned edgesRow103[] = {0};+unsigned edgesRow104[] = {0};+unsigned edgesRow105[] = {0};+unsigned edgesRow106[] = {0};+unsigned edgesRow107[] = {10, 5, 7, 10, 7, 4, 10, 4, 2, 12, 7, 4};+unsigned edgesRow108[] = {0};+unsigned edgesRow109[] = {12, 4, 7, 4, 1, 7, 7, 1, 5};+unsigned edgesRow110[] = {9, 5, 12, 12, 5, 7};+unsigned edgesRow111[] = {7, 11, 8, 8, 11, 9, 11, 10, 9};+unsigned edgesRow112[] = {1, 11, 10, 1, 8, 11, 1, 4, 8, 7, 8, 11};+unsigned edgesRow113[] = {8, 1, 9, 8, 7, 2, 8, 2, 1, 2, 7, 11};+unsigned edgesRow114[] = {7, 11, 8, 11, 2, 8, 8, 2, 4};+unsigned edgesRow115[] = {2, 10, 9, 2, 9, 7, 2, 7, 3, 9, 8, 7};+unsigned edgesRow116[] = {0};+unsigned edgesRow117[] = {9, 8, 1, 8, 7, 1, 1, 7, 3};+unsigned edgesRow118[] = {8, 7, 4, 4, 7, 3};+unsigned edgesRow119[] = {0};+unsigned edgesRow120[] = {0};+unsigned edgesRow121[] = {0};+unsigned edgesRow122[] = {0};+unsigned edgesRow123[] = {0};+unsigned edgesRow124[] = {0};+unsigned edgesRow125[] = {0};+unsigned edgesRow126[] = {8, 7, 12};+unsigned edgesRow127[] = {8, 12, 7};+unsigned edgesRow128[] = {0};+unsigned edgesRow129[] = {0};+unsigned edgesRow130[] = {0};+unsigned edgesRow131[] = {0};+unsigned edgesRow132[] = {0};+unsigned edgesRow133[] = {0};+unsigned edgesRow134[] = {0};+unsigned edgesRow135[] = {8, 4, 7, 4, 3, 7};+unsigned edgesRow136[] = {9, 1, 8, 8, 1, 7, 1, 3, 7};+unsigned edgesRow137[] = {0};+unsigned edgesRow138[] = {2, 9, 10, 2, 7, 9, 2, 3, 7, 9, 7, 8};+unsigned edgesRow139[] = {7, 8, 11, 11, 8, 2, 8, 4, 2};+unsigned edgesRow140[] = {8, 9, 1, 8, 2, 7, 8, 1, 2, 2, 11, 7};+unsigned edgesRow141[] = {1, 10, 11, 1, 11, 8, 1, 8, 4, 7, 11, 8};+unsigned edgesRow142[] = {7, 8, 11, 8, 9, 11, 11, 9, 10};+unsigned edgesRow143[] = {9, 12, 5, 12, 7, 5};+unsigned edgesRow144[] = {12, 7, 4, 4, 7, 1, 7, 5, 1};+unsigned edgesRow145[] = {0};+unsigned edgesRow146[] = {10, 7, 5, 10, 4, 7, 10, 2, 4, 12, 4, 7};+unsigned edgesRow147[] = {0};+unsigned edgesRow148[] = {0};+unsigned edgesRow149[] = {0};+unsigned edgesRow150[] = {0};+unsigned edgesRow151[] = {4, 3, 9, 9, 3, 5, 3, 7, 5};+unsigned edgesRow152[] = {3, 5, 1, 3, 7, 5};+unsigned edgesRow153[] = {0};+unsigned edgesRow154[] = {10, 2, 5, 2, 3, 5, 5, 3, 7};+unsigned edgesRow155[] = {9, 7, 5, 9, 2, 7, 9, 4, 2, 7, 2, 11};+unsigned edgesRow156[] = {2, 11, 1, 11, 7, 1, 1, 7, 5};+unsigned edgesRow157[] = {0};+unsigned edgesRow158[] = {5, 10, 7, 7, 10, 11};+unsigned edgesRow159[] = {0};+unsigned edgesRow160[] = {0};+unsigned edgesRow161[] = {0};+unsigned edgesRow162[] = {0};+unsigned edgesRow163[] = {0};+unsigned edgesRow164[] = {0};+unsigned edgesRow165[] = {0};+unsigned edgesRow166[] = {0};+unsigned edgesRow167[] = {0};+unsigned edgesRow168[] = {0};+unsigned edgesRow169[] = {0};+unsigned edgesRow170[] = {0};+unsigned edgesRow171[] = {0};+unsigned edgesRow172[] = {0};+unsigned edgesRow173[] = {0};+unsigned edgesRow174[] = {0};+unsigned edgesRow175[] = {6, 10, 7, 7, 10, 12, 10, 9, 12};+unsigned edgesRow176[] = {4, 12, 7, 4, 6, 1, 4, 7, 6, 6, 10, 1};+unsigned edgesRow177[] = {1, 9, 12, 1, 12, 6, 1, 6, 2, 12, 7, 6};+unsigned edgesRow178[] = {12, 7, 4, 7, 6, 4, 4, 6, 2};+unsigned edgesRow179[] = {0};+unsigned edgesRow180[] = {0};+unsigned edgesRow181[] = {0};+unsigned edgesRow182[] = {0};+unsigned edgesRow183[] = {6, 3, 7, 6, 9, 3, 6, 10, 9, 4, 9, 3};+unsigned edgesRow184[] = {6, 10, 7, 10, 1, 7, 7, 1, 3};+unsigned edgesRow185[] = {0};+unsigned edgesRow186[] = {2, 3, 6, 6, 3, 7};+unsigned edgesRow187[] = {0};+unsigned edgesRow188[] = {0};+unsigned edgesRow189[] = {0};+unsigned edgesRow190[] = {6, 11, 7};+unsigned edgesRow191[] = {12, 11, 8, 11, 6, 8};+unsigned edgesRow192[] = {0};+unsigned edgesRow193[] = {0};+unsigned edgesRow194[] = {0};+unsigned edgesRow195[] = {3, 2, 12, 12, 2, 8, 2, 6, 8};+unsigned edgesRow196[] = {0};+unsigned edgesRow197[] = {10, 3, 1, 10, 8, 3, 10, 6, 8, 3, 8, 12};+unsigned edgesRow198[] = {0};+unsigned edgesRow199[] = {11, 6, 3, 3, 6, 4, 6, 8, 4};+unsigned edgesRow200[] = {9, 1, 3, 9, 3, 6, 9, 6, 8, 11, 3, 6};+unsigned edgesRow201[] = {0};+unsigned edgesRow202[] = {0};+unsigned edgesRow203[] = {4, 6, 8, 4, 2, 6};+unsigned edgesRow204[] = {9, 1, 8, 1, 2, 8, 8, 2, 6};+unsigned edgesRow205[] = {1, 10, 4, 10, 6, 4, 4, 6, 8};+unsigned edgesRow206[] = {8, 9, 6, 6, 9, 10};+unsigned edgesRow207[] = {5, 9, 6, 6, 9, 11, 9, 12, 11};+unsigned edgesRow208[] = {6, 5, 1, 6, 1, 12, 6, 12, 11, 1, 4, 12};+unsigned edgesRow209[] = {0};+unsigned edgesRow210[] = {0};+unsigned edgesRow211[] = {3, 2, 6, 3, 6, 9, 3, 9, 12, 5, 6, 9};+unsigned edgesRow212[] = {0};+unsigned edgesRow213[] = {0};+unsigned edgesRow214[] = {0};+unsigned edgesRow215[] = {6, 5, 9, 6, 4, 11, 6, 9, 4, 4, 3, 11};+unsigned edgesRow216[] = {11, 6, 3, 6, 5, 3, 3, 5, 1};+unsigned edgesRow217[] = {0};+unsigned edgesRow218[] = {0};+unsigned edgesRow219[] = {5, 9, 6, 9, 4, 6, 6, 4, 2};+unsigned edgesRow220[] = {1, 2, 5, 5, 2, 6};+unsigned edgesRow221[] = {0};+unsigned edgesRow222[] = {6, 5, 10};+unsigned edgesRow223[] = {8, 12, 5, 5, 12, 10, 12, 11, 10};+unsigned edgesRow224[] = {0};+unsigned edgesRow225[] = {2, 12, 11, 2, 5, 12, 2, 1, 5, 8, 5, 12};+unsigned edgesRow226[] = {0};+unsigned edgesRow227[] = {2, 10, 5, 2, 8, 3, 2, 5, 8, 8, 12, 3};+unsigned edgesRow228[] = {0};+unsigned edgesRow229[] = {8, 12, 5, 12, 3, 5, 5, 3, 1};+unsigned edgesRow230[] = {0};+unsigned edgesRow231[] = {3, 8, 4, 3, 10, 8, 3, 11, 10, 8, 10, 5};+unsigned edgesRow232[] = {0};+unsigned edgesRow233[] = {0};+unsigned edgesRow234[] = {0};+unsigned edgesRow235[] = {10, 5, 2, 5, 8, 2, 2, 8, 4};+unsigned edgesRow236[] = {0};+unsigned edgesRow237[] = {5, 8, 1, 1, 8, 4};+unsigned edgesRow238[] = {8, 9, 5};+unsigned edgesRow239[] = {11, 10, 9, 11, 9, 12};+unsigned edgesRow240[] = {1, 4, 10, 4, 12, 10, 10, 12, 11};+unsigned edgesRow241[] = {2, 1, 11, 1, 9, 11, 11, 9, 12};+unsigned edgesRow242[] = {4, 12, 2, 2, 12, 11};+unsigned edgesRow243[] = {3, 2, 12, 2, 10, 12, 12, 10, 9};+unsigned edgesRow244[] = {0};+unsigned edgesRow245[] = {12, 3, 9, 9, 3, 1};+unsigned edgesRow246[] = {12, 3, 4};+unsigned edgesRow247[] = {4, 3, 9, 3, 11, 9, 9, 11, 10};+unsigned edgesRow248[] = {10, 1, 11, 11, 1, 3};+unsigned edgesRow249[] = {0};+unsigned edgesRow250[] = {2, 3, 11};+unsigned edgesRow251[] = {9, 4, 10, 10, 4, 2};+unsigned edgesRow252[] = {1, 2, 10};+unsigned edgesRow253[] = {4, 1, 9};+unsigned (*Edges[])[] = {&edgesRow0, &edgesRow1, &edgesRow2, &edgesRow3, &edgesRow4, &edgesRow5, &edgesRow6, &edgesRow7, &edgesRow8, &edgesRow9, &edgesRow10, &edgesRow11, &edgesRow12, &edgesRow13, &edgesRow14, &edgesRow15, &edgesRow16, &edgesRow17, &edgesRow18, &edgesRow19, &edgesRow20, &edgesRow21, &edgesRow22, &edgesRow23, &edgesRow24, &edgesRow25, &edgesRow26, &edgesRow27, &edgesRow28, &edgesRow29, &edgesRow30, &edgesRow31, &edgesRow32, &edgesRow33, &edgesRow34, &edgesRow35, &edgesRow36,+   &edgesRow37, &edgesRow38, &edgesRow39, &edgesRow40, &edgesRow41, &edgesRow42, &edgesRow43, &edgesRow44, &edgesRow45, &edgesRow46, &edgesRow47, &edgesRow48, &edgesRow49, &edgesRow50, &edgesRow51, &edgesRow52, &edgesRow53, &edgesRow54, &edgesRow55, &edgesRow56, &edgesRow57, &edgesRow58, &edgesRow59, &edgesRow60, &edgesRow61, &edgesRow62, &edgesRow63, &edgesRow64, &edgesRow65, &edgesRow66, &edgesRow67, &edgesRow68, &edgesRow69, &edgesRow70, &edgesRow71, &edgesRow72, &edgesRow73, &edgesRow74,+   &edgesRow75, &edgesRow76, &edgesRow77, &edgesRow78, &edgesRow79, &edgesRow80, &edgesRow81, &edgesRow82, &edgesRow83, &edgesRow84, &edgesRow85, &edgesRow86, &edgesRow87, &edgesRow88, &edgesRow89, &edgesRow90, &edgesRow91, &edgesRow92, &edgesRow93, &edgesRow94, &edgesRow95, &edgesRow96, &edgesRow97, &edgesRow98, &edgesRow99, &edgesRow100, &edgesRow101, &edgesRow102, &edgesRow103, &edgesRow104, &edgesRow105, &edgesRow106, &edgesRow107, &edgesRow108, &edgesRow109, &edgesRow110, &edgesRow111,+   &edgesRow112, &edgesRow113, &edgesRow114, &edgesRow115, &edgesRow116, &edgesRow117, &edgesRow118, &edgesRow119, &edgesRow120, &edgesRow121, &edgesRow122, &edgesRow123, &edgesRow124, &edgesRow125, &edgesRow126, &edgesRow127, &edgesRow128, &edgesRow129, &edgesRow130, &edgesRow131, &edgesRow132, &edgesRow133, &edgesRow134, &edgesRow135, &edgesRow136, &edgesRow137, &edgesRow138, &edgesRow139, &edgesRow140, &edgesRow141, &edgesRow142, &edgesRow143, &edgesRow144, &edgesRow145, &edgesRow146,+   &edgesRow147, &edgesRow148, &edgesRow149, &edgesRow150, &edgesRow151, &edgesRow152, &edgesRow153, &edgesRow154, &edgesRow155, &edgesRow156, &edgesRow157, &edgesRow158, &edgesRow159, &edgesRow160, &edgesRow161, &edgesRow162, &edgesRow163, &edgesRow164, &edgesRow165, &edgesRow166, &edgesRow167, &edgesRow168, &edgesRow169, &edgesRow170, &edgesRow171, &edgesRow172, &edgesRow173, &edgesRow174, &edgesRow175, &edgesRow176, &edgesRow177, &edgesRow178, &edgesRow179, &edgesRow180, &edgesRow181,+   &edgesRow182, &edgesRow183, &edgesRow184, &edgesRow185, &edgesRow186, &edgesRow187, &edgesRow188, &edgesRow189, &edgesRow190, &edgesRow191, &edgesRow192, &edgesRow193, &edgesRow194, &edgesRow195, &edgesRow196, &edgesRow197, &edgesRow198, &edgesRow199, &edgesRow200, &edgesRow201, &edgesRow202, &edgesRow203, &edgesRow204, &edgesRow205, &edgesRow206, &edgesRow207, &edgesRow208, &edgesRow209, &edgesRow210, &edgesRow211, &edgesRow212, &edgesRow213, &edgesRow214, &edgesRow215, &edgesRow216,+   &edgesRow217, &edgesRow218, &edgesRow219, &edgesRow220, &edgesRow221, &edgesRow222, &edgesRow223, &edgesRow224, &edgesRow225, &edgesRow226, &edgesRow227, &edgesRow228, &edgesRow229, &edgesRow230, &edgesRow231, &edgesRow232, &edgesRow233, &edgesRow234, &edgesRow235, &edgesRow236, &edgesRow237, &edgesRow238, &edgesRow239, &edgesRow240, &edgesRow241, &edgesRow242, &edgesRow243, &edgesRow244, &edgesRow245, &edgesRow246, &edgesRow247, &edgesRow248, &edgesRow249, &edgesRow250, &edgesRow251,+   &edgesRow252, &edgesRow253};++/* Special cases */+int facesRow0[] = {0};+int facesRow1[] = {0};+int facesRow2[] = {0};+int facesRow3[] = {0};+int facesRow4[] = {5};+int facesRow5[] = {0};+int facesRow6[] = {0};+int facesRow7[] = {0};+int facesRow8[] = {0};+int facesRow9[] = {-5};+int facesRow10[] = {0};+int facesRow11[] = {0};+int facesRow12[] = {0};+int facesRow13[] = {0};+int facesRow14[] = {0};+int facesRow15[] = {0};+int facesRow16[] = {0};+int facesRow17[] = {-1};+int facesRow18[] = {0};+int facesRow19[] = {7};+int facesRow20[] = {5, 7};+int facesRow21[] = {-1, 7};+int facesRow22[] = {0};+int facesRow23[] = {-4};+int facesRow24[] = {0};+int facesRow25[] = {-5, -1, -4, 7};+int facesRow26[] = {0};+int facesRow27[] = {-4, 7};+int facesRow28[] = {0};+int facesRow29[] = {-4, -1, 7};+int facesRow30[] = {0};+int facesRow31[] = {0};+int facesRow32[] = {1};+int facesRow33[] = {0};+int facesRow34[] = {0};+int facesRow35[] = {-2};+int facesRow36[] = {-2, 1, 5, -7};+int facesRow37[] = {0};+int facesRow38[] = {0};+int facesRow39[] = {-7};+int facesRow40[] = {1, -7};+int facesRow41[] = {-5, -7};+int facesRow42[] = {0};+int facesRow43[] = {-2, -7};+int facesRow44[] = {-2, 1, -7};+int facesRow45[] = {0};+int facesRow46[] = {0};+int facesRow47[] = {0};+int facesRow48[] = {0};+int facesRow49[] = {0};+int facesRow50[] = {0};+int facesRow51[] = {-2, 7};+int facesRow52[] = {5, -2, 7};+int facesRow53[] = {0};+int facesRow54[] = {0};+int facesRow55[] = {-4, -7};+int facesRow56[] = {0};+int facesRow57[] = {-4, -5, -7};+int facesRow58[] = {0};+int facesRow59[] = {-2, -4, 7};+int facesRow60[] = {2, -7};+int facesRow61[] = {4, 7};+int facesRow62[] = {0};+int facesRow63[] = {0};+int facesRow64[] = {7};+int facesRow65[] = {2};+int facesRow66[] = {2, 7};+int facesRow67[] = {0};+int facesRow68[] = {5, 7};+int facesRow69[] = {0};+int facesRow70[] = {0};+int facesRow71[] = {3};+int facesRow72[] = {3, 7};+int facesRow73[] = {-5, 3, 2, 7};+int facesRow74[] = {2, 3, 7};+int facesRow75[] = {0};+int facesRow76[] = {0};+int facesRow77[] = {0};+int facesRow78[] = {0};+int facesRow79[] = {6};+int facesRow80[] = {6, 7};+int facesRow81[] = {6, -1, 2, -7};+int facesRow82[] = {2, 6, 7};+int facesRow83[] = {6, 7};+int facesRow84[] = {6, 5, 7};+int facesRow85[] = {6, -1, 7};+int facesRow86[] = {-6, -7};+int facesRow87[] = {6, 3, -4, -7};+int facesRow88[] = {6, 3, 7};+int facesRow89[] = {-5, -1, 6, 3, 2, -4, 7};+int facesRow90[] = {-3, -6, -2, 7};+int facesRow91[] = {-4, 6, 7};+int facesRow92[] = {-6, -7};+int facesRow93[] = {-6, 4, 1, 7};+int facesRow94[] = {-6};+int facesRow95[] = {0};+int facesRow96[] = {1, 7};+int facesRow97[] = {0};+int facesRow98[] = {0};+int facesRow99[] = {0};+int facesRow100[] = {5, 1, 7};+int facesRow101[] = {0};+int facesRow102[] = {0};+int facesRow103[] = {3, -7};+int facesRow104[] = {1, 3, 7};+int facesRow105[] = {-5, 3, -7};+int facesRow106[] = {-3, 7};+int facesRow107[] = {0};+int facesRow108[] = {-1, -7};+int facesRow109[] = {0};+int facesRow110[] = {0};+int facesRow111[] = {0};+int facesRow112[] = {0};+int facesRow113[] = {0};+int facesRow114[] = {0};+int facesRow115[] = {0};+int facesRow116[] = {-5, -7};+int facesRow117[] = {0};+int facesRow118[] = {0};+int facesRow119[] = {-4, 3, -7};+int facesRow120[] = {-3, -7};+int facesRow121[] = {4, -3, 5, -7};+int facesRow122[] = {-3};+int facesRow123[] = {4, -7};+int facesRow124[] = {-7};+int facesRow125[] = {4};+int facesRow126[] = {0};+int facesRow127[] = {0};+int facesRow128[] = {4};+int facesRow129[] = {-7};+int facesRow130[] = {4, -7};+int facesRow131[] = {-3};+int facesRow132[] = {4, -3, 5, -7};+int facesRow133[] = {-3, -7};+int facesRow134[] = {4, -3, -7};+int facesRow135[] = {0};+int facesRow136[] = {0};+int facesRow137[] = {-5, -7};+int facesRow138[] = {0};+int facesRow139[] = {0};+int facesRow140[] = {0};+int facesRow141[] = {0};+int facesRow142[] = {0};+int facesRow143[] = {0};+int facesRow144[] = {0};+int facesRow145[] = {-1, -7};+int facesRow146[] = {0};+int facesRow147[] = {-3, 7};+int facesRow148[] = {5, -3, 7};+int facesRow149[] = {-1, -3, 7};+int facesRow150[] = {3, -7};+int facesRow151[] = {0};+int facesRow152[] = {0};+int facesRow153[] = {-5, -1, -7};+int facesRow154[] = {0};+int facesRow155[] = {0};+int facesRow156[] = {0};+int facesRow157[] = {1, 7};+int facesRow158[] = {0};+int facesRow159[] = {-6};+int facesRow160[] = {-6, 4, 1, 7};+int facesRow161[] = {-6, -7};+int facesRow162[] = {4, -6, -7};+int facesRow163[] = {-3, -6, -2, 7};+int facesRow164[] = {5, 1, -6, -3, -2, 4, 7};+int facesRow165[] = {-6, -3, -7};+int facesRow166[] = {6, 3, -4, -7};+int facesRow167[] = {-6, -7};+int facesRow168[] = {-6, 1, -7};+int facesRow169[] = {-6, -5, -7};+int facesRow170[] = {6, 7};+int facesRow171[] = {-2, -6, -7};+int facesRow172[] = {6, -1, 2, -7};+int facesRow173[] = {6, 7};+int facesRow174[] = {6};+int facesRow175[] = {0};+int facesRow176[] = {0};+int facesRow177[] = {0};+int facesRow178[] = {0};+int facesRow179[] = {-2, -3, 7};+int facesRow180[] = {-5, 3, 2, 7};+int facesRow181[] = {3, 7};+int facesRow182[] = {3};+int facesRow183[] = {0};+int facesRow184[] = {0};+int facesRow185[] = {5, 7};+int facesRow186[] = {0};+int facesRow187[] = {2, 7};+int facesRow188[] = {2};+int facesRow189[] = {7};+int facesRow190[] = {0};+int facesRow191[] = {0};+int facesRow192[] = {4, 7};+int facesRow193[] = {2, -7};+int facesRow194[] = {2, 4, 7};+int facesRow195[] = {0};+int facesRow196[] = {4, 5, 7};+int facesRow197[] = {0};+int facesRow198[] = {-4, -7};+int facesRow199[] = {0};+int facesRow200[] = {0};+int facesRow201[] = {-5, 2, -7};+int facesRow202[] = {-2, 7};+int facesRow203[] = {0};+int facesRow204[] = {0};+int facesRow205[] = {0};+int facesRow206[] = {0};+int facesRow207[] = {0};+int facesRow208[] = {0};+int facesRow209[] = {2, -1, -7};+int facesRow210[] = {-2, -7};+int facesRow211[] = {0};+int facesRow212[] = {-5, -7};+int facesRow213[] = {1, -7};+int facesRow214[] = {-7};+int facesRow215[] = {0};+int facesRow216[] = {0};+int facesRow217[] = {-2, 1, 5, -7};+int facesRow218[] = {-2};+int facesRow219[] = {0};+int facesRow220[] = {0};+int facesRow221[] = {1};+int facesRow222[] = {0};+int facesRow223[] = {0};+int facesRow224[] = {4, 1, 7};+int facesRow225[] = {0};+int facesRow226[] = {-4, 7};+int facesRow227[] = {0};+int facesRow228[] = {-5, -1, -4, 7};+int facesRow229[] = {0};+int facesRow230[] = {-4};+int facesRow231[] = {0};+int facesRow232[] = {-1, 7};+int facesRow233[] = {5, 7};+int facesRow234[] = {7};+int facesRow235[] = {0};+int facesRow236[] = {-1};+int facesRow237[] = {0};+int facesRow238[] = {0};+int facesRow239[] = {0};+int facesRow240[] = {0};+int facesRow241[] = {0};+int facesRow242[] = {0};+int facesRow243[] = {0};+int facesRow244[] = {-5};+int facesRow245[] = {0};+int facesRow246[] = {0};+int facesRow247[] = {0};+int facesRow248[] = {0};+int facesRow249[] = {5};+int facesRow250[] = {0};+int facesRow251[] = {0};+int facesRow252[] = {0};+int facesRow253[] = {0};+int (*Faces[])[] = {&facesRow0, &facesRow1, &facesRow2, &facesRow3, &facesRow4, &facesRow5, &facesRow6, &facesRow7, &facesRow8, &facesRow9, &facesRow10, &facesRow11, &facesRow12, &facesRow13, &facesRow14, &facesRow15, &facesRow16, &facesRow17, &facesRow18, &facesRow19, &facesRow20, &facesRow21, &facesRow22, &facesRow23, &facesRow24, &facesRow25, &facesRow26, &facesRow27, &facesRow28, &facesRow29, &facesRow30, &facesRow31, &facesRow32, &facesRow33, &facesRow34, &facesRow35, &facesRow36,+   &facesRow37, &facesRow38, &facesRow39, &facesRow40, &facesRow41, &facesRow42, &facesRow43, &facesRow44, &facesRow45, &facesRow46, &facesRow47, &facesRow48, &facesRow49, &facesRow50, &facesRow51, &facesRow52, &facesRow53, &facesRow54, &facesRow55, &facesRow56, &facesRow57, &facesRow58, &facesRow59, &facesRow60, &facesRow61, &facesRow62, &facesRow63, &facesRow64, &facesRow65, &facesRow66, &facesRow67, &facesRow68, &facesRow69, &facesRow70, &facesRow71, &facesRow72, &facesRow73, &facesRow74,+   &facesRow75, &facesRow76, &facesRow77, &facesRow78, &facesRow79, &facesRow80, &facesRow81, &facesRow82, &facesRow83, &facesRow84, &facesRow85, &facesRow86, &facesRow87, &facesRow88, &facesRow89, &facesRow90, &facesRow91, &facesRow92, &facesRow93, &facesRow94, &facesRow95, &facesRow96, &facesRow97, &facesRow98, &facesRow99, &facesRow100, &facesRow101, &facesRow102, &facesRow103, &facesRow104, &facesRow105, &facesRow106, &facesRow107, &facesRow108, &facesRow109, &facesRow110, &facesRow111,+   &facesRow112, &facesRow113, &facesRow114, &facesRow115, &facesRow116, &facesRow117, &facesRow118, &facesRow119, &facesRow120, &facesRow121, &facesRow122, &facesRow123, &facesRow124, &facesRow125, &facesRow126, &facesRow127, &facesRow128, &facesRow129, &facesRow130, &facesRow131, &facesRow132, &facesRow133, &facesRow134, &facesRow135, &facesRow136, &facesRow137, &facesRow138, &facesRow139, &facesRow140, &facesRow141, &facesRow142, &facesRow143, &facesRow144, &facesRow145, &facesRow146,+   &facesRow147, &facesRow148, &facesRow149, &facesRow150, &facesRow151, &facesRow152, &facesRow153, &facesRow154, &facesRow155, &facesRow156, &facesRow157, &facesRow158, &facesRow159, &facesRow160, &facesRow161, &facesRow162, &facesRow163, &facesRow164, &facesRow165, &facesRow166, &facesRow167, &facesRow168, &facesRow169, &facesRow170, &facesRow171, &facesRow172, &facesRow173, &facesRow174, &facesRow175, &facesRow176, &facesRow177, &facesRow178, &facesRow179, &facesRow180, &facesRow181,+   &facesRow182, &facesRow183, &facesRow184, &facesRow185, &facesRow186, &facesRow187, &facesRow188, &facesRow189, &facesRow190, &facesRow191, &facesRow192, &facesRow193, &facesRow194, &facesRow195, &facesRow196, &facesRow197, &facesRow198, &facesRow199, &facesRow200, &facesRow201, &facesRow202, &facesRow203, &facesRow204, &facesRow205, &facesRow206, &facesRow207, &facesRow208, &facesRow209, &facesRow210, &facesRow211, &facesRow212, &facesRow213, &facesRow214, &facesRow215, &facesRow216,+   &facesRow217, &facesRow218, &facesRow219, &facesRow220, &facesRow221, &facesRow222, &facesRow223, &facesRow224, &facesRow225, &facesRow226, &facesRow227, &facesRow228, &facesRow229, &facesRow230, &facesRow231, &facesRow232, &facesRow233, &facesRow234, &facesRow235, &facesRow236, &facesRow237, &facesRow238, &facesRow239, &facesRow240, &facesRow241, &facesRow242, &facesRow243, &facesRow244, &facesRow245, &facesRow246, &facesRow247, &facesRow248, &facesRow249, &facesRow250, &facesRow251,+   &facesRow252, &facesRow253};++unsigned FacesSizes[254] = {1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,+                            2, 2, 1, 1, 1, 4, 1, 2, 1, 3, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 2,+                            2, 1, 2, 3, 1, 1, 1, 1, 1, 1, 2, 3, 1, 1, 2, 1, 3, 1, 3, 2, 2,+                            1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 2, 4, 3, 1, 1, 1, 1, 1, 2, 4, 3,+                            2, 3, 3, 2, 4, 3, 7, 4, 3, 2, 4, 1, 1, 2, 1, 1, 1, 3, 1, 1, 2,+                            3, 3, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 3, 2, 4, 1, 2, 1,+                            1, 1, 1, 1, 1, 2, 1, 4, 2, 3, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2,+                            1, 2, 3, 3, 2, 1, 1, 3, 1, 1, 1, 2, 1, 1, 4, 2, 3, 4, 7, 3, 4,+                            2, 3, 3, 2, 3, 4, 2, 1, 1, 1, 1, 1, 3, 4, 2, 1, 1, 1, 2, 1, 2,+                            1, 1, 1, 1, 2, 2, 3, 1, 3, 1, 2, 1, 1, 3, 2, 1, 1, 1, 1, 1, 1,+                            3, 2, 1, 2, 2, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 3, 1, 2, 1, 4, 1,+                            1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,+                            1, 1, 1};++unsigned FacePoints[6][6] =+  {{1, 1, 6, 2, 5, 60},+   {2, 2, 7, 3, 6, 252},+   {3, 4, 7, 3, 8, 672},+   {4, 1, 8, 4, 5, 160},+   {5, 1, 3, 2, 4, 24},+   {6, 5, 7, 6, 8, 1680}};++unsigned edges2Row0[] = {0};+unsigned edges2Row1[] = {0};+unsigned edges2Row2[] = {0};+unsigned edges2Row3[] = {0};+unsigned edges2Row4[] = {2, 11, 3, 4, 9, 1, 2, 11, 1, 11, 9, 1, 11, 4, 9, 11, 3, 4};+unsigned edges2Row5[] = {0};+unsigned edges2Row6[] = {0};+unsigned edges2Row7[] = {0};+unsigned edges2Row8[] = {0};+unsigned edges2Row9[] = {1, 10, 2, 3, 12, 4, 1, 10, 4, 10, 12, 4, 10, 3, 12, 10, 2, 3};+unsigned edges2Row10[] = {0};+unsigned edges2Row11[] = {0};+unsigned edges2Row12[] = {0};+unsigned edges2Row13[] = {0};+unsigned edges2Row14[] = {0};+unsigned edges2Row15[] = {0};+unsigned edges2Row16[] = {0};+unsigned edges2Row17[] = {10, 2, 1, 9, 8, 5, 10, 2, 5, 2, 8, 5, 2, 9, 8, 2, 1, 9};+unsigned edges2Row18[] = {0};+unsigned edges2Row19[] = {2, 11, 3, 9, 8, 5, 2, 9, 3, 9, 8, 3, 11, 3, 8, 2, 11, 5, 2, 5, 9, 11, 8, 5};+unsigned edges2Row20[] = {4, 8, 1, 8, 5, 1, 2, 11, 3, 4, 8, 3, 8, 11, 3, 8, 5, 11, 2, 11, 5, 2, 5, 1, 4, 1, 2, 4, 2, 3, 3, 2, 11, 4, 8, 3, 8, 11, 3, 8, 5, 11, 2, 11, 5, 2, 5, 1};+unsigned edges2Row21[] = {10, 11, 1, 11, 3, 1, 9, 8, 5, 10, 11, 5, 11, 8, 5, 11, 3, 8, 9, 8, 3, 9, 3, 1, 10, 1, 9, 10, 9, 5, 5, 9, 8, 10, 11, 5, 11, 8, 5, 11, 3, 8, 9, 8, 3, 9, 3, 1};+unsigned edges2Row22[] = {0};+unsigned edges2Row23[] = {8, 5, 9, 4, 3, 12, 8, 5, 12, 5, 3, 12, 5, 4, 3, 5, 9, 4};+unsigned edges2Row24[] = {0};+unsigned edges2Row25[] = {2, 1, 10, 4, 3, 12, 9, 8, 5, 2, 1, 10, 4, 5, 9, 4, 3, 5, 8, 5, 3, 8, 3, 12, 4, 3, 12, 2, 1, 9, 2, 9, 8, 2, 5, 10, 2, 8, 5, 9, 8, 5, 2, 3, 12, 2, 12, 10, 10, 12, 4, 10, 4, 1, 2, 13, 10, 10, 13, 5, 8, 5, 13, 8, 13, 12, 12, 13, 3, 4, 3, 13, 4, 13, 9, 9, 13, 1, 2, 1, 13, 2, 13, 10, 2, 13, 3, 12, 13, 3, 8, 13, 12, 8, 5, 13, 9, 5, 13, 4, 13, 9, 4, 1, 13, 10, 13, 1, 8, 5, 13, 10, 13, 5, 2, 13, 10, 2, 3, 13, 12, 13, 3, 4, 13, 12, 4, 1, 13, 9, 13, 1, 9, 8, 13, 2, 3, 12, 2, 12,+   5, 2, 5, 10, 8, 5, 12, 4, 1, 9, 2, 3, 1, 4, 1, 3, 4, 3, 12, 4, 12, 8, 4, 8, 9, 9, 8, 5, 9, 5, 1, 10, 1, 5, 2, 1, 10};+unsigned edges2Row26[] = {0};+unsigned edges2Row27[] = {4, 2, 12, 2, 11, 12, 8, 5, 9, 4, 2, 9, 2, 5, 9, 2, 11, 5, 8, 5, 11, 8, 11, 12, 4, 12, 8, 4, 8, 9, 9, 8, 5, 4, 2, 9, 2, 5, 9, 2, 11, 5, 8, 5, 11, 8, 11, 12};+unsigned edges2Row28[] = {0};+unsigned edges2Row29[] = {8, 5, 12, 5, 10, 12, 4, 1, 9, 12, 10, 11, 8, 4, 9, 8, 12, 4, 4, 11, 12, 4, 1, 11, 5, 10, 1, 5, 1, 9, 8, 5, 9, 1, 10, 11, 8, 5, 13, 5, 9, 13, 4, 13, 9, 1, 10, 13, 4, 1, 13, 11, 13, 10, 12, 13, 11, 8, 12, 13, 1, 9, 13, 8, 13, 9, 8, 5, 13, 5, 10, 13, 11, 13, 10, 12, 11, 13, 4, 13, 12, 4, 1, 13, 8, 5, 9, 4, 1, 10, 4, 10, 12, 12, 10, 11};+unsigned edges2Row30[] = {0};+unsigned edges2Row31[] = {0};+unsigned edges2Row32[] = {5, 6, 10, 1, 4, 9, 5, 6, 9, 6, 4, 9, 6, 1, 4, 6, 10, 1};+unsigned edges2Row33[] = {0};+unsigned edges2Row34[] = {0};+unsigned edges2Row35[] = {10, 5, 6, 11, 3, 2, 10, 5, 2, 5, 3, 2, 5, 11, 3, 5, 6, 11};+unsigned edges2Row36[] = {6, 10, 5, 2, 11, 3, 1, 4, 9, 6, 10, 5, 2, 9, 1, 2, 11, 9, 4, 9, 11, 4, 11, 3, 2, 11, 3, 6, 10, 1, 6, 1, 4, 6, 9, 5, 6, 4, 9, 1, 4, 9, 6, 11, 3, 6, 3, 5, 5, 3, 2, 5, 2, 10, 6, 13, 5, 5, 13, 9, 4, 9, 13, 4, 13, 3, 3, 13, 11, 2, 11, 13, 2, 13, 1, 1, 13, 10, 6, 10, 13, 6, 13, 5, 6, 13, 11, 3, 13, 11, 4, 13, 3, 4, 9, 13, 1, 9, 13, 2, 13, 1, 2, 10, 13, 5, 13, 10, 4, 9, 13, 5, 13, 9, 6, 13, 5, 6, 11, 13, 3, 13, 11, 2, 13, 3, 2, 10, 13, 1, 13, 10, 1, 4, 13, 6, 11, 3, 6, 3, 9,+   6, 9, 5, 4, 9, 3, 2, 10, 1, 6, 11, 10, 2, 10, 11, 2, 11, 3, 2, 3, 4, 2, 4, 1, 1, 4, 9, 1, 9, 10, 5, 10, 9, 6, 10, 5};+unsigned edges2Row37[] = {0};+unsigned edges2Row38[] = {0};+unsigned edges2Row39[] = {6, 10, 5, 12, 4, 3, 6, 12, 5, 12, 4, 5, 10, 5, 4, 6, 10, 3, 6, 3, 12, 10, 4, 3};+unsigned edges2Row40[] = {1, 3, 9, 3, 12, 9, 5, 6, 10, 1, 3, 10, 3, 6, 10, 3, 12, 6, 5, 6, 12, 5, 12, 9, 1, 9, 5, 1, 5, 10, 10, 5, 6, 1, 3, 10, 3, 6, 10, 3, 12, 6, 5, 6, 12, 5, 12, 9};+unsigned edges2Row41[] = {1, 5, 2, 5, 6, 2, 3, 12, 4, 1, 5, 4, 5, 12, 4, 5, 6, 12, 3, 12, 6, 3, 6, 2, 1, 2, 3, 1, 3, 4, 4, 3, 12, 1, 5, 4, 5, 12, 4, 5, 6, 12, 3, 12, 6, 3, 6, 2};+unsigned edges2Row42[] = {0};+unsigned edges2Row43[] = {11, 12, 2, 12, 4, 2, 10, 5, 6, 11, 12, 6, 12, 5, 6, 12, 4, 5, 10, 5, 4, 10, 4, 2, 11, 2, 10, 11, 10, 6, 6, 10, 5, 11, 12, 6, 12, 5, 6, 12, 4, 5, 10, 5, 4, 10, 4, 2};+unsigned edges2Row44[] = {2, 11, 1, 1, 11, 9, 6, 10, 5, 11, 12, 9, 2, 10, 6, 2, 6, 11, 6, 11, 12, 6, 12, 5, 1, 5, 9, 1, 10, 5, 2, 10, 1, 5, 12, 9, 2, 13, 1, 1, 13, 10, 6, 10, 13, 5, 13, 9, 6, 13, 5, 12, 9, 13, 11, 12, 13, 2, 13, 11, 5, 13, 10, 2, 10, 13, 2, 13, 1, 1, 13, 9, 12, 9, 13, 11, 13, 12, 6, 11, 13, 6, 13, 5, 2, 10, 1, 6, 9, 5, 6, 11, 9, 11, 12, 9};+unsigned edges2Row45[] = {0};+unsigned edges2Row46[] = {0};+unsigned edges2Row47[] = {0};+unsigned edges2Row48[] = {0};+unsigned edges2Row49[] = {0};+unsigned edges2Row50[] = {0};+unsigned edges2Row51[] = {10, 9, 6, 9, 8, 6, 11, 3, 2, 10, 9, 2, 9, 3, 2, 9, 8, 3, 11, 3, 8, 11, 8, 6, 10, 6, 11, 10, 11, 2, 2, 11, 3, 10, 9, 2, 9, 3, 2, 9, 8, 3, 11, 3, 8, 11, 8, 6};+unsigned edges2Row52[] = {1, 4, 10, 10, 4, 6, 3, 2, 11, 4, 8, 6, 1, 2, 3, 1, 3, 4, 3, 4, 8, 3, 8, 11, 10, 11, 6, 10, 2, 11, 1, 2, 10, 11, 8, 6, 1, 13, 10, 10, 13, 2, 3, 2, 13, 11, 13, 6, 3, 13, 11, 8, 6, 13, 4, 8, 13, 1, 13, 4, 11, 13, 2, 1, 2, 13, 1, 13, 10, 10, 13, 6, 8, 6, 13, 4, 13, 8, 3, 4, 13, 3, 13, 11, 1, 2, 10, 3, 6, 11, 3, 4, 6, 4, 8, 6};+unsigned edges2Row53[] = {0};+unsigned edges2Row54[] = {0};+unsigned edges2Row55[] = {8, 6, 9, 6, 10, 9, 4, 3, 12, 8, 6, 12, 6, 3, 12, 6, 10, 3, 4, 3, 10, 4, 10, 9, 8, 9, 4, 8, 4, 12, 12, 4, 3, 8, 6, 12, 6, 3, 12, 6, 10, 3, 4, 3, 10, 4, 10, 9};+unsigned edges2Row56[] = {0};+unsigned edges2Row57[] = {9, 8, 1, 1, 8, 2, 12, 4, 3, 8, 6, 2, 9, 4, 12, 9, 12, 8, 12, 8, 6, 12, 6, 3, 1, 3, 2, 1, 4, 3, 9, 4, 1, 3, 6, 2, 9, 13, 1, 1, 13, 4, 12, 4, 13, 3, 13, 2, 12, 13, 3, 6, 2, 13, 8, 6, 13, 9, 13, 8, 3, 13, 4, 9, 4, 13, 9, 13, 1, 1, 13, 2, 6, 2, 13, 8, 13, 6, 12, 8, 13, 12, 13, 3, 9, 4, 1, 12, 2, 3, 12, 8, 2, 8, 6, 2};+unsigned edges2Row58[] = {0};+unsigned edges2Row59[] = {2, 11, 4, 10, 9, 6, 6, 9, 8, 11, 12, 4, 2, 11, 6, 2, 6, 10, 6, 11, 12, 6, 12, 8, 4, 12, 9, 9, 12, 8, 2, 10, 9, 2, 9, 4, 2, 10, 13, 2, 13, 4, 4, 13, 12, 11, 12, 13, 6, 11, 13, 6, 13, 8, 9, 8, 13, 10, 9, 13, 4, 13, 9, 2, 13, 4, 2, 11, 13, 11, 12, 13, 8, 13, 12, 6, 13, 8, 10, 13, 6, 10, 9, 13, 2, 10, 4, 10, 9, 4, 6, 11, 12, 6, 12, 8};+unsigned edges2Row60[] = {6, 11, 8, 8, 11, 12, 2, 10, 1, 6, 10, 8, 8, 10, 1, 8, 1, 12, 2, 12, 1, 2, 11, 12, 6, 2, 11, 6, 10, 2, 10, 1, 2, 6, 10, 8, 8, 10, 1, 8, 1, 12, 2, 12, 1, 2, 11, 12};+unsigned edges2Row61[] = {12, 8, 11, 11, 8, 6, 9, 4, 1, 12, 4, 11, 11, 4, 1, 11, 1, 6, 9, 6, 1, 9, 8, 6, 12, 9, 8, 12, 4, 9, 4, 1, 9, 12, 4, 11, 11, 4, 1, 11, 1, 6, 9, 6, 1, 9, 8, 6};+unsigned edges2Row62[] = {0};+unsigned edges2Row63[] = {0};+unsigned edges2Row64[] = {6, 7, 11, 1, 4, 9, 6, 1, 11, 1, 4, 11, 7, 11, 4, 6, 7, 9, 6, 9, 1, 7, 4, 9};+unsigned edges2Row65[] = {6, 7, 11, 2, 1, 10, 6, 7, 10, 7, 1, 10, 7, 2, 1, 7, 11, 2};+unsigned edges2Row66[] = {2, 4, 10, 4, 9, 10, 6, 7, 11, 2, 4, 11, 4, 7, 11, 4, 9, 7, 6, 7, 9, 6, 9, 10, 2, 10, 6, 2, 6, 11, 11, 6, 7, 2, 4, 11, 4, 7, 11, 4, 9, 7, 6, 7, 9, 6, 9, 10};+unsigned edges2Row67[] = {0};+unsigned edges2Row68[] = {2, 6, 3, 6, 7, 3, 4, 9, 1, 2, 6, 1, 6, 9, 1, 6, 7, 9, 4, 9, 7, 4, 7, 3, 2, 3, 4, 2, 4, 1, 1, 4, 9, 2, 6, 1, 6, 9, 1, 6, 7, 9, 4, 9, 7, 4, 7, 3};+unsigned edges2Row69[] = {0};+unsigned edges2Row70[] = {0};+unsigned edges2Row71[] = {12, 4, 3, 11, 6, 7, 12, 4, 7, 4, 6, 7, 4, 11, 6, 4, 3, 11};+unsigned edges2Row72[] = {12, 9, 3, 9, 1, 3, 11, 6, 7, 12, 9, 7, 9, 6, 7, 9, 1, 6, 11, 6, 1, 11, 1, 3, 12, 3, 11, 12, 11, 7, 7, 11, 6, 12, 9, 7, 9, 6, 7, 9, 1, 6, 11, 6, 1, 11, 1, 3};+unsigned edges2Row73[] = {4, 3, 12, 2, 1, 10, 11, 6, 7, 4, 3, 12, 2, 7, 11, 2, 1, 7, 6, 7, 1, 6, 1, 10, 2, 1, 10, 4, 3, 11, 4, 11, 6, 4, 7, 12, 4, 6, 7, 11, 6, 7, 4, 1, 10, 4, 10, 12, 12, 10, 2, 12, 2, 3, 4, 13, 12, 12, 13, 7, 6, 7, 13, 6, 13, 10, 10, 13, 1, 2, 1, 13, 2, 13, 11, 11, 13, 3, 4, 3, 13, 4, 13, 12, 4, 13, 1, 10, 13, 1, 6, 13, 10, 6, 7, 13, 11, 7, 13, 2, 13, 11, 2, 3, 13, 12, 13, 3, 6, 7, 13, 12, 13, 7, 4, 13, 12, 4, 1, 13, 10, 13, 1, 2, 13, 10, 2, 3, 13, 11, 13, 3, 11, 6, 13, 4, 1,+   10, 4, 10, 7, 4, 7, 12, 6, 7, 10, 2, 3, 11, 4, 1, 3, 2, 3, 1, 2, 1, 10, 2, 10, 6, 2, 6, 11, 11, 6, 7, 11, 7, 3, 12, 3, 7, 4, 3, 12};+unsigned edges2Row74[] = {6, 7, 10, 7, 12, 10, 2, 3, 11, 10, 12, 9, 6, 2, 11, 6, 10, 2, 2, 9, 10, 2, 3, 9, 7, 12, 3, 7, 3, 11, 6, 7, 11, 3, 12, 9, 6, 7, 13, 7, 11, 13, 2, 13, 11, 3, 12, 13, 2, 3, 13, 9, 13, 12, 10, 13, 9, 6, 10, 13, 3, 11, 13, 6, 13, 11, 6, 7, 13, 7, 12, 13, 9, 13, 12, 10, 9, 13, 2, 13, 10, 2, 3, 13, 6, 7, 11, 2, 3, 12, 2, 12, 10, 10, 12, 9};+unsigned edges2Row75[] = {0};+unsigned edges2Row76[] = {0};+unsigned edges2Row77[] = {0};+unsigned edges2Row78[] = {0};+unsigned edges2Row79[] = {7, 11, 6, 5, 9, 8, 7, 11, 8, 11, 9, 8, 11, 5, 9, 11, 6, 5};+unsigned edges2Row80[] = {5, 1, 8, 1, 4, 8, 7, 11, 6, 5, 1, 6, 1, 11, 6, 1, 4, 11, 7, 11, 4, 7, 4, 8, 5, 8, 7, 5, 7, 6, 6, 7, 11, 5, 1, 6, 1, 11, 6, 1, 4, 11, 7, 11, 4, 7, 4, 8};+unsigned edges2Row81[] = {8, 5, 9, 6, 7, 11, 10, 2, 1, 8, 5, 9, 6, 1, 10, 6, 7, 1, 2, 1, 7, 2, 7, 11, 6, 7, 11, 8, 5, 10, 8, 10, 2, 8, 1, 9, 8, 2, 1, 10, 2, 1, 8, 7, 11, 8, 11, 9, 9, 11, 6, 9, 6, 5, 8, 13, 9, 9, 13, 1, 2, 1, 13, 2, 13, 11, 11, 13, 7, 6, 7, 13, 6, 13, 10, 10, 13, 5, 8, 5, 13, 8, 13, 9, 8, 13, 7, 11, 13, 7, 2, 13, 11, 2, 1, 13, 10, 1, 13, 6, 13, 10, 6, 5, 13, 9, 13, 5, 2, 1, 13, 9, 13, 1, 8, 13, 9, 8, 7, 13, 11, 13, 7, 6, 13, 11, 6, 5, 13, 10, 13, 5, 10, 2, 13, 8, 7, 11, 8, 11,+   1, 8, 1, 9, 2, 1, 11, 6, 5, 10, 8, 7, 5, 6, 5, 7, 6, 7, 11, 6, 11, 2, 6, 2, 10, 10, 2, 1, 10, 1, 5, 9, 5, 1, 8, 5, 9};+unsigned edges2Row82[] = {10, 2, 5, 5, 2, 8, 11, 6, 7, 2, 4, 8, 10, 6, 11, 10, 11, 2, 11, 2, 4, 11, 4, 7, 5, 7, 8, 5, 6, 7, 10, 6, 5, 7, 4, 8, 10, 13, 5, 5, 13, 6, 11, 6, 13, 7, 13, 8, 11, 13, 7, 4, 8, 13, 2, 4, 13, 10, 13, 2, 7, 13, 6, 10, 6, 13, 10, 13, 5, 5, 13, 8, 4, 8, 13, 2, 13, 4, 11, 2, 13, 11, 13, 7, 10, 6, 5, 11, 8, 7, 11, 2, 8, 2, 4, 8};+unsigned edges2Row83[] = {7, 3, 6, 3, 2, 6, 5, 9, 8, 7, 3, 8, 3, 9, 8, 3, 2, 9, 5, 9, 2, 5, 2, 6, 7, 6, 5, 7, 5, 8, 8, 5, 9, 7, 3, 8, 3, 9, 8, 3, 2, 9, 5, 9, 2, 5, 2, 6};+unsigned edges2Row84[] = {6, 7, 2, 5, 1, 8, 8, 1, 4, 7, 3, 2, 6, 7, 8, 6, 8, 5, 8, 7, 3, 8, 3, 4, 2, 3, 1, 1, 3, 4, 6, 5, 1, 6, 1, 2, 6, 5, 13, 6, 13, 2, 2, 13, 3, 7, 3, 13, 8, 7, 13, 8, 13, 4, 1, 4, 13, 5, 1, 13, 2, 13, 1, 6, 13, 2, 6, 7, 13, 7, 3, 13, 4, 13, 3, 8, 13, 4, 5, 13, 8, 5, 1, 13, 6, 5, 2, 5, 1, 2, 8, 7, 3, 8, 3, 4};+unsigned edges2Row85[] = {6, 7, 10, 10, 7, 1, 8, 5, 9, 7, 3, 1, 6, 5, 8, 6, 8, 7, 8, 7, 3, 8, 3, 9, 10, 9, 1, 10, 5, 9, 6, 5, 10, 9, 3, 1, 6, 13, 10, 10, 13, 5, 8, 5, 13, 9, 13, 1, 8, 13, 9, 3, 1, 13, 7, 3, 13, 6, 13, 7, 9, 13, 5, 6, 5, 13, 6, 13, 10, 10, 13, 1, 3, 1, 13, 7, 13, 3, 8, 7, 13, 8, 13, 9, 6, 5, 10, 8, 1, 9, 8, 7, 1, 7, 3, 1};+unsigned edges2Row86[] = {8, 7, 4, 4, 7, 3, 6, 5, 10, 8, 5, 4, 4, 5, 10, 4, 10, 3, 6, 3, 10, 6, 7, 3, 8, 6, 7, 8, 5, 6, 5, 10, 6, 8, 5, 4, 4, 5, 10, 4, 10, 3, 6, 3, 10, 6, 7, 3};+unsigned edges2Row87[] = {6, 7, 11, 8, 5, 9, 12, 4, 3, 6, 7, 11, 8, 3, 12, 8, 5, 3, 4, 3, 5, 4, 5, 9, 8, 5, 9, 6, 7, 12, 6, 12, 4, 6, 3, 11, 6, 4, 3, 12, 4, 3, 6, 5, 9, 6, 9, 11, 11, 9, 8, 11, 8, 7, 6, 13, 11, 11, 13, 3, 4, 3, 13, 4, 13, 9, 9, 13, 5, 8, 5, 13, 8, 13, 12, 12, 13, 7, 6, 7, 13, 6, 13, 11, 6, 13, 5, 9, 13, 5, 4, 13, 9, 4, 3, 13, 12, 3, 13, 8, 13, 12, 8, 7, 13, 11, 13, 7, 4, 3, 13, 11, 13, 3, 6, 13, 11, 6, 5, 13, 9, 13, 5, 8, 13, 9, 8, 7, 13, 12, 13, 7, 12, 4, 13, 6, 5, 9, 6, 9, 3,+   6, 3, 11, 4, 3, 9, 8, 7, 12, 6, 5, 7, 8, 7, 5, 8, 5, 9, 8, 9, 4, 8, 4, 12, 12, 4, 3, 12, 3, 7, 11, 7, 3, 6, 7, 11};+unsigned edges2Row88[] = {8, 5, 12, 12, 5, 3, 6, 7, 11, 5, 1, 3, 8, 7, 6, 8, 6, 5, 6, 5, 1, 6, 1, 11, 12, 11, 3, 12, 7, 11, 8, 7, 12, 11, 1, 3, 8, 13, 12, 12, 13, 7, 6, 7, 13, 11, 13, 3, 6, 13, 11, 1, 3, 13, 5, 1, 13, 8, 13, 5, 11, 13, 7, 8, 7, 13, 8, 13, 12, 12, 13, 3, 1, 3, 13, 5, 13, 1, 6, 5, 13, 6, 13, 11, 8, 7, 12, 6, 3, 11, 6, 5, 3, 5, 1, 3};+unsigned edges2Row89[] = {8, 7, 12, 10, 5, 6, 2, 3, 11, 4, 1, 9, 10, 5, 6, 8, 7, 12, 3, 11, 4, 4, 9, 11, 11, 2, 9, 9, 1, 2, 2, 3, 11, 8, 7, 12, 1, 4, 10, 10, 6, 4, 4, 9, 6, 6, 5, 9, 4, 1, 9, 2, 3, 11, 5, 10, 8, 8, 12, 10, 10, 6, 12, 12, 7, 6, 4, 1, 9, 10, 5, 6, 7, 8, 11, 11, 2, 8, 8, 12, 2, 2, 3, 12, 4, 1, 9, 8, 7, 12, 10, 2, 5, 5, 3, 2, 5, 6, 3, 3, 11, 6, 10, 5, 6, 11, 3, 2, 4, 12, 1, 1, 7, 12, 1, 9, 7, 7, 8, 9, 8, 7, 12, 13, 5, 6, 13, 9, 5, 13, 4, 9, 13, 3, 4, 13, 11, 3, 13, 2, 11, 13, 1, 2,+   13, 10, 1, 13, 6, 10, 10, 5, 6, 13, 1, 9, 13, 2, 1, 13, 11, 2, 13, 7, 11, 13, 8, 7, 13, 12, 8, 13, 3, 12, 13, 4, 3, 13, 9, 4, 8, 7, 12, 13, 4, 9, 13, 3, 4, 13, 11, 3, 13, 6, 11, 13, 5, 6, 13, 10, 5, 13, 2, 10, 13, 1, 2, 13, 9, 1, 10, 5, 6, 13, 2, 11, 13, 1, 2, 13, 9, 1, 13, 8, 9, 13, 7, 8, 13, 12, 7, 13, 4, 12, 13, 3, 4, 13, 11, 3, 2, 3, 11, 13, 7, 12, 13, 6, 7, 13, 10, 6, 13, 1, 10, 13, 4, 1, 13, 9, 4, 13, 5, 9, 13, 8, 5, 13, 12, 8, 8, 7, 12, 13, 9, 4, 13, 5, 9, 13, 6, 5, 13, 11, 6, 13, 3,+   11, 13, 2, 3, 13, 10, 2, 13, 1, 10, 13, 4, 1, 2, 3, 11, 13, 10, 6, 13, 1, 10, 13, 4, 1, 13, 12, 4, 13, 7, 12, 13, 8, 7, 13, 9, 8, 13, 5, 9, 13, 6, 5, 4, 1, 9, 13, 3, 2, 13, 12, 3, 13, 8, 12, 13, 5, 8, 13, 10, 5, 13, 6, 10, 13, 7, 6, 13, 11, 7, 13, 2, 11, 4, 1, 9, 13, 8, 12, 13, 5, 8, 13, 10, 5, 13, 2, 10, 13, 3, 2, 13, 11, 3, 13, 6, 11, 13, 7, 6, 13, 12, 7, 2, 3, 11, 13, 6, 10, 13, 7, 6, 13, 12, 7, 13, 4, 12, 13, 1, 4, 13, 9, 1, 13, 8, 9, 13, 5, 8, 13, 10, 5, 4, 1, 9, 13, 12, 8, 13, 3, 12,+   13, 2, 3, 13, 10, 2, 13, 5, 10, 13, 6, 5, 13, 11, 6, 13, 7, 11, 13, 8, 7, 10, 5, 6, 13, 11, 2, 13, 7, 11, 13, 8, 7, 13, 9, 8, 13, 1, 9, 13, 4, 1, 13, 12, 4, 13, 3, 12, 13, 2, 3, 13, 5, 6, 13, 9, 5, 13, 8, 9, 13, 7, 8, 13, 12, 7, 13, 4, 12, 13, 3, 4, 13, 11, 3, 13, 2, 11, 13, 1, 2, 13, 10, 1, 13, 6, 10, 13, 7, 11, 13, 8, 7, 13, 12, 8, 13, 3, 12, 13, 4, 3, 13, 9, 4, 13, 1, 9, 13, 2, 1, 13, 10, 2, 13, 5, 10, 13, 6, 5, 13, 11, 6, 13, 1, 10, 13, 4, 1, 13, 9, 4, 13, 5, 9, 13, 8, 5, 13, 12, 8,+   13, 7, 12, 13, 6, 7, 13, 11, 6, 13, 3, 11, 13, 2, 3, 13, 10, 2, 13, 3, 2, 13, 12, 3, 13, 4, 12, 13, 1, 4, 13, 9, 1, 13, 8, 9, 13, 5, 8, 13, 10, 5, 13, 6, 10, 13, 7, 6, 13, 11, 7, 13, 2, 11, 8, 7, 12, 10, 1, 2, 6, 5, 11, 11, 9, 5, 11, 3, 9, 9, 4, 3, 6, 5, 10, 12, 3, 4, 11, 7, 2, 2, 8, 7, 2, 1, 8, 8, 9, 1, 2, 3, 11, 9, 5, 8, 10, 1, 6, 6, 4, 1, 6, 7, 4, 4, 12, 7, 4, 1, 9, 11, 7, 6, 2, 3, 10, 10, 12, 3, 10, 5, 12, 12, 8, 5, 2, 1, 10, 12, 3, 4, 12, 11, 3, 12, 7, 11, 11, 6, 7, 7, 5, 6, 7, 8, 5,+   5, 9, 8, 8, 12, 9, 9, 4, 12, 12, 3, 4, 6, 7, 11, 6, 8, 7, 6, 5, 8, 8, 9, 5, 5, 1, 9, 5, 10, 1, 1, 2, 10, 10, 6, 2, 2, 11, 6, 9, 5, 8, 2, 1, 10, 2, 4, 1, 2, 3, 4, 4, 12, 3, 3, 7, 12, 3, 11, 7, 7, 6, 11, 11, 2, 6, 6, 10, 2, 6, 7, 11, 9, 5, 8, 9, 10, 5, 9, 1, 10, 10, 2, 1, 1, 3, 2, 1, 4, 3, 4, 12, 3, 4, 9, 12, 12, 8, 9};+unsigned edges2Row90[] = {12, 8, 7, 11, 2, 3, 6, 5, 10, 12, 8, 7, 11, 6, 5, 11, 5, 3, 10, 3, 5, 10, 2, 3, 11, 2, 3, 12, 6, 7, 12, 10, 6, 12, 8, 5, 12, 5, 10, 6, 5, 10, 12, 2, 3, 12, 8, 2, 8, 11, 2, 8, 7, 11, 12, 8, 13, 8, 5, 13, 10, 13, 5, 10, 2, 13, 2, 3, 13, 11, 13, 3, 11, 6, 13, 6, 7, 13, 12, 13, 7, 12, 8, 13, 12, 3, 13, 2, 3, 13, 10, 2, 13, 10, 13, 5, 6, 13, 5, 11, 6, 13, 11, 13, 7, 8, 7, 13, 10, 13, 5, 8, 5, 13, 12, 8, 13, 12, 13, 3, 2, 3, 13, 11, 2, 13, 11, 13, 7, 6, 7, 13, 6, 13, 10, 12,+   2, 3, 12, 5, 2, 12, 8, 5, 10, 2, 5, 11, 6, 7, 12, 7, 3, 11, 3, 7, 11, 2, 3, 11, 10, 2, 11, 6, 10, 6, 5, 10, 6, 7, 5, 8, 5, 7, 12, 8, 7};+unsigned edges2Row91[] = {12, 4, 7, 7, 4, 6, 9, 8, 5, 4, 2, 6, 12, 8, 9, 12, 9, 4, 9, 4, 2, 9, 2, 5, 7, 5, 6, 7, 8, 5, 12, 8, 7, 5, 2, 6, 12, 13, 7, 7, 13, 8, 9, 8, 13, 5, 13, 6, 9, 13, 5, 2, 6, 13, 4, 2, 13, 12, 13, 4, 5, 13, 8, 12, 8, 13, 12, 13, 7, 7, 13, 6, 2, 6, 13, 4, 13, 2, 9, 4, 13, 9, 13, 5, 12, 8, 7, 9, 6, 5, 9, 4, 6, 4, 2, 6};+unsigned edges2Row92[] = {6, 5, 2, 2, 5, 1, 8, 7, 12, 6, 7, 2, 2, 7, 12, 2, 12, 1, 8, 1, 12, 8, 5, 1, 6, 8, 5, 6, 7, 8, 7, 12, 8, 6, 7, 2, 2, 7, 12, 2, 12, 1, 8, 1, 12, 8, 5, 1};+unsigned edges2Row93[] = {7, 12, 8, 5, 10, 6, 9, 4, 1, 7, 12, 8, 5, 9, 4, 5, 4, 6, 1, 6, 4, 1, 10, 6, 5, 10, 6, 7, 9, 8, 7, 1, 9, 7, 12, 4, 7, 4, 1, 9, 4, 1, 7, 10, 6, 7, 12, 10, 12, 5, 10, 12, 8, 5, 7, 12, 13, 12, 4, 13, 1, 13, 4, 1, 10, 13, 10, 6, 13, 5, 13, 6, 5, 9, 13, 9, 8, 13, 7, 13, 8, 7, 12, 13, 7, 6, 13, 10, 6, 13, 1, 10, 13, 1, 13, 4, 9, 13, 4, 5, 9, 13, 5, 13, 8, 12, 8, 13, 1, 13, 4, 12, 4, 13, 7, 12, 13, 7, 13, 6, 10, 6, 13, 5, 10, 13, 5, 13, 8, 9, 8, 13, 9, 13, 1, 7, 10, 6, 7, 4,+   10, 7, 12, 4, 1, 10, 4, 5, 9, 8, 7, 8, 6, 5, 6, 8, 5, 10, 6, 5, 1, 10, 5, 9, 1, 9, 4, 1, 9, 8, 4, 12, 4, 8, 7, 12, 8};+unsigned edges2Row94[] = {8, 7, 12, 6, 5, 10, 8, 5, 12, 12, 5, 10, 12, 10, 6, 12, 6, 7};+unsigned edges2Row95[] = {0};+unsigned edges2Row96[] = {5, 7, 10, 7, 11, 10, 1, 4, 9, 5, 7, 9, 7, 4, 9, 7, 11, 4, 1, 4, 11, 1, 11, 10, 5, 10, 1, 5, 1, 9, 9, 1, 4, 5, 7, 9, 7, 4, 9, 7, 11, 4, 1, 4, 11, 1, 11, 10};+unsigned edges2Row97[] = {0};+unsigned edges2Row98[] = {0};+unsigned edges2Row99[] = {0};+unsigned edges2Row100[] = {4, 9, 3, 9, 5, 3, 2, 10, 1, 3, 5, 7, 4, 2, 1, 4, 3, 2, 2, 7, 3, 2, 10, 7, 9, 5, 10, 9, 10, 1, 4, 9, 1, 10, 5, 7, 4, 9, 13, 9, 1, 13, 2, 13, 1, 10, 5, 13, 2, 10, 13, 7, 13, 5, 3, 13, 7, 4, 3, 13, 10, 1, 13, 4, 13, 1, 4, 9, 13, 9, 5, 13, 7, 13, 5, 3, 7, 13, 2, 13, 3, 2, 10, 13, 4, 9, 1, 2, 10, 5, 2, 5, 3, 3, 5, 7};+unsigned edges2Row101[] = {0};+unsigned edges2Row102[] = {0};+unsigned edges2Row103[] = {11, 10, 7, 10, 5, 7, 12, 4, 3, 11, 10, 3, 10, 4, 3, 10, 5, 4, 12, 4, 5, 12, 5, 7, 11, 7, 12, 11, 12, 3, 3, 12, 4, 11, 10, 3, 10, 4, 3, 10, 5, 4, 12, 4, 5, 12, 5, 7};+unsigned edges2Row104[] = {9, 1, 12, 5, 7, 10, 10, 7, 11, 1, 3, 12, 9, 1, 10, 9, 10, 5, 10, 1, 3, 10, 3, 11, 12, 3, 7, 7, 3, 11, 9, 5, 7, 9, 7, 12, 9, 5, 13, 9, 13, 12, 12, 13, 3, 1, 3, 13, 10, 1, 13, 10, 13, 11, 7, 11, 13, 5, 7, 13, 12, 13, 7, 9, 13, 12, 9, 1, 13, 1, 3, 13, 11, 13, 3, 10, 13, 11, 5, 13, 10, 5, 7, 13, 9, 5, 12, 5, 7, 12, 10, 1, 3, 10, 3, 11};+unsigned edges2Row105[] = {2, 1, 11, 11, 1, 7, 4, 3, 12, 1, 5, 7, 2, 3, 4, 2, 4, 1, 4, 1, 5, 4, 5, 12, 11, 12, 7, 11, 3, 12, 2, 3, 11, 12, 5, 7, 2, 13, 11, 11, 13, 3, 4, 3, 13, 12, 13, 7, 4, 13, 12, 5, 7, 13, 1, 5, 13, 2, 13, 1, 12, 13, 3, 2, 3, 13, 2, 13, 11, 11, 13, 7, 5, 7, 13, 1, 13, 5, 4, 1, 13, 4, 13, 12, 2, 3, 11, 4, 7, 12, 4, 1, 7, 1, 5, 7};+unsigned edges2Row106[] = {7, 12, 5, 5, 12, 9, 3, 11, 2, 7, 11, 5, 5, 11, 2, 5, 2, 9, 3, 9, 2, 3, 12, 9, 7, 3, 12, 7, 11, 3, 11, 2, 3, 7, 11, 5, 5, 11, 2, 5, 2, 9, 3, 9, 2, 3, 12, 9};+unsigned edges2Row107[] = {0};+unsigned edges2Row108[] = {9, 5, 12, 12, 5, 7, 10, 1, 2, 9, 1, 12, 12, 1, 2, 12, 2, 7, 10, 7, 2, 10, 5, 7, 9, 10, 5, 9, 1, 10, 1, 2, 10, 9, 1, 12, 12, 1, 2, 12, 2, 7, 10, 7, 2, 10, 5, 7};+unsigned edges2Row109[] = {0};+unsigned edges2Row110[] = {0};+unsigned edges2Row111[] = {0};+unsigned edges2Row112[] = {0};+unsigned edges2Row113[] = {0};+unsigned edges2Row114[] = {0};+unsigned edges2Row115[] = {0};+unsigned edges2Row116[] = {3, 4, 7, 7, 4, 8, 1, 2, 10, 3, 2, 7, 7, 2, 10, 7, 10, 8, 1, 8, 10, 1, 4, 8, 3, 1, 4, 3, 2, 1, 2, 10, 1, 3, 2, 7, 7, 2, 10, 7, 10, 8, 1, 8, 10, 1, 4, 8};+unsigned edges2Row117[] = {0};+unsigned edges2Row118[] = {0};+unsigned edges2Row119[] = {4, 3, 9, 3, 11, 9, 8, 7, 12, 9, 11, 10, 4, 8, 12, 4, 9, 8, 8, 10, 9, 8, 7, 10, 3, 11, 7, 3, 7, 12, 4, 3, 12, 7, 11, 10, 4, 3, 13, 3, 12, 13, 8, 13, 12, 7, 11, 13, 8, 7, 13, 10, 13, 11, 9, 13, 10, 4, 9, 13, 7, 12, 13, 4, 13, 12, 4, 3, 13, 3, 11, 13, 10, 13, 11, 9, 10, 13, 8, 13, 9, 8, 7, 13, 4, 3, 12, 8, 7, 11, 8, 11, 9, 9, 11, 10};+unsigned edges2Row120[] = {3, 11, 1, 1, 11, 10, 7, 12, 8, 3, 12, 1, 1, 12, 8, 1, 8, 10, 7, 10, 8, 7, 11, 10, 3, 7, 11, 3, 12, 7, 12, 8, 7, 3, 12, 1, 1, 12, 8, 1, 8, 10, 7, 10, 8, 7, 11, 10};+unsigned edges2Row121[] = {8, 7, 12, 4, 1, 9, 3, 11, 2, 8, 7, 12, 4, 3, 11, 4, 11, 9, 2, 9, 11, 2, 1, 9, 4, 1, 9, 8, 3, 12, 8, 2, 3, 8, 7, 11, 8, 11, 2, 3, 11, 2, 8, 1, 9, 8, 7, 1, 7, 4, 1, 7, 12, 4, 8, 7, 13, 7, 11, 13, 2, 13, 11, 2, 1, 13, 1, 9, 13, 4, 13, 9, 4, 3, 13, 3, 12, 13, 8, 13, 12, 8, 7, 13, 8, 9, 13, 1, 9, 13, 2, 1, 13, 2, 13, 11, 3, 13, 11, 4, 3, 13, 4, 13, 12, 7, 12, 13, 2, 13, 11, 7, 11, 13, 8, 7, 13, 8, 13, 9, 1, 9, 13, 4, 1, 13, 4, 13, 12, 3, 12, 13, 3, 13, 2, 8, 1, 9, 8, 11,+   1, 8, 7, 11, 2, 1, 11, 4, 3, 12, 8, 12, 9, 4, 9, 12, 4, 1, 9, 4, 2, 1, 4, 3, 2, 3, 11, 2, 3, 12, 11, 7, 11, 12, 8, 7, 12};+unsigned edges2Row122[] = {3, 11, 2, 7, 12, 8, 3, 12, 2, 2, 12, 8, 2, 8, 7, 2, 7, 11};+unsigned edges2Row123[] = {9, 4, 10, 10, 4, 2, 12, 8, 7, 9, 8, 10, 10, 8, 7, 10, 7, 2, 12, 2, 7, 12, 4, 2, 9, 12, 4, 9, 8, 12, 8, 7, 12, 9, 8, 10, 10, 8, 7, 10, 7, 2, 12, 2, 7, 12, 4, 2};+unsigned edges2Row124[] = {8, 7, 12, 10, 1, 2, 8, 7, 10, 10, 7, 2, 12, 2, 7, 8, 1, 12, 8, 10, 1, 12, 1, 2};+unsigned edges2Row125[] = {9, 4, 1, 12, 8, 7, 9, 8, 1, 1, 8, 7, 1, 7, 12, 1, 12, 4};+unsigned edges2Row126[] = {0};+unsigned edges2Row127[] = {0};+unsigned edges2Row128[] = {9, 1, 4, 12, 7, 8, 9, 1, 8, 1, 7, 8, 1, 12, 7, 1, 4, 12};+unsigned edges2Row129[] = {8, 12, 7, 10, 2, 1, 8, 10, 7, 10, 2, 7, 12, 7, 2, 8, 12, 1, 8, 1, 10, 12, 2, 1};+unsigned edges2Row130[] = {9, 10, 4, 10, 2, 4, 12, 7, 8, 9, 10, 8, 10, 7, 8, 10, 2, 7, 12, 7, 2, 12, 2, 4, 9, 4, 12, 9, 12, 8, 8, 12, 7, 9, 10, 8, 10, 7, 8, 10, 2, 7, 12, 7, 2, 12, 2, 4};+unsigned edges2Row131[] = {3, 2, 11, 7, 8, 12, 3, 2, 12, 2, 8, 12, 2, 7, 8, 2, 11, 7};+unsigned edges2Row132[] = {8, 12, 7, 4, 9, 1, 3, 2, 11, 8, 12, 7, 4, 11, 3, 4, 9, 11, 2, 11, 9, 2, 9, 1, 4, 9, 1, 8, 12, 3, 8, 3, 2, 8, 11, 7, 8, 2, 11, 3, 2, 11, 8, 9, 1, 8, 1, 7, 7, 1, 4, 7, 4, 12, 8, 13, 7, 7, 13, 11, 2, 11, 13, 2, 13, 1, 1, 13, 9, 4, 9, 13, 4, 13, 3, 3, 13, 12, 8, 12, 13, 8, 13, 7, 8, 13, 9, 1, 13, 9, 2, 13, 1, 2, 11, 13, 3, 11, 13, 4, 13, 3, 4, 12, 13, 7, 13, 12, 2, 11, 13, 7, 13, 11, 8, 13, 7, 8, 9, 13, 1, 13, 9, 4, 13, 1, 4, 12, 13, 3, 13, 12, 3, 2, 13, 8, 9, 1, 8, 1,+   11, 8, 11, 7, 2, 11, 1, 4, 12, 3, 8, 9, 12, 4, 12, 9, 4, 9, 1, 4, 1, 2, 4, 2, 3, 3, 2, 11, 3, 11, 12, 7, 12, 11, 8, 12, 7};+unsigned edges2Row133[] = {3, 1, 11, 1, 10, 11, 7, 8, 12, 3, 1, 12, 1, 8, 12, 1, 10, 8, 7, 8, 10, 7, 10, 11, 3, 11, 7, 3, 7, 12, 12, 7, 8, 3, 1, 12, 1, 8, 12, 1, 10, 8, 7, 8, 10, 7, 10, 11};+unsigned edges2Row134[] = {4, 9, 3, 3, 9, 11, 8, 12, 7, 9, 10, 11, 4, 12, 8, 4, 8, 9, 8, 9, 10, 8, 10, 7, 3, 7, 11, 3, 12, 7, 4, 12, 3, 7, 10, 11, 4, 13, 3, 3, 13, 12, 8, 12, 13, 7, 13, 11, 8, 13, 7, 10, 11, 13, 9, 10, 13, 4, 13, 9, 7, 13, 12, 4, 12, 13, 4, 13, 3, 3, 13, 11, 10, 11, 13, 9, 13, 10, 8, 9, 13, 8, 13, 7, 4, 12, 3, 8, 11, 7, 8, 9, 11, 9, 10, 11};+unsigned edges2Row135[] = {0};+unsigned edges2Row136[] = {0};+unsigned edges2Row137[] = {3, 7, 4, 7, 8, 4, 1, 10, 2, 3, 7, 2, 7, 10, 2, 7, 8, 10, 1, 10, 8, 1, 8, 4, 3, 4, 1, 3, 1, 2, 2, 1, 10, 3, 7, 2, 7, 10, 2, 7, 8, 10, 1, 10, 8, 1, 8, 4};+unsigned edges2Row138[] = {0};+unsigned edges2Row139[] = {0};+unsigned edges2Row140[] = {0};+unsigned edges2Row141[] = {0};+unsigned edges2Row142[] = {0};+unsigned edges2Row143[] = {0};+unsigned edges2Row144[] = {0};+unsigned edges2Row145[] = {9, 12, 5, 12, 7, 5, 10, 2, 1, 9, 12, 1, 12, 2, 1, 12, 7, 2, 10, 2, 7, 10, 7, 5, 9, 5, 10, 9, 10, 1, 1, 10, 2, 9, 12, 1, 12, 2, 1, 12, 7, 2, 10, 2, 7, 10, 7, 5};+unsigned edges2Row146[] = {0};+unsigned edges2Row147[] = {7, 5, 12, 5, 9, 12, 3, 2, 11, 7, 5, 11, 5, 2, 11, 5, 9, 2, 3, 2, 9, 3, 9, 12, 7, 12, 3, 7, 3, 11, 11, 3, 2, 7, 5, 11, 5, 2, 11, 5, 9, 2, 3, 2, 9, 3, 9, 12};+unsigned edges2Row148[] = {2, 11, 1, 11, 7, 1, 4, 12, 3, 1, 7, 5, 2, 4, 3, 2, 1, 4, 4, 5, 1, 4, 12, 5, 11, 7, 12, 11, 12, 3, 2, 11, 3, 12, 7, 5, 2, 11, 13, 11, 3, 13, 4, 13, 3, 12, 7, 13, 4, 12, 13, 5, 13, 7, 1, 13, 5, 2, 1, 13, 12, 3, 13, 2, 13, 3, 2, 11, 13, 11, 7, 13, 5, 13, 7, 1, 5, 13, 4, 13, 1, 4, 12, 13, 2, 11, 3, 4, 12, 7, 4, 7, 1, 1, 7, 5};+unsigned edges2Row149[] = {9, 12, 1, 5, 10, 7, 10, 11, 7, 1, 12, 3, 9, 10, 1, 9, 5, 10, 10, 3, 1, 10, 11, 3, 12, 7, 3, 7, 11, 3, 9, 7, 5, 9, 12, 7, 9, 13, 5, 9, 12, 13, 12, 3, 13, 1, 13, 3, 10, 13, 1, 10, 11, 13, 7, 13, 11, 5, 13, 7, 12, 7, 13, 9, 12, 13, 9, 13, 1, 1, 13, 3, 11, 3, 13, 10, 11, 13, 5, 10, 13, 5, 13, 7, 9, 12, 5, 5, 12, 7, 10, 3, 1, 10, 11, 3};+unsigned edges2Row150[] = {11, 7, 10, 10, 7, 5, 12, 3, 4, 11, 3, 10, 10, 3, 4, 10, 4, 5, 12, 5, 4, 12, 7, 5, 11, 12, 7, 11, 3, 12, 3, 4, 12, 11, 3, 10, 10, 3, 4, 10, 4, 5, 12, 5, 4, 12, 7, 5};+unsigned edges2Row151[] = {0};+unsigned edges2Row152[] = {0};+unsigned edges2Row153[] = {4, 3, 9, 9, 3, 5, 2, 1, 10, 3, 7, 5, 4, 1, 2, 4, 2, 3, 2, 3, 7, 2, 7, 10, 9, 10, 5, 9, 1, 10, 4, 1, 9, 10, 7, 5, 4, 13, 9, 9, 13, 1, 2, 1, 13, 10, 13, 5, 2, 13, 10, 7, 5, 13, 3, 7, 13, 4, 13, 3, 10, 13, 1, 4, 1, 13, 4, 13, 9, 9, 13, 5, 7, 5, 13, 3, 13, 7, 2, 3, 13, 2, 13, 10, 4, 1, 9, 2, 5, 10, 2, 3, 5, 3, 7, 5};+unsigned edges2Row154[] = {0};+unsigned edges2Row155[] = {0};+unsigned edges2Row156[] = {0};+unsigned edges2Row157[] = {5, 10, 7, 7, 10, 11, 1, 9, 4, 5, 9, 7, 7, 9, 4, 7, 4, 11, 1, 11, 4, 1, 10, 11, 5, 1, 10, 5, 9, 1, 9, 4, 1, 5, 9, 7, 7, 9, 4, 7, 4, 11, 1, 11, 4, 1, 10, 11};+unsigned edges2Row158[] = {0};+unsigned edges2Row159[] = {8, 12, 7, 6, 10, 5, 8, 12, 5, 12, 10, 5, 12, 6, 10, 12, 7, 6};+unsigned edges2Row160[] = {7, 8, 12, 5, 6, 10, 9, 1, 4, 7, 8, 12, 5, 4, 9, 5, 6, 4, 1, 4, 6, 1, 6, 10, 5, 6, 10, 7, 8, 9, 7, 9, 1, 7, 4, 12, 7, 1, 4, 9, 1, 4, 7, 6, 10, 7, 10, 12, 12, 10, 5, 12, 5, 8, 7, 13, 12, 12, 13, 4, 1, 4, 13, 1, 13, 10, 10, 13, 6, 5, 6, 13, 5, 13, 9, 9, 13, 8, 7, 8, 13, 7, 13, 12, 7, 13, 6, 10, 13, 6, 1, 13, 10, 1, 4, 13, 9, 4, 13, 5, 13, 9, 5, 8, 13, 12, 13, 8, 1, 4, 13, 12, 13, 4, 7, 13, 12, 7, 6, 13, 10, 13, 6, 5, 13, 10, 5, 8, 13, 9, 13, 8, 9, 1, 13, 7, 6, 10, 7, 10,+   4, 7, 4, 12, 1, 4, 10, 5, 8, 9, 7, 6, 8, 5, 8, 6, 5, 6, 10, 5, 10, 1, 5, 1, 9, 9, 1, 4, 9, 4, 8, 12, 8, 4, 7, 8, 12};+unsigned edges2Row161[] = {6, 2, 5, 2, 1, 5, 8, 12, 7, 6, 2, 7, 2, 12, 7, 2, 1, 12, 8, 12, 1, 8, 1, 5, 6, 5, 8, 6, 8, 7, 7, 8, 12, 6, 2, 7, 2, 12, 7, 2, 1, 12, 8, 12, 1, 8, 1, 5};+unsigned edges2Row162[] = {12, 7, 4, 7, 6, 4, 9, 5, 8, 4, 6, 2, 12, 9, 8, 12, 4, 9, 9, 2, 4, 9, 5, 2, 7, 6, 5, 7, 5, 8, 12, 7, 8, 5, 6, 2, 12, 7, 13, 7, 8, 13, 9, 13, 8, 5, 6, 13, 9, 5, 13, 2, 13, 6, 4, 13, 2, 12, 4, 13, 5, 8, 13, 12, 13, 8, 12, 7, 13, 7, 6, 13, 2, 13, 6, 4, 2, 13, 9, 13, 4, 9, 5, 13, 12, 7, 8, 9, 5, 6, 9, 6, 4, 4, 6, 2};+unsigned edges2Row163[] = {12, 7, 8, 11, 3, 2, 6, 10, 5, 12, 7, 8, 11, 5, 6, 11, 3, 5, 10, 5, 3, 10, 3, 2, 11, 3, 2, 12, 7, 6, 12, 6, 10, 12, 5, 8, 12, 10, 5, 6, 10, 5, 12, 3, 2, 12, 2, 8, 8, 2, 11, 8, 11, 7, 12, 13, 8, 8, 13, 5, 10, 5, 13, 10, 13, 2, 2, 13, 3, 11, 3, 13, 11, 13, 6, 6, 13, 7, 12, 7, 13, 12, 13, 8, 12, 13, 3, 2, 13, 3, 10, 13, 2, 10, 5, 13, 6, 5, 13, 11, 13, 6, 11, 7, 13, 8, 13, 7, 10, 5, 13, 8, 13, 5, 12, 13, 8, 12, 3, 13, 2, 13, 3, 11, 13, 2, 11, 7, 13, 6, 13, 7, 6, 10, 13,+   12, 3, 2, 12, 2, 5, 12, 5, 8, 10, 5, 2, 11, 7, 6, 12, 3, 7, 11, 7, 3, 11, 3, 2, 11, 2, 10, 11, 10, 6, 6, 10, 5, 6, 5, 7, 8, 7, 5, 12, 7, 8};+unsigned edges2Row164[] = {8, 12, 7, 10, 6, 5, 2, 11, 3, 4, 9, 1, 10, 6, 5, 8, 12, 7, 3, 4, 11, 4, 11, 9, 11, 9, 2, 9, 2, 1, 2, 11, 3, 8, 12, 7, 1, 10, 4, 10, 4, 6, 4, 6, 9, 6, 9, 5, 4, 9, 1, 2, 11, 3, 5, 8, 10, 8, 10, 12, 10, 12, 6, 12, 6, 7, 4, 9, 1, 10, 6, 5, 7, 11, 8, 11, 8, 2, 8, 2, 12, 2, 12, 3, 4, 9, 1, 8, 12, 7, 10, 5, 2, 5, 2, 3, 5, 3, 6, 3, 6, 11, 10, 6, 5, 11, 2, 3, 4, 1, 12, 1, 12, 7, 1, 7, 9, 7, 9, 8, 8, 12, 7, 13, 6, 5, 13, 5, 9, 13, 9, 4, 13, 4, 3, 13, 3, 11, 13, 11, 2, 13, 2, 1,+   13, 1, 10, 13, 10, 6, 10, 6, 5, 13, 9, 1, 13, 1, 2, 13, 2, 11, 13, 11, 7, 13, 7, 8, 13, 8, 12, 13, 12, 3, 13, 3, 4, 13, 4, 9, 8, 12, 7, 13, 9, 4, 13, 4, 3, 13, 3, 11, 13, 11, 6, 13, 6, 5, 13, 5, 10, 13, 10, 2, 13, 2, 1, 13, 1, 9, 10, 6, 5, 13, 11, 2, 13, 2, 1, 13, 1, 9, 13, 9, 8, 13, 8, 7, 13, 7, 12, 13, 12, 4, 13, 4, 3, 13, 3, 11, 2, 11, 3, 13, 12, 7, 13, 7, 6, 13, 6, 10, 13, 10, 1, 13, 1, 4, 13, 4, 9, 13, 9, 5, 13, 5, 8, 13, 8, 12, 8, 12, 7, 13, 4, 9, 13, 9, 5, 13, 5, 6, 13, 6, 11, 13, 11,+   3, 13, 3, 2, 13, 2, 10, 13, 10, 1, 13, 1, 4, 2, 11, 3, 13, 6, 10, 13, 10, 1, 13, 1, 4, 13, 4, 12, 13, 12, 7, 13, 7, 8, 13, 8, 9, 13, 9, 5, 13, 5, 6, 4, 9, 1, 13, 2, 3, 13, 3, 12, 13, 12, 8, 13, 8, 5, 13, 5, 10, 13, 10, 6, 13, 6, 7, 13, 7, 11, 13, 11, 2, 4, 9, 1, 13, 12, 8, 13, 8, 5, 13, 5, 10, 13, 10, 2, 13, 2, 3, 13, 3, 11, 13, 11, 6, 13, 6, 7, 13, 7, 12, 2, 11, 3, 13, 10, 6, 13, 6, 7, 13, 7, 12, 13, 12, 4, 13, 4, 1, 13, 1, 9, 13, 9, 8, 13, 8, 5, 13, 5, 10, 4, 9, 1, 13, 8, 12, 13, 12, 3,+   13, 3, 2, 13, 2, 10, 13, 10, 5, 13, 5, 6, 13, 6, 11, 13, 11, 7, 13, 7, 8, 10, 6, 5, 13, 2, 11, 13, 11, 7, 13, 7, 8, 13, 8, 9, 13, 9, 1, 13, 1, 4, 13, 4, 12, 13, 12, 3, 13, 3, 2, 13, 6, 5, 13, 5, 9, 13, 9, 8, 13, 8, 7, 13, 7, 12, 13, 12, 4, 13, 4, 3, 13, 3, 11, 13, 11, 2, 13, 2, 1, 13, 1, 10, 13, 10, 6, 13, 11, 7, 13, 7, 8, 13, 8, 12, 13, 12, 3, 13, 3, 4, 13, 4, 9, 13, 9, 1, 13, 1, 2, 13, 2, 10, 13, 10, 5, 13, 5, 6, 13, 6, 11, 13, 10, 1, 13, 1, 4, 13, 4, 9, 13, 9, 5, 13, 5, 8, 13, 8, 12,+   13, 12, 7, 13, 7, 6, 13, 6, 11, 13, 11, 3, 13, 3, 2, 13, 2, 10, 13, 2, 3, 13, 3, 12, 13, 12, 4, 13, 4, 1, 13, 1, 9, 13, 9, 8, 13, 8, 5, 13, 5, 10, 13, 10, 6, 13, 6, 7, 13, 7, 11, 13, 11, 2, 8, 12, 7, 10, 2, 1, 6, 11, 5, 11, 5, 9, 11, 9, 3, 9, 3, 4, 6, 10, 5, 12, 4, 3, 11, 2, 7, 2, 7, 8, 2, 8, 1, 8, 1, 9, 2, 11, 3, 9, 8, 5, 10, 6, 1, 6, 1, 4, 6, 4, 7, 4, 7, 12, 4, 9, 1, 11, 6, 7, 2, 10, 3, 10, 3, 12, 10, 12, 5, 12, 5, 8, 2, 10, 1, 12, 4, 3, 12, 3, 11, 12, 11, 7, 11, 7, 6, 7, 6, 5, 7, 5, 8,+   5, 8, 9, 8, 9, 12, 9, 12, 4, 12, 4, 3, 6, 11, 7, 6, 7, 8, 6, 8, 5, 8, 5, 9, 5, 9, 1, 5, 1, 10, 1, 10, 2, 10, 2, 6, 2, 6, 11, 9, 8, 5, 2, 10, 1, 2, 1, 4, 2, 4, 3, 4, 3, 12, 3, 12, 7, 3, 7, 11, 7, 11, 6, 11, 6, 2, 6, 2, 10, 6, 11, 7, 9, 8, 5, 9, 5, 10, 9, 10, 1, 10, 1, 2, 1, 2, 3, 1, 3, 4, 4, 3, 12, 4, 12, 9, 12, 9, 8};+unsigned edges2Row165[] = {8, 12, 5, 12, 3, 5, 6, 11, 7, 5, 3, 1, 8, 6, 7, 8, 5, 6, 6, 1, 5, 6, 11, 1, 12, 3, 11, 12, 11, 7, 8, 12, 7, 11, 3, 1, 8, 12, 13, 12, 7, 13, 6, 13, 7, 11, 3, 13, 6, 11, 13, 1, 13, 3, 5, 13, 1, 8, 5, 13, 11, 7, 13, 8, 13, 7, 8, 12, 13, 12, 3, 13, 1, 13, 3, 5, 1, 13, 6, 13, 5, 6, 11, 13, 8, 12, 7, 6, 11, 3, 6, 3, 5, 5, 3, 1};+unsigned edges2Row166[] = {6, 11, 7, 8, 9, 5, 12, 3, 4, 6, 11, 7, 8, 12, 3, 8, 3, 5, 4, 5, 3, 4, 9, 5, 8, 9, 5, 6, 12, 7, 6, 4, 12, 6, 11, 3, 6, 3, 4, 12, 3, 4, 6, 9, 5, 6, 11, 9, 11, 8, 9, 11, 7, 8, 6, 11, 13, 11, 3, 13, 4, 13, 3, 4, 9, 13, 9, 5, 13, 8, 13, 5, 8, 12, 13, 12, 7, 13, 6, 13, 7, 6, 11, 13, 6, 5, 13, 9, 5, 13, 4, 9, 13, 4, 13, 3, 12, 13, 3, 8, 12, 13, 8, 13, 7, 11, 7, 13, 4, 13, 3, 11, 3, 13, 6, 11, 13, 6, 13, 5, 9, 5, 13, 8, 9, 13, 8, 13, 7, 12, 7, 13, 12, 13, 4, 6, 9, 5, 6, 3, 9,+   6, 11, 3, 4, 9, 3, 8, 12, 7, 6, 7, 5, 8, 5, 7, 8, 9, 5, 8, 4, 9, 8, 12, 4, 12, 3, 4, 12, 7, 3, 11, 3, 7, 6, 11, 7};+unsigned edges2Row167[] = {8, 4, 7, 4, 3, 7, 6, 10, 5, 8, 4, 5, 4, 10, 5, 4, 3, 10, 6, 10, 3, 6, 3, 7, 8, 7, 6, 8, 6, 5, 5, 6, 10, 8, 4, 5, 4, 10, 5, 4, 3, 10, 6, 10, 3, 6, 3, 7};+unsigned edges2Row168[] = {6, 10, 7, 10, 1, 7, 8, 9, 5, 7, 1, 3, 6, 8, 5, 6, 7, 8, 8, 3, 7, 8, 9, 3, 10, 1, 9, 10, 9, 5, 6, 10, 5, 9, 1, 3, 6, 10, 13, 10, 5, 13, 8, 13, 5, 9, 1, 13, 8, 9, 13, 3, 13, 1, 7, 13, 3, 6, 7, 13, 9, 5, 13, 6, 13, 5, 6, 10, 13, 10, 1, 13, 3, 13, 1, 7, 3, 13, 8, 13, 7, 8, 9, 13, 6, 10, 5, 8, 9, 1, 8, 1, 7, 7, 1, 3};+unsigned edges2Row169[] = {6, 2, 7, 5, 8, 1, 8, 4, 1, 7, 2, 3, 6, 8, 7, 6, 5, 8, 8, 3, 7, 8, 4, 3, 2, 1, 3, 1, 4, 3, 6, 1, 5, 6, 2, 1, 6, 13, 5, 6, 2, 13, 2, 3, 13, 7, 13, 3, 8, 13, 7, 8, 4, 13, 1, 13, 4, 5, 13, 1, 2, 1, 13, 6, 2, 13, 6, 13, 7, 7, 13, 3, 4, 3, 13, 8, 4, 13, 5, 8, 13, 5, 13, 1, 6, 2, 5, 5, 2, 1, 8, 3, 7, 8, 4, 3};+unsigned edges2Row170[] = {7, 6, 3, 3, 6, 2, 5, 8, 9, 7, 8, 3, 3, 8, 9, 3, 9, 2, 5, 2, 9, 5, 6, 2, 7, 5, 6, 7, 8, 5, 8, 9, 5, 7, 8, 3, 3, 8, 9, 3, 9, 2, 5, 2, 9, 5, 6, 2};+unsigned edges2Row171[] = {10, 5, 2, 5, 8, 2, 11, 7, 6, 2, 8, 4, 10, 11, 6, 10, 2, 11, 11, 4, 2, 11, 7, 4, 5, 8, 7, 5, 7, 6, 10, 5, 6, 7, 8, 4, 10, 5, 13, 5, 6, 13, 11, 13, 6, 7, 8, 13, 11, 7, 13, 4, 13, 8, 2, 13, 4, 10, 2, 13, 7, 6, 13, 10, 13, 6, 10, 5, 13, 5, 8, 13, 4, 13, 8, 2, 4, 13, 11, 13, 2, 11, 7, 13, 10, 5, 6, 11, 7, 8, 11, 8, 2, 2, 8, 4};+unsigned edges2Row172[] = {8, 9, 5, 6, 11, 7, 10, 1, 2, 8, 9, 5, 6, 10, 1, 6, 1, 7, 2, 7, 1, 2, 11, 7, 6, 11, 7, 8, 10, 5, 8, 2, 10, 8, 9, 1, 8, 1, 2, 10, 1, 2, 8, 11, 7, 8, 9, 11, 9, 6, 11, 9, 5, 6, 8, 9, 13, 9, 1, 13, 2, 13, 1, 2, 11, 13, 11, 7, 13, 6, 13, 7, 6, 10, 13, 10, 5, 13, 8, 13, 5, 8, 9, 13, 8, 7, 13, 11, 7, 13, 2, 11, 13, 2, 13, 1, 10, 13, 1, 6, 10, 13, 6, 13, 5, 9, 5, 13, 2, 13, 1, 9, 1, 13, 8, 9, 13, 8, 13, 7, 11, 7, 13, 6, 11, 13, 6, 13, 5, 10, 5, 13, 10, 13, 2, 8, 11, 7, 8, 1,+   11, 8, 9, 1, 2, 11, 1, 6, 10, 5, 8, 5, 7, 6, 7, 5, 6, 11, 7, 6, 2, 11, 6, 10, 2, 10, 1, 2, 10, 5, 1, 9, 1, 5, 8, 9, 5};+unsigned edges2Row173[] = {5, 8, 1, 1, 8, 4, 7, 6, 11, 5, 6, 1, 1, 6, 11, 1, 11, 4, 7, 4, 11, 7, 8, 4, 5, 7, 8, 5, 6, 7, 6, 11, 7, 5, 6, 1, 1, 6, 11, 1, 11, 4, 7, 4, 11, 7, 8, 4};+unsigned edges2Row174[] = {7, 6, 11, 5, 8, 9, 7, 8, 11, 11, 8, 9, 11, 9, 5, 11, 5, 6};+unsigned edges2Row175[] = {0};+unsigned edges2Row176[] = {0};+unsigned edges2Row177[] = {0};+unsigned edges2Row178[] = {0};+unsigned edges2Row179[] = {6, 10, 7, 7, 10, 12, 2, 11, 3, 10, 9, 12, 6, 11, 2, 6, 2, 10, 2, 10, 9, 2, 9, 3, 7, 3, 12, 7, 11, 3, 6, 11, 7, 3, 9, 12, 6, 13, 7, 7, 13, 11, 2, 11, 13, 3, 13, 12, 2, 13, 3, 9, 12, 13, 10, 9, 13, 6, 13, 10, 3, 13, 11, 6, 11, 13, 6, 13, 7, 7, 13, 12, 9, 12, 13, 10, 13, 9, 2, 10, 13, 2, 13, 3, 6, 11, 7, 2, 12, 3, 2, 10, 12, 10, 9, 12};+unsigned edges2Row180[] = {4, 12, 3, 2, 10, 1, 11, 7, 6, 4, 12, 3, 2, 11, 7, 2, 7, 1, 6, 1, 7, 6, 10, 1, 2, 10, 1, 4, 11, 3, 4, 6, 11, 4, 12, 7, 4, 7, 6, 11, 7, 6, 4, 10, 1, 4, 12, 10, 12, 2, 10, 12, 3, 2, 4, 12, 13, 12, 7, 13, 6, 13, 7, 6, 10, 13, 10, 1, 13, 2, 13, 1, 2, 11, 13, 11, 3, 13, 4, 13, 3, 4, 12, 13, 4, 1, 13, 10, 1, 13, 6, 10, 13, 6, 13, 7, 11, 13, 7, 2, 11, 13, 2, 13, 3, 12, 3, 13, 6, 13, 7, 12, 7, 13, 4, 12, 13, 4, 13, 1, 10, 1, 13, 2, 10, 13, 2, 13, 3, 11, 3, 13, 11, 13, 6, 4,+   10, 1, 4, 7, 10, 4, 12, 7, 6, 10, 7, 2, 11, 3, 4, 3, 1, 2, 1, 3, 2, 10, 1, 2, 6, 10, 2, 11, 6, 11, 7, 6, 11, 3, 7, 12, 7, 3, 4, 12, 3};+unsigned edges2Row181[] = {12, 3, 9, 9, 3, 1, 11, 7, 6, 12, 7, 9, 9, 7, 6, 9, 6, 1, 11, 1, 6, 11, 3, 1, 12, 11, 3, 12, 7, 11, 7, 6, 11, 12, 7, 9, 9, 7, 6, 9, 6, 1, 11, 1, 6, 11, 3, 1};+unsigned edges2Row182[] = {12, 3, 4, 11, 7, 6, 12, 7, 4, 4, 7, 6, 4, 6, 11, 4, 11, 3};+unsigned edges2Row183[] = {0};+unsigned edges2Row184[] = {0};+unsigned edges2Row185[] = {2, 3, 6, 6, 3, 7, 4, 1, 9, 2, 1, 6, 6, 1, 9, 6, 9, 7, 4, 7, 9, 4, 3, 7, 2, 4, 3, 2, 1, 4, 1, 9, 4, 2, 1, 6, 6, 1, 9, 6, 9, 7, 4, 7, 9, 4, 3, 7};+unsigned edges2Row186[] = {0};+unsigned edges2Row187[] = {2, 10, 4, 4, 10, 9, 6, 11, 7, 2, 11, 4, 4, 11, 7, 4, 7, 9, 6, 9, 7, 6, 10, 9, 2, 6, 10, 2, 11, 6, 11, 7, 6, 2, 11, 4, 4, 11, 7, 4, 7, 9, 6, 9, 7, 6, 10, 9};+unsigned edges2Row188[] = {6, 11, 7, 2, 10, 1, 6, 10, 7, 7, 10, 1, 7, 1, 2, 7, 2, 11};+unsigned edges2Row189[] = {6, 11, 7, 1, 9, 4, 6, 11, 1, 1, 11, 4, 7, 4, 11, 6, 9, 7, 6, 1, 9, 7, 9, 4};+unsigned edges2Row190[] = {0};+unsigned edges2Row191[] = {0};+unsigned edges2Row192[] = {12, 11, 8, 11, 6, 8, 9, 1, 4, 12, 11, 4, 11, 1, 4, 11, 6, 1, 9, 1, 6, 9, 6, 8, 12, 8, 9, 12, 9, 4, 4, 9, 1, 12, 11, 4, 11, 1, 4, 11, 6, 1, 9, 1, 6, 9, 6, 8};+unsigned edges2Row193[] = {6, 8, 11, 8, 12, 11, 2, 1, 10, 6, 8, 10, 8, 1, 10, 8, 12, 1, 2, 1, 12, 2, 12, 11, 6, 11, 2, 6, 2, 10, 10, 2, 1, 6, 8, 10, 8, 1, 10, 8, 12, 1, 2, 1, 12, 2, 12, 11};+unsigned edges2Row194[] = {2, 4, 11, 10, 6, 9, 6, 8, 9, 11, 4, 12, 2, 6, 11, 2, 10, 6, 6, 12, 11, 6, 8, 12, 4, 9, 12, 9, 8, 12, 2, 9, 10, 2, 4, 9, 2, 13, 10, 2, 4, 13, 4, 12, 13, 11, 13, 12, 6, 13, 11, 6, 8, 13, 9, 13, 8, 10, 13, 9, 4, 9, 13, 2, 4, 13, 2, 13, 11, 11, 13, 12, 8, 12, 13, 6, 8, 13, 10, 6, 13, 10, 13, 9, 2, 4, 10, 10, 4, 9, 6, 12, 11, 6, 8, 12};+unsigned edges2Row195[] = {0};+unsigned edges2Row196[] = {9, 1, 8, 1, 2, 8, 12, 3, 4, 8, 2, 6, 9, 12, 4, 9, 8, 12, 12, 6, 8, 12, 3, 6, 1, 2, 3, 1, 3, 4, 9, 1, 4, 3, 2, 6, 9, 1, 13, 1, 4, 13, 12, 13, 4, 3, 2, 13, 12, 3, 13, 6, 13, 2, 8, 13, 6, 9, 8, 13, 3, 4, 13, 9, 13, 4, 9, 1, 13, 1, 2, 13, 6, 13, 2, 8, 6, 13, 12, 13, 8, 12, 3, 13, 9, 1, 4, 12, 3, 2, 12, 2, 8, 8, 2, 6};+unsigned edges2Row197[] = {0};+unsigned edges2Row198[] = {8, 9, 6, 6, 9, 10, 4, 12, 3, 8, 12, 6, 6, 12, 3, 6, 3, 10, 4, 10, 3, 4, 9, 10, 8, 4, 9, 8, 12, 4, 12, 3, 4, 8, 12, 6, 6, 12, 3, 6, 3, 10, 4, 10, 3, 4, 9, 10};+unsigned edges2Row199[] = {0};+unsigned edges2Row200[] = {0};+unsigned edges2Row201[] = {1, 10, 4, 10, 6, 4, 3, 11, 2, 4, 6, 8, 1, 3, 2, 1, 4, 3, 3, 8, 4, 3, 11, 8, 10, 6, 11, 10, 11, 2, 1, 10, 2, 11, 6, 8, 1, 10, 13, 10, 2, 13, 3, 13, 2, 11, 6, 13, 3, 11, 13, 8, 13, 6, 4, 13, 8, 1, 4, 13, 11, 2, 13, 1, 13, 2, 1, 10, 13, 10, 6, 13, 8, 13, 6, 4, 8, 13, 3, 13, 4, 3, 11, 13, 1, 10, 2, 3, 11, 6, 3, 6, 4, 4, 6, 8};+unsigned edges2Row202[] = {10, 6, 9, 9, 6, 8, 11, 2, 3, 10, 2, 9, 9, 2, 3, 9, 3, 8, 11, 8, 3, 11, 6, 8, 10, 11, 6, 10, 2, 11, 2, 3, 11, 10, 2, 9, 9, 2, 3, 9, 3, 8, 11, 8, 3, 11, 6, 8};+unsigned edges2Row203[] = {0};+unsigned edges2Row204[] = {0};+unsigned edges2Row205[] = {0};+unsigned edges2Row206[] = {0};+unsigned edges2Row207[] = {0};+unsigned edges2Row208[] = {0};+unsigned edges2Row209[] = {2, 1, 11, 1, 9, 11, 6, 5, 10, 11, 9, 12, 2, 6, 10, 2, 11, 6, 6, 12, 11, 6, 5, 12, 1, 9, 5, 1, 5, 10, 2, 1, 10, 5, 9, 12, 2, 1, 13, 1, 10, 13, 6, 13, 10, 5, 9, 13, 6, 5, 13, 12, 13, 9, 11, 13, 12, 2, 11, 13, 5, 10, 13, 2, 13, 10, 2, 1, 13, 1, 9, 13, 12, 13, 9, 11, 12, 13, 6, 13, 11, 6, 5, 13, 2, 1, 10, 6, 5, 9, 6, 9, 11, 11, 9, 12};+unsigned edges2Row210[] = {11, 2, 12, 12, 2, 4, 10, 6, 5, 11, 6, 12, 12, 6, 5, 12, 5, 4, 10, 4, 5, 10, 2, 4, 11, 10, 2, 11, 6, 10, 6, 5, 10, 11, 6, 12, 12, 6, 5, 12, 5, 4, 10, 4, 5, 10, 2, 4};+unsigned edges2Row211[] = {0};+unsigned edges2Row212[] = {1, 2, 5, 5, 2, 6, 3, 4, 12, 1, 4, 5, 5, 4, 12, 5, 12, 6, 3, 6, 12, 3, 2, 6, 1, 3, 2, 1, 4, 3, 4, 12, 3, 1, 4, 5, 5, 4, 12, 5, 12, 6, 3, 6, 12, 3, 2, 6};+unsigned edges2Row213[] = {1, 9, 3, 3, 9, 12, 5, 10, 6, 1, 10, 3, 3, 10, 6, 3, 6, 12, 5, 12, 6, 5, 9, 12, 1, 5, 9, 1, 10, 5, 10, 6, 5, 1, 10, 3, 3, 10, 6, 3, 6, 12, 5, 12, 6, 5, 9, 12};+unsigned edges2Row214[] = {6, 5, 10, 12, 3, 4, 6, 5, 12, 12, 5, 4, 10, 4, 5, 6, 3, 10, 6, 12, 3, 10, 3, 4};+unsigned edges2Row215[] = {0};+unsigned edges2Row216[] = {0};+unsigned edges2Row217[] = {6, 5, 10, 2, 3, 11, 1, 9, 4, 6, 5, 10, 2, 1, 9, 2, 9, 11, 4, 11, 9, 4, 3, 11, 2, 3, 11, 6, 1, 10, 6, 4, 1, 6, 5, 9, 6, 9, 4, 1, 9, 4, 6, 3, 11, 6, 5, 3, 5, 2, 3, 5, 10, 2, 6, 5, 13, 5, 9, 13, 4, 13, 9, 4, 3, 13, 3, 11, 13, 2, 13, 11, 2, 1, 13, 1, 10, 13, 6, 13, 10, 6, 5, 13, 6, 11, 13, 3, 11, 13, 4, 3, 13, 4, 13, 9, 1, 13, 9, 2, 1, 13, 2, 13, 10, 5, 10, 13, 4, 13, 9, 5, 9, 13, 6, 5, 13, 6, 13, 11, 3, 11, 13, 2, 3, 13, 2, 13, 10, 1, 10, 13, 1, 13, 4, 6, 3, 11, 6, 9, 3,+   6, 5, 9, 4, 3, 9, 2, 1, 10, 6, 10, 11, 2, 11, 10, 2, 3, 11, 2, 4, 3, 2, 1, 4, 1, 9, 4, 1, 10, 9, 5, 9, 10, 6, 5, 10};+unsigned edges2Row218[] = {10, 6, 5, 11, 2, 3, 10, 2, 5, 5, 2, 3, 5, 3, 11, 5, 11, 6};+unsigned edges2Row219[] = {0};+unsigned edges2Row220[] = {0};+unsigned edges2Row221[] = {5, 10, 6, 1, 9, 4, 5, 9, 6, 6, 9, 4, 6, 4, 1, 6, 1, 10};+unsigned edges2Row222[] = {0};+unsigned edges2Row223[] = {0};+unsigned edges2Row224[] = {8, 12, 5, 5, 12, 10, 4, 9, 1, 12, 11, 10, 8, 9, 4, 8, 4, 12, 4, 12, 11, 4, 11, 1, 5, 1, 10, 5, 9, 1, 8, 9, 5, 1, 11, 10, 8, 13, 5, 5, 13, 9, 4, 9, 13, 1, 13, 10, 4, 13, 1, 11, 10, 13, 12, 11, 13, 8, 13, 12, 1, 13, 9, 8, 9, 13, 8, 13, 5, 5, 13, 10, 11, 10, 13, 12, 13, 11, 4, 12, 13, 4, 13, 1, 8, 9, 5, 4, 10, 1, 4, 12, 10, 12, 11, 10};+unsigned edges2Row225[] = {0};+unsigned edges2Row226[] = {4, 12, 2, 2, 12, 11, 8, 9, 5, 4, 9, 2, 2, 9, 5, 2, 5, 11, 8, 11, 5, 8, 12, 11, 4, 8, 12, 4, 9, 8, 9, 5, 8, 4, 9, 2, 2, 9, 5, 2, 5, 11, 8, 11, 5, 8, 12, 11};+unsigned edges2Row227[] = {0};+unsigned edges2Row228[] = {2, 10, 1, 4, 12, 3, 9, 5, 8, 2, 10, 1, 4, 9, 5, 4, 5, 3, 8, 3, 5, 8, 12, 3, 4, 12, 3, 2, 9, 1, 2, 8, 9, 2, 10, 5, 2, 5, 8, 9, 5, 8, 2, 12, 3, 2, 10, 12, 10, 4, 12, 10, 1, 4, 2, 10, 13, 10, 5, 13, 8, 13, 5, 8, 12, 13, 12, 3, 13, 4, 13, 3, 4, 9, 13, 9, 1, 13, 2, 13, 1, 2, 10, 13, 2, 3, 13, 12, 3, 13, 8, 12, 13, 8, 13, 5, 9, 13, 5, 4, 9, 13, 4, 13, 1, 10, 1, 13, 8, 13, 5, 10, 5, 13, 2, 10, 13, 2, 13, 3, 12, 3, 13, 4, 12, 13, 4, 13, 1, 9, 1, 13, 9, 13, 8, 2, 12, 3, 2, 5,+   12, 2, 10, 5, 8, 12, 5, 4, 9, 1, 2, 1, 3, 4, 3, 1, 4, 12, 3, 4, 8, 12, 4, 9, 8, 9, 5, 8, 9, 1, 5, 10, 5, 1, 2, 10, 1};+unsigned edges2Row229[] = {0};+unsigned edges2Row230[] = {8, 9, 5, 4, 12, 3, 8, 12, 5, 5, 12, 3, 5, 3, 4, 5, 4, 9};+unsigned edges2Row231[] = {0};+unsigned edges2Row232[] = {10, 1, 11, 11, 1, 3, 9, 5, 8, 10, 5, 11, 11, 5, 8, 11, 8, 3, 9, 3, 8, 9, 1, 3, 10, 9, 1, 10, 5, 9, 5, 8, 9, 10, 5, 11, 11, 5, 8, 11, 8, 3, 9, 3, 8, 9, 1, 3};+unsigned edges2Row233[] = {4, 1, 8, 8, 1, 5, 2, 3, 11, 4, 3, 8, 8, 3, 11, 8, 11, 5, 2, 5, 11, 2, 1, 5, 4, 2, 1, 4, 3, 2, 3, 11, 2, 4, 3, 8, 8, 3, 11, 8, 11, 5, 2, 5, 11, 2, 1, 5};+unsigned edges2Row234[] = {2, 3, 11, 9, 5, 8, 2, 3, 9, 9, 3, 8, 11, 8, 3, 2, 5, 11, 2, 9, 5, 11, 5, 8};+unsigned edges2Row235[] = {0};+unsigned edges2Row236[] = {10, 1, 2, 9, 5, 8, 10, 5, 2, 2, 5, 8, 2, 8, 9, 2, 9, 1};+unsigned edges2Row237[] = {0};+unsigned edges2Row238[] = {0};+unsigned edges2Row239[] = {0};+unsigned edges2Row240[] = {0};+unsigned edges2Row241[] = {0};+unsigned edges2Row242[] = {0};+unsigned edges2Row243[] = {0};+unsigned edges2Row244[] = {1, 2, 10, 3, 4, 12, 1, 4, 10, 10, 4, 12, 10, 12, 3, 10, 3, 2};+unsigned edges2Row245[] = {0};+unsigned edges2Row246[] = {0};+unsigned edges2Row247[] = {0};+unsigned edges2Row248[] = {0};+unsigned edges2Row249[] = {2, 3, 11, 4, 1, 9, 2, 1, 11, 11, 1, 9, 11, 9, 4, 11, 4, 3};+unsigned edges2Row250[] = {0};+unsigned edges2Row251[] = {0};+unsigned edges2Row252[] = {0};+unsigned edges2Row253[] = {0};+unsigned (*Edges2[])[] = {&edges2Row0, &edges2Row1, &edges2Row2, &edges2Row3, &edges2Row4, &edges2Row5, &edges2Row6, &edges2Row7, &edges2Row8, &edges2Row9, &edges2Row10, &edges2Row11, &edges2Row12, &edges2Row13, &edges2Row14, &edges2Row15, &edges2Row16, &edges2Row17, &edges2Row18, &edges2Row19, &edges2Row20, &edges2Row21, &edges2Row22, &edges2Row23, &edges2Row24, &edges2Row25, &edges2Row26, &edges2Row27, &edges2Row28, &edges2Row29, &edges2Row30, &edges2Row31, &edges2Row32, &edges2Row33,+   &edges2Row34, &edges2Row35, &edges2Row36, &edges2Row37, &edges2Row38, &edges2Row39, &edges2Row40, &edges2Row41, &edges2Row42, &edges2Row43, &edges2Row44, &edges2Row45, &edges2Row46, &edges2Row47, &edges2Row48, &edges2Row49, &edges2Row50, &edges2Row51, &edges2Row52, &edges2Row53, &edges2Row54, &edges2Row55, &edges2Row56, &edges2Row57, &edges2Row58, &edges2Row59, &edges2Row60, &edges2Row61, &edges2Row62, &edges2Row63, &edges2Row64, &edges2Row65, &edges2Row66, &edges2Row67, &edges2Row68,+   &edges2Row69, &edges2Row70, &edges2Row71, &edges2Row72, &edges2Row73, &edges2Row74, &edges2Row75, &edges2Row76, &edges2Row77, &edges2Row78, &edges2Row79, &edges2Row80, &edges2Row81, &edges2Row82, &edges2Row83, &edges2Row84, &edges2Row85, &edges2Row86, &edges2Row87, &edges2Row88, &edges2Row89, &edges2Row90, &edges2Row91, &edges2Row92, &edges2Row93, &edges2Row94, &edges2Row95, &edges2Row96, &edges2Row97, &edges2Row98, &edges2Row99, &edges2Row100, &edges2Row101, &edges2Row102, &edges2Row103,+   &edges2Row104, &edges2Row105, &edges2Row106, &edges2Row107, &edges2Row108, &edges2Row109, &edges2Row110, &edges2Row111, &edges2Row112, &edges2Row113, &edges2Row114, &edges2Row115, &edges2Row116, &edges2Row117, &edges2Row118, &edges2Row119, &edges2Row120, &edges2Row121, &edges2Row122, &edges2Row123, &edges2Row124, &edges2Row125, &edges2Row126, &edges2Row127, &edges2Row128, &edges2Row129, &edges2Row130, &edges2Row131, &edges2Row132, &edges2Row133, &edges2Row134, &edges2Row135, &edges2Row136,+   &edges2Row137, &edges2Row138, &edges2Row139, &edges2Row140, &edges2Row141, &edges2Row142, &edges2Row143, &edges2Row144, &edges2Row145, &edges2Row146, &edges2Row147, &edges2Row148, &edges2Row149, &edges2Row150, &edges2Row151, &edges2Row152, &edges2Row153, &edges2Row154, &edges2Row155, &edges2Row156, &edges2Row157, &edges2Row158, &edges2Row159, &edges2Row160, &edges2Row161, &edges2Row162, &edges2Row163, &edges2Row164, &edges2Row165, &edges2Row166, &edges2Row167, &edges2Row168,+   &edges2Row169, &edges2Row170, &edges2Row171, &edges2Row172, &edges2Row173, &edges2Row174, &edges2Row175, &edges2Row176, &edges2Row177, &edges2Row178, &edges2Row179, &edges2Row180, &edges2Row181, &edges2Row182, &edges2Row183, &edges2Row184, &edges2Row185, &edges2Row186, &edges2Row187, &edges2Row188, &edges2Row189, &edges2Row190, &edges2Row191, &edges2Row192, &edges2Row193, &edges2Row194, &edges2Row195, &edges2Row196, &edges2Row197, &edges2Row198, &edges2Row199, &edges2Row200,+   &edges2Row201, &edges2Row202, &edges2Row203, &edges2Row204, &edges2Row205, &edges2Row206, &edges2Row207, &edges2Row208, &edges2Row209, &edges2Row210, &edges2Row211, &edges2Row212, &edges2Row213, &edges2Row214, &edges2Row215, &edges2Row216, &edges2Row217, &edges2Row218, &edges2Row219, &edges2Row220, &edges2Row221, &edges2Row222, &edges2Row223, &edges2Row224, &edges2Row225, &edges2Row226, &edges2Row227, &edges2Row228, &edges2Row229, &edges2Row230, &edges2Row231, &edges2Row232,+   &edges2Row233, &edges2Row234, &edges2Row235, &edges2Row236, &edges2Row237, &edges2Row238, &edges2Row239, &edges2Row240, &edges2Row241, &edges2Row242, &edges2Row243, &edges2Row244, &edges2Row245, &edges2Row246, &edges2Row247, &edges2Row248, &edges2Row249, &edges2Row250, &edges2Row251, &edges2Row252, &edges2Row253};++unsigned Edges2Sizes[254] = { 1, 1, 1, 1, 18, 1, 1, 1, 1, 18, 1, 1, 1, 1, 1, 1, 1, 18, 1,+                              24, 48, 48, 1, 18, 1, 177, 1, 48, 1, 96, 1, 1, 18, 1, 1, 18,+                              177, 1, 1, 24, 48, 48, 1, 48, 96, 1, 1, 1, 1, 1, 1, 48, 96, 1,+                              1, 48, 1, 96, 1, 96, 48, 48, 1, 1, 24, 18, 48, 1, 48, 1, 1, 18,+                              48, 177, 96, 1, 1, 1, 1, 18, 48, 177, 96, 48, 96, 96, 48, 177,+                              96, 816, 177, 96, 48, 177, 18, 1, 48, 1, 1, 1, 96, 1, 1, 48,+                              96, 96, 48, 1, 48, 1, 1, 1, 1, 1, 1, 1, 48, 1, 1, 96, 48, 177,+                              18, 48, 24, 18, 1, 1, 18, 24, 48, 18, 177, 48, 96, 1, 1, 48,+                              1, 1, 1, 1, 1, 1, 1, 48, 1, 48, 96, 96, 48, 1, 1, 96, 1, 1, 1,+                              48, 1, 18, 177, 48, 96, 177, 816, 96, 177, 48, 96, 96, 48, 96,+                              177, 48, 18, 1, 1, 1, 1, 96, 177, 48, 18, 1, 1, 48, 1, 48, 18,+                              24, 1, 1, 48, 48, 96, 1, 96, 1, 48, 1, 1, 96, 48, 1, 1, 1, 1,+                              1, 1, 96, 48, 1, 48, 48, 24, 1, 1, 177, 18, 1, 1, 18, 1, 1, 96,+                              1, 48, 1, 177, 1, 18, 1, 48, 48, 24, 1, 18, 1, 1, 1, 1, 1, 1,+                              1, 18, 1, 1, 1, 1, 18, 1, 1, 1, 1};+++unsigned special_name[7] = {3, 4, 6, 7, 10, 12, 13};+unsigned special_nface[7] = {1, 1, 2, 4, 3, 3, 7};+unsigned special_nedge[7] = {18, 24, 48, 177, 96, 96, 816};+unsigned special_indSizes[7] = {2, 2, 4, 16, 8, 8, 31};++unsigned specialindRow0[] = {0, 1};+unsigned specialindRow1[] = {0, 1};+unsigned specialindRow2[] = {0, 2, 1, 3};+unsigned specialindRow3[] = {0, 8, 4, 12, 2, 10, 1, 9, 6, 14, 5, 13, 3, 11, 15, 7};+unsigned specialindRow4[] = {0, 4, 1, 5, 2, 6, 3, 7};+unsigned specialindRow5[] = {0, 4, 2, 6, 1, 5, 3, 7};+unsigned specialindRow6[] = {0, 1, 2, 4, 8, 16, 32, 3, 9, 17, 33, 6, 18, 34, 12, 20, 36, 24, 40, 35, 25, 22, 44, 19, 41, 38, 28, 83, 105, 102, 92};+unsigned (*special_ind[])[] = {&specialindRow0, &specialindRow1, &specialindRow2, &specialindRow3, &specialindRow4, &specialindRow5, &specialindRow6};++unsigned specialpos0Row0[] = {1, 2, 3, 4, 5, 6};+unsigned specialpos0Row1[] = {7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18};+unsigned (*specialpos0[])[] = {&specialpos0Row0, &specialpos0Row1};+unsigned specialpos1Row0[] = {1, 2, 3, 4, 5, 6};+unsigned specialpos1Row1[] = {7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24};+unsigned (*specialpos1[])[] = {&specialpos1Row0, &specialpos1Row1};+unsigned specialpos2Row0[] = {1, 2, 3, 4, 5, 6, 7, 8, 9};+unsigned specialpos2Row1[] = {10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33};+unsigned specialpos2Row2[] = {34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48};+unsigned specialpos2Row3[] = {34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48};+unsigned (*specialpos2[])[] = {&specialpos2Row0, &specialpos2Row1, &specialpos2Row2, &specialpos2Row3};+unsigned specialpos3Row0[] = {1, 2, 3, 4, 5, 6, 7, 8, 9};+unsigned specialpos3Row1[] = {1, 2, 3, 4, 5, 6, 7, 8, 9};+unsigned specialpos3Row2[] = {10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24};+unsigned specialpos3Row3[] = {10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24};+unsigned specialpos3Row4[] = {25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39};+unsigned specialpos3Row5[] = {25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39};+unsigned specialpos3Row6[] = {40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54};+unsigned specialpos3Row7[] = {40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54};+unsigned specialpos3Row8[] = {55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81};+unsigned specialpos3Row9[] = {55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81};+unsigned specialpos3Row10[] = {82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108};+unsigned specialpos3Row11[] = {82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108};+unsigned specialpos3Row12[] = {109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135};+unsigned specialpos3Row13[] = {109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135};+unsigned specialpos3Row14[] = {136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150};+unsigned specialpos3Row15[] = {151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 176, 177};+unsigned (*specialpos3[])[] = {&specialpos3Row0, &specialpos3Row1, &specialpos3Row2, &specialpos3Row3, &specialpos3Row4, &specialpos3Row5, &specialpos3Row6, &specialpos3Row7, &specialpos3Row8, &specialpos3Row9, &specialpos3Row10, &specialpos3Row11, &specialpos3Row12, &specialpos3Row13, &specialpos3Row14, &specialpos3Row15};+unsigned specialpos4Row0[] = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12};+unsigned specialpos4Row1[] = {13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36};+unsigned specialpos4Row2[] = {37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60};+unsigned specialpos4Row3[] = {37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60};+unsigned specialpos4Row4[] = {61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84};+unsigned specialpos4Row5[] = {61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84};+unsigned specialpos4Row6[] = {85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96};+unsigned specialpos4Row7[] = {85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96};+unsigned (*specialpos4[])[] = {&specialpos4Row0, &specialpos4Row1, &specialpos4Row2, &specialpos4Row3, &specialpos4Row4, &specialpos4Row5, &specialpos4Row6, &specialpos4Row7};+unsigned specialpos5Row0[] = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12};+unsigned specialpos5Row1[] = {13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36};+unsigned specialpos5Row2[] = {37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60};+unsigned specialpos5Row3[] = {37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60};+unsigned specialpos5Row4[] = {61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84};+unsigned specialpos5Row5[] = {61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84};+unsigned specialpos5Row6[] = {85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96};+unsigned specialpos5Row7[] = {85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96};+unsigned (*specialpos5[])[] = {&specialpos5Row0, &specialpos5Row1, &specialpos5Row2, &specialpos5Row3, &specialpos5Row4, &specialpos5Row5, &specialpos5Row6, &specialpos5Row7};+unsigned specialpos6Row0[] = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12};+unsigned specialpos6Row1[] = {13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30};+unsigned specialpos6Row2[] = {31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48};+unsigned specialpos6Row3[] = {49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66};+unsigned specialpos6Row4[] = {67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84};+unsigned specialpos6Row5[] = {85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102};+unsigned specialpos6Row6[] = {103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120};+unsigned specialpos6Row7[] = {121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150};+unsigned specialpos6Row8[] = {151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180};+unsigned specialpos6Row9[] = {181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204, 205, 206, 207, 208, 209, 210};+unsigned specialpos6Row10[] = {211, 212, 213, 214, 215, 216, 217, 218, 219, 220, 221, 222, 223, 224, 225, 226, 227, 228, 229, 230, 231, 232, 233, 234, 235, 236, 237, 238, 239, 240};+unsigned specialpos6Row11[] = {241, 242, 243, 244, 245, 246, 247, 248, 249, 250, 251, 252, 253, 254, 255, 256, 257, 258, 259, 260, 261, 262, 263, 264, 265, 266, 267, 268, 269, 270};+unsigned specialpos6Row12[] = {271, 272, 273, 274, 275, 276, 277, 278, 279, 280, 281, 282, 283, 284, 285, 286, 287, 288, 289, 290, 291, 292, 293, 294, 295, 296, 297, 298, 299, 300};+unsigned specialpos6Row13[] = {301, 302, 303, 304, 305, 306, 307, 308, 309, 310, 311, 312, 313, 314, 315, 316, 317, 318, 319, 320, 321, 322, 323, 324, 325, 326, 327, 328, 329, 330};+unsigned specialpos6Row14[] = {331, 332, 333, 334, 335, 336, 337, 338, 339, 340, 341, 342, 343, 344, 345, 346, 347, 348, 349, 350, 351, 352, 353, 354, 355, 356, 357, 358, 359, 360};+unsigned specialpos6Row15[] = {361, 362, 363, 364, 365, 366, 367, 368, 369, 370, 371, 372, 373, 374, 375, 376, 377, 378, 379, 380, 381, 382, 383, 384, 385, 386, 387, 388, 389, 390};+unsigned specialpos6Row16[] = {391, 392, 393, 394, 395, 396, 397, 398, 399, 400, 401, 402, 403, 404, 405, 406, 407, 408, 409, 410, 411, 412, 413, 414, 415, 416, 417, 418, 419, 420};+unsigned specialpos6Row17[] = {421, 422, 423, 424, 425, 426, 427, 428, 429, 430, 431, 432, 433, 434, 435, 436, 437, 438, 439, 440, 441, 442, 443, 444, 445, 446, 447, 448, 449, 450};+unsigned specialpos6Row18[] = {451, 452, 453, 454, 455, 456, 457, 458, 459, 460, 461, 462, 463, 464, 465, 466, 467, 468, 469, 470, 471, 472, 473, 474, 475, 476, 477, 478, 479, 480};+unsigned specialpos6Row19[] = {481, 482, 483, 484, 485, 486, 487, 488, 489, 490, 491, 492, 493, 494, 495, 496, 497, 498, 499, 500, 501, 502, 503, 504, 505, 506, 507, 508, 509, 510, 511, 512, 513, 514, 515, 516};+unsigned specialpos6Row20[] = {517, 518, 519, 520, 521, 522, 523, 524, 525, 526, 527, 528, 529, 530, 531, 532, 533, 534, 535, 536, 537, 538, 539, 540, 541, 542, 543, 544, 545, 546, 547, 548, 549, 550, 551, 552};+unsigned specialpos6Row21[] = {553, 554, 555, 556, 557, 558, 559, 560, 561, 562, 563, 564, 565, 566, 567, 568, 569, 570, 571, 572, 573, 574, 575, 576, 577, 578, 579, 580, 581, 582, 583, 584, 585, 586, 587, 588};+unsigned specialpos6Row22[] = {589, 590, 591, 592, 593, 594, 595, 596, 597, 598, 599, 600, 601, 602, 603, 604, 605, 606, 607, 608, 609, 610, 611, 612, 613, 614, 615, 616, 617, 618, 619, 620, 621, 622, 623, 624};+unsigned specialpos6Row23[] = {625, 626, 627, 628, 629, 630, 631, 632, 633, 634, 635, 636, 637, 638, 639, 640, 641, 642};+unsigned specialpos6Row24[] = {643, 644, 645, 646, 647, 648, 649, 650, 651, 652, 653, 654, 655, 656, 657, 658, 659, 660};+unsigned specialpos6Row25[] = {661, 662, 663, 664, 665, 666, 667, 668, 669, 670, 671, 672, 673, 674, 675, 676, 677, 678};+unsigned specialpos6Row26[] = {679, 680, 681, 682, 683, 684, 685, 686, 687, 688, 689, 690, 691, 692, 693, 694, 695, 696};+unsigned specialpos6Row27[] = {697, 698, 699, 700, 701, 702, 703, 704, 705, 706, 707, 708, 709, 710, 711, 712, 713, 714, 715, 716, 717, 718, 719, 720, 721, 722, 723, 724, 725, 726};+unsigned specialpos6Row28[] = {727, 728, 729, 730, 731, 732, 733, 734, 735, 736, 737, 738, 739, 740, 741, 742, 743, 744, 745, 746, 747, 748, 749, 750, 751, 752, 753, 754, 755, 756};+unsigned specialpos6Row29[] = {757, 758, 759, 760, 761, 762, 763, 764, 765, 766, 767, 768, 769, 770, 771, 772, 773, 774, 775, 776, 777, 778, 779, 780, 781, 782, 783, 784, 785, 786};+unsigned specialpos6Row30[] = {787, 788, 789, 790, 791, 792, 793, 794, 795, 796, 797, 798, 799, 800, 801, 802, 803, 804, 805, 806, 807, 808, 809, 810, 811, 812, 813, 814, 815, 816};+unsigned (*specialpos6[])[] = {&specialpos6Row0, &specialpos6Row1, &specialpos6Row2, &specialpos6Row3, &specialpos6Row4, &specialpos6Row5, &specialpos6Row6, &specialpos6Row7, &specialpos6Row8, &specialpos6Row9, &specialpos6Row10, &specialpos6Row11, &specialpos6Row12, &specialpos6Row13, &specialpos6Row14, &specialpos6Row15, &specialpos6Row16, &specialpos6Row17, &specialpos6Row18, &specialpos6Row19, &specialpos6Row20, &specialpos6Row21, &specialpos6Row22, &specialpos6Row23, &specialpos6Row24,+   &specialpos6Row25, &specialpos6Row26, &specialpos6Row27, &specialpos6Row28, &specialpos6Row29, &specialpos6Row30};+unsigned (*(*special_pos[])[])[] = {&specialpos0, &specialpos1, &specialpos2, &specialpos3, &specialpos4, &specialpos5, &specialpos6};++unsigned specialposSize0Row0 = 6;+unsigned specialposSize0Row1 = 12;+unsigned* specialposSize0[] = {&specialposSize0Row0, &specialposSize0Row1};+unsigned specialposSize1Row0 = 6;+unsigned specialposSize1Row1 = 18;+unsigned* specialposSize1[] = {&specialposSize1Row0, &specialposSize1Row1};+unsigned specialposSize2Row0 = 9;+unsigned specialposSize2Row1 = 24;+unsigned specialposSize2Row2 = 15;+unsigned specialposSize2Row3 = 15;+unsigned* specialposSize2[] = {&specialposSize2Row0, &specialposSize2Row1, &specialposSize2Row2, &specialposSize2Row3};+unsigned specialposSize3Row0 = 9;+unsigned specialposSize3Row1 = 9;+unsigned specialposSize3Row2 = 15;+unsigned specialposSize3Row3 = 15;+unsigned specialposSize3Row4 = 15;+unsigned specialposSize3Row5 = 15;+unsigned specialposSize3Row6 = 15;+unsigned specialposSize3Row7 = 15;+unsigned specialposSize3Row8 = 27;+unsigned specialposSize3Row9 = 27;+unsigned specialposSize3Row10 = 27;+unsigned specialposSize3Row11 = 27;+unsigned specialposSize3Row12 = 27;+unsigned specialposSize3Row13 = 27;+unsigned specialposSize3Row14 = 15;+unsigned specialposSize3Row15 = 27;+unsigned* specialposSize3[] = {&specialposSize3Row0, &specialposSize3Row1, &specialposSize3Row2, &specialposSize3Row3, &specialposSize3Row4, &specialposSize3Row5, &specialposSize3Row6, &specialposSize3Row7, &specialposSize3Row8, &specialposSize3Row9, &specialposSize3Row10, &specialposSize3Row11, &specialposSize3Row12, &specialposSize3Row13, &specialposSize3Row14, &specialposSize3Row15};+unsigned specialposSize4Row0 = 12;+unsigned specialposSize4Row1 = 24;+unsigned specialposSize4Row2 = 24;+unsigned specialposSize4Row3 = 24;+unsigned specialposSize4Row4 = 24;+unsigned specialposSize4Row5 = 24;+unsigned specialposSize4Row6 = 12;+unsigned specialposSize4Row7 = 12;+unsigned* specialposSize4[] = {&specialposSize4Row0, &specialposSize4Row1, &specialposSize4Row2, &specialposSize4Row3, &specialposSize4Row4, &specialposSize4Row5, &specialposSize4Row6, &specialposSize4Row7};+unsigned specialposSize5Row0 = 12;+unsigned specialposSize5Row1 = 24;+unsigned specialposSize5Row2 = 24;+unsigned specialposSize5Row3 = 24;+unsigned specialposSize5Row4 = 24;+unsigned specialposSize5Row5 = 24;+unsigned specialposSize5Row6 = 12;+unsigned specialposSize5Row7 = 12;+unsigned* specialposSize5[] = {&specialposSize5Row0, &specialposSize5Row1, &specialposSize5Row2, &specialposSize5Row3, &specialposSize5Row4, &specialposSize5Row5, &specialposSize5Row6, &specialposSize5Row7};+unsigned specialposSize6Row0 = 12;+unsigned specialposSize6Row1 = 18;+unsigned specialposSize6Row2 = 18;+unsigned specialposSize6Row3 = 18;+unsigned specialposSize6Row4 = 18;+unsigned specialposSize6Row5 = 18;+unsigned specialposSize6Row6 = 18;+unsigned specialposSize6Row7 = 30;+unsigned specialposSize6Row8 = 30;+unsigned specialposSize6Row9 = 30;+unsigned specialposSize6Row10 = 30;+unsigned specialposSize6Row11 = 30;+unsigned specialposSize6Row12 = 30;+unsigned specialposSize6Row13 = 30;+unsigned specialposSize6Row14 = 30;+unsigned specialposSize6Row15 = 30;+unsigned specialposSize6Row16 = 30;+unsigned specialposSize6Row17 = 30;+unsigned specialposSize6Row18 = 30;+unsigned specialposSize6Row19 = 36;+unsigned specialposSize6Row20 = 36;+unsigned specialposSize6Row21 = 36;+unsigned specialposSize6Row22 = 36;+unsigned specialposSize6Row23 = 18;+unsigned specialposSize6Row24 = 18;+unsigned specialposSize6Row25 = 18;+unsigned specialposSize6Row26 = 18;+unsigned specialposSize6Row27 = 30;+unsigned specialposSize6Row28 = 30;+unsigned specialposSize6Row29 = 30;+unsigned specialposSize6Row30 = 30;+unsigned* specialposSize6[] = {&specialposSize6Row0, &specialposSize6Row1, &specialposSize6Row2, &specialposSize6Row3, &specialposSize6Row4, &specialposSize6Row5, &specialposSize6Row6, &specialposSize6Row7, &specialposSize6Row8, &specialposSize6Row9, &specialposSize6Row10, &specialposSize6Row11, &specialposSize6Row12, &specialposSize6Row13, &specialposSize6Row14, &specialposSize6Row15, &specialposSize6Row16, &specialposSize6Row17, &specialposSize6Row18, &specialposSize6Row19, &specialposSize6Row20, &specialposSize6Row21, &specialposSize6Row22, &specialposSize6Row23, &specialposSize6Row24, &specialposSize6Row25, &specialposSize6Row26, &specialposSize6Row27, &specialposSize6Row28, &specialposSize6Row29, &specialposSize6Row30};+unsigned* (*(special_posSize)[])[] = {&specialposSize0, &specialposSize1, &specialposSize2, &specialposSize3, &specialposSize4, &specialposSize5, &specialposSize6};

@@ -0,0 +1,407 @@+#ifndef __HEADER__+#define __HEADER__+#include "marchingCubes.h"+#endif++void displayMatrix_u(unsigned** M, unsigned nx, unsigned ny) {+  for(unsigned i = 0; i < nx; i++) {+    for(unsigned j = 0; j < ny; j++) {+      printf("M[%u][%u]=%u ", i, j, M[i][j]);+    }+    printf("\n");+  }+}++void displayMatrix(double** M, unsigned nx, unsigned ny) {+  for(unsigned i = 0; i < nx; i++) {+    for(unsigned j = 0; j < ny; j++) {+      printf("M[%u][%u]=%f ", i, j, M[i][j]);+    }+    printf("\n");+  }+}++void freeMatrix_u(unsigned** M, unsigned nx) {+  for(unsigned i = 0; i < nx; i++) {+    free(M[i]);+  }+  free(M);+}++void freeMatrix_s(size_t** M, size_t nx) {+  for(size_t i = 0; i < nx; i++) {+    free(M[i]);+  }+  free(M);+}++void freeMatrix_d(double** M, unsigned nx) {+  for(unsigned i = 0; i < nx; i++) {+    free(M[i]);+  }+  free(M);+}++size_t** minorMatrix(size_t** M,+                     unsigned m,+                     unsigned n,+                     unsigned row,+                     unsigned col) {+  size_t** out = malloc((m - 1) * sizeof(size_t*));+  for(unsigned i = 0; i < (m - 1); i++) {+    out[i] = malloc((n - 1) * sizeof(size_t));+    unsigned ii = i < row ? i : i + 1;+    for(unsigned j = 0; j < (n - 1); j++) {+      unsigned jj = j < col ? j : j + 1;+      out[i][j] = M[ii][jj];+    }+  }+  return out;+}++size_t** scaleMatrix(size_t k, size_t** M, unsigned m, unsigned n) {+  size_t** out = malloc(m * sizeof(size_t*));+  for(unsigned i = 0; i < m; i++) {+    out[i] = malloc(n * sizeof(size_t));+    for(unsigned j = 0; j < n; j++) {+      out[i][j] = k * M[i][j];+    }+  }+  return out;+}++size_t** scaleMinorMatrix(size_t k,+                          size_t** M,+                          unsigned m,+                          unsigned n,+                          unsigned row,+                          unsigned col) {+  size_t** out = malloc((m - 1) * sizeof(size_t*));+  for(unsigned i = 0; i < (m - 1); i++) {+    out[i] = malloc((n - 1) * sizeof(size_t));+    unsigned ii = i < row ? i : i + 1;+    for(unsigned j = 0; j < (n - 1); j++) {+      unsigned jj = j < col ? j : j + 1;+      out[i][j] = k * M[ii][jj];+    }+  }+  return out;+}++size_t** matricialSum(size_t** M1, size_t** M2, size_t m, unsigned n) {+  size_t** out = malloc(m * sizeof(size_t*));+  for(size_t i = 0; i < m; i++) {+    out[i] = malloc(n * sizeof(size_t));+    for(unsigned j = 0; j < n; j++) {+      out[i][j] = M1[i][j] + M2[i][j];+    }+  }+  return out;+}++size_t** levelMatrix(double** M,+                     unsigned m,+                     unsigned n,+                     double level,+                     unsigned strict) {+  size_t** out = malloc(m * sizeof(size_t*));+  for(unsigned i = 0; i < m; i++) {+    out[i] = malloc(n * sizeof(size_t));+    for(unsigned j = 0; j < n; j++) {+      out[i][j] = (strict ? (M[i][j] > level) : (M[i][j] >= level)) ? 1 : 0;+    }+  }+  return out;+}++double** toMatrix(double*** A, unsigned m, unsigned n, unsigned k) {+  double** out = malloc(m * sizeof(double*));+  for(unsigned i = 0; i < m; i++) {+    out[i] = malloc(n * sizeof(double));+    for(unsigned j = 0; j < n; j++) {+      out[i][j] = A[i][j][k];+    }+  }+  return out;+}++unsigned** whichIndicesAndItems(size_t** M,+                                unsigned m,+                                unsigned n,+                                unsigned* outlength) {+  unsigned** out = malloc(2 * sizeof(unsigned*));+  out[0] = malloc(m * n * sizeof(unsigned));+  out[1] = malloc(m * n * sizeof(unsigned));+  unsigned k = 0;+  unsigned ij = 0;+  for(unsigned j = 0; j < n; j++) {+    for(unsigned i = 0; i < m; i++) {+      if(M[i][j] > 0 && M[i][j] < 255) {+        out[0][k] = ij;+        out[1][k] = M[i][j];+        k++;+      }+      ij++;+    }+  }+  *outlength = k;+  return out;+}++void freeArray(double*** A, unsigned nx, unsigned ny) {+  for(unsigned i = 0; i < nx; i++) {+    for(unsigned j = 0; j < ny; j++) {+      free(A[i][j]);+    }+    free(A[i]);+  }+  free(A);+}++size_t** kro1(size_t** M, unsigned nx, unsigned ny, size_t n) {+  size_t** out = malloc((nx * n + 1) * sizeof(size_t*));+  for(size_t i = 0; i < nx * n; i++) {+    out[i] = malloc(ny * sizeof(size_t));+    for(unsigned j = 0; j < ny; j++) {+      out[i][j] = M[i % nx][j];+    }+  }+  out[nx * n] = malloc(ny * sizeof(size_t));+  for(unsigned j = 0; j < ny; j++) {+    out[nx * n][j] = 0;+  }+  return out;+}++size_t** kro2(size_t** M, size_t nx, unsigned ny, unsigned n) {+  size_t** out = malloc((nx * n + 1) * sizeof(size_t*));+  unsigned* replicates = malloc((nx * n) * sizeof(unsigned));+  unsigned count = 0;+  for(size_t i = 0; i < nx; i++) {+    for(unsigned m = 0; m < n; m++) {+      replicates[count] = i;+      count++;+    }+  }+  for(size_t i = 0; i < nx * n; i++) {+    out[i] = malloc(ny * sizeof(size_t));+    for(unsigned j = 0; j < ny; j++) {+      out[i][j] = M[replicates[i]][j];+    }+  }+  free(replicates);+  out[nx * n] = malloc(ny * sizeof(size_t));+  for(unsigned j = 0; j < ny; j++) {+    out[nx * n][j] = 0;+  }+  return out;+}++size_t* replicate(size_t* x, unsigned* counts, size_t n) {+  unsigned total = 0;+  for(size_t i = 0; i < n; i++) {+    total += counts[i];+  }+  size_t* out = malloc(total * sizeof(size_t));+  unsigned count = 0;+  for(size_t i = 0; i < n; i++) {+    for(unsigned j = 0; j < counts[i]; j++) {+      out[count] = x[i];+      count++;+    }+  }+  return out;+}++// concatenate rows of a jagged array+int* unlist(int (*jagged[])[],+            unsigned* jaggedsizes,+            unsigned* rows,+            unsigned nrows,+            unsigned* outlength) {+  unsigned count = 0;+  for(unsigned i = 0; i < nrows; i++) {+    count += jaggedsizes[rows[i]];+  }+  *outlength = count;+  unsigned* out = malloc(count * sizeof(unsigned));+  count = 0;+  for(unsigned i = 0; i < nrows; i++) {+    for(unsigned j = 0; j < jaggedsizes[rows[i]]; j++) {+      out[count] = (*jagged[rows[i]])[j];+      count++;+    }+  }+  return out;+}++unsigned* unlist_u(unsigned (*jagged[])[],+                   unsigned* jaggedsizes,+                   unsigned* rows,+                   unsigned nrows,+                   unsigned* outlength) {+  unsigned count = 0;+  for(unsigned i = 0; i < nrows; i++) {+    count += jaggedsizes[rows[i]];+  }+  *outlength = count;+  unsigned* out = malloc(count * sizeof(unsigned));+  count = 0;+  for(unsigned i = 0; i < nrows; i++) {+    for(unsigned j = 0; j < jaggedsizes[rows[i]]; j++) {+      out[count] = (*jagged[rows[i]])[j];+      count++;+    }+  }+  return out;+}++unsigned upow(unsigned base, unsigned exp) {+  unsigned result = 1;+  while(exp) {+    if(exp & 1)+      result *= base;+    exp >>= 1;+    base *= base;+  }+  return result;+}++// R: matrix(M, ncol=ncol, byrow=TRUE)+unsigned** vector2matrix(unsigned* M, unsigned lengthM, unsigned ncol) {+  unsigned nrow = lengthM / ncol;+  unsigned** out = malloc(nrow * sizeof(unsigned*));+  for(unsigned i = 0; i < nrow; i++) {+    out[i] = malloc(ncol * sizeof(unsigned));+    for(unsigned j = 0; j < ncol; j++) {+      out[i][j] = M[i * ncol + j];+    }+  }+  return out;+}++// cbind(matrix, vector, vector)+unsigned** cbind(unsigned** M,+                 size_t* v1,+                 unsigned* v2,+                 unsigned nrow,+                 unsigned ncol) {+  unsigned** out = malloc(nrow * sizeof(unsigned*));+  for(unsigned i = 0; i < nrow; i++) {+    out[i] = malloc((ncol + 2) * sizeof(unsigned));+    for(unsigned j = 0; j < ncol; j++) {+      out[i][j] = M[i][j];+    }+    out[i][ncol] = v1[i];+    out[i][ncol + 1] = v2[i];+  }+  return out;+}++unsigned** subsetMatrix(unsigned** M,+                        unsigned* rows,+                        unsigned* cols,+                        unsigned nrows,+                        unsigned ncols) {+  unsigned** out = malloc(nrows * sizeof(unsigned*));+  for(unsigned i = 0; i < nrows; i++) {+    out[i] = malloc(ncols * sizeof(unsigned));+    for(unsigned j = 0; j < ncols; j++) {+      out[i][j] = M[rows[i]][cols[j]];+    }+  }+  return out;+}++unsigned* whichEqual(unsigned* v,+                     unsigned x,+                     unsigned lv,+                     unsigned* outlength) {+  short* flag = malloc(lv * sizeof(short));+  unsigned count = 0;+  for(unsigned i = 0; i < lv; i++) {+    if(v[i] == x) {+      flag[i] = 1;+      count++;+    } else {+      flag[i] = 0;+    }+  }+  *outlength = count;+  unsigned* out = malloc(count * sizeof(unsigned));+  count = 0;+  for(unsigned i = 0; i < lv; i++) {+    if(flag[i]) {+      out[count] = i;+      count++;+    }+  }+  free(flag);+  return (out);+}++unsigned* matrix2vectorMinusFirstColumn(unsigned** M,+                                        unsigned nrow,+                                        unsigned ncol) {+  unsigned* out = malloc(nrow * (ncol - 1) * sizeof(unsigned));+  for(size_t i = 0; i < nrow; i++) {+    for(size_t j = 1; j < ncol; j++) {+      out[i * (ncol - 1) + j - 1] = M[i][j];+    }+  }+  return out;+}++size_t* repeach(size_t* x, unsigned times, size_t n) {+  size_t* out = malloc(n * times * sizeof(size_t));+  unsigned count = 0;+  for(size_t i = 0; i < n; i++) {+    for(unsigned j = 0; j < times; j++) {+      out[count] = x[i];+      count++;+    }+  }+  return out;+}++int* jthColumn(int* M, unsigned lengthM, unsigned ncol, unsigned j) {+  unsigned nrow = lengthM / ncol;+  int* out = malloc(nrow * sizeof(unsigned));+  for(unsigned i = 0; i < nrow; i++) {+    out[i] = M[i * ncol + j];+  }+  return out;+}++double** transpose(double** M, size_t nrow, size_t ncol) {+  double** out = malloc(ncol * sizeof(double*));+  for(size_t i = 0; i < ncol; i++) {+    out[i] = malloc(nrow * sizeof(double));+    for(size_t j = 0; j < nrow; j++) {+      out[i][j] = M[j][i];+    }+  }+  return out;+}++unsigned* vectorialSum(unsigned* v1, unsigned* v2, unsigned l) {+  unsigned* out = malloc(l * sizeof(unsigned));+  for(unsigned i = 0; i < l; i++) {+    out[i] = v1[i] + v2[i];+  }+  return out;+}++double*** vector2array(double* v, unsigned nx, unsigned ny, unsigned nz) {+  double*** out = malloc(nx * sizeof(double**));+  for(size_t i = 0; i < nx; i++) {+    out[i] = malloc(ny * sizeof(double*));+    for(size_t j = 0; j < ny; j++) {+      out[i][j] = malloc(nz * sizeof(double));+      for(size_t k = 0; k < nz; k++) {+        out[i][j][k] = v[i * nz * ny + j * nz + k];+      }+    }+  }+  return out;+}

@@ -0,0 +1,8 @@+# Changelog for `marching-cubes2`+++## 0.1.0.0 - 2023-02-23++First release.++

@@ -0,0 +1,30 @@+Copyright Stéphane Laurent (c) 2023++All rights reserved.++Redistribution and use in source and binary forms, with or without+modification, are permitted provided that the following conditions are met:++    * Redistributions of source code must retain the above copyright+      notice, this list of conditions and the following disclaimer.++    * Redistributions in binary form must reproduce the above+      copyright notice, this list of conditions and the following+      disclaimer in the documentation and/or other materials provided+      with the distribution.++    * Neither the name of Stéphane Laurent nor the names of other+      contributors may be used to endorse or promote products derived+      from this software without specific prior written permission.++THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS+"AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT+LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR+A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT+OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,+SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT+LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,+DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY+THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT+(INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE+OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.

@@ -0,0 +1,12 @@+# marching-cubes2++<!-- badges: start -->+[![Stack-lts](https://github.com/stla/marching-cubes2/actions/workflows/Stack-lts.yml/badge.svg)](https://github.com/stla/marching-cubes2/actions/workflows/Stack-lts.yml)+[![Stack-nightly](https://github.com/stla/marching-cubes2/actions/workflows/Stack-nightly.yml/badge.svg)](https://github.com/stla/marching-cubes2/actions/workflows/Stack-nightly.yml)+<!-- badges: end -->+++Haskell and C implementation of the marching cubes algorithm. +Faster than the 'marching-cubes' library. +See [some examples here](https://github.com/stla/isosurfaces2).+

@@ -0,0 +1,2 @@+import Distribution.Simple+main = defaultMain

@@ -0,0 +1,50 @@+name:                marching-cubes2+version:             0.1.0.0+synopsis:            Marching Cubes+description:         The marching cubes algorithm.+homepage:            https://github.com/stla/marching-cubes2#readme+license:             BSD3+license-file:        LICENSE+author:              Stéphane Laurent+maintainer:          laurent_step@outlook.fr+copyright:           2023 Stéphane Laurent+category:            Computational geometry+build-type:          Simple+extra-source-files:  README.md+                     CHANGELOG.md+cabal-version:       >=1.10++library+  hs-source-dirs:      src+  exposed-modules:     MarchingCubes+  other-modules:       MarchingCubes.MarchingCubes+                     , MarchingCubes.Mesh+  other-extensions:    ForeignFunctionInterface+                     , FlexibleContexts+                     , TupleSections+  include-dirs:        C+  C-sources:           C/utils.c+                     , C/tables.c+                     , C/marchingCubes.c+  install-includes:    C/marchingCubes.h+  build-depends:       base >= 4.7 && < 5+                     , vector >= 0.12.3.1+                     , matrix >= 0.3.6.1+                     , extra >= 1.7.9+                     , split >= 0.2.3.4+                     , linear >= 1.21.6+                     , containers >= 0.6.4.1+  default-language:    Haskell2010+  ghc-options:         -Wall+                       -Wcompat+                       -Widentities+                       -Wincomplete-record-updates+                       -Wincomplete-uni-patterns+                       -Wmissing-export-lists+                       -Wmissing-home-modules+                       -Wpartial-fields+                       -Wredundant-constraints++source-repository head+  type:     git+  location: https://github.com/stla/marching-cubes2

@@ -0,0 +1,12 @@+module MarchingCubes+  ( module X+  ) where+import           MarchingCubes.MarchingCubes   as X+                                                ( Voxel+                                                , XYZ+                                                , makeVoxel+                                                )+import           MarchingCubes.Mesh            as X+                                                ( Mesh (..)+                                                , makeMesh+                                                )

@@ -0,0 +1,117 @@+{-# LANGUAGE ForeignFunctionInterface #-}+{-# LANGUAGE FlexibleContexts         #-}+module MarchingCubes.MarchingCubes+  ( Voxel+  , XYZ+  , marchingCubes+  , makeVoxel+  ) where+import           Control.Monad                 (when)+import           Data.List                     (transpose)+import           Foreign.C.Types+import           Foreign.Marshal.Alloc         (free, mallocBytes)+import           Foreign.Marshal.Array         (peekArray, pokeArray)+import           Foreign.Ptr                   (Ptr)+import           Foreign.Storable              (peek, sizeOf)+import           Data.Matrix                   ( Matrix+                                               , mapCol+                                               , fromLists )++type Bounds a = ((a, a), (a, a), (a, a))+type Dims = (Int, Int, Int)+type Voxel a = ([a], (Bounds a, Dims))+type XYZ a = (a, a, a)++-- | Make the voxel. +makeVoxel+  :: RealFloat a+  => (XYZ a -> a) -- ^ the function defining the isosurface+  -> Bounds a     -- ^ bounds of the grid+  -> Dims         -- ^ numbers of subdivisions of the grid+  -> Voxel a+makeVoxel fun bds@((xm, xM), (ym, yM), (zm, zM)) dims@(nx, ny, nz) =+  (values, (bds, dims))+  where+    x_ = [ xm + (xM - xm) * fracx i | i <- [0 .. nx - 1] ]+    fracx p = realToFrac p / (realToFrac nx - 1)+    y_ = [ ym + (yM - ym) * fracy i | i <- [0 .. ny - 1] ]+    fracy p = realToFrac p / (realToFrac ny - 1)+    z_ = [ zm + (zM - zm) * fracz i | i <- [0 .. nz - 1] ]+    fracz p = realToFrac p / (realToFrac nz - 1)+    values = [ fun (x, y, z) | x <- x_, y <- y_, z <- z_ ]++voxelMax :: RealFloat a => Voxel a -> a+voxelMax (values, _) = maximum (filter (not . isNaN) values)++rescale :: Fractional a => (a, a) -> Int -> a -> a+rescale (minmm, maxmm) n w = minmm + (maxmm - minmm) * w / fromIntegral (n - 1)++rescaleMatrix :: Fractional a => Matrix a -> Bounds a -> Dims -> Matrix a+rescaleMatrix mtrx (xbds, ybds, zbds) (nx, ny, nz) = mtrx'''+ where+  mtrx'   = mapCol (\_ w -> rescale xbds nx (w - 1)) 1 mtrx+  mtrx''  = mapCol (\_ w -> rescale ybds ny (w - 1)) 2 mtrx'+  mtrx''' = mapCol (\_ w -> rescale zbds nz (w - 1)) 3 mtrx''++foreign import ccall unsafe "computeContour3d" c_computeContour3d+    :: Ptr CDouble+    -> CUInt+    -> CUInt+    -> CUInt+    -> CDouble+    -> CDouble+    -> Ptr CSize+    -> IO (Ptr (Ptr CDouble))++computeContour3d :: +  Voxel Double -> Double -> IO (Ptr (Ptr CDouble), Int)+computeContour3d voxel level = do +  let voxmax = voxelMax voxel+      (nx, ny, nz) = snd $ snd voxel +      voxel' = map realToFrac (fst voxel)+  nrowsPtr <- mallocBytes (sizeOf (undefined :: CSize))+  putStrLn "Allocating voxel"+  voxelPtr <- mallocBytes (nx*ny*nz * sizeOf (undefined :: CDouble))+  putStrLn "Poking voxel..."+  pokeArray voxelPtr voxel'+  putStrLn "Run C code..."+  result <- c_computeContour3d voxelPtr+            (fromIntegral nx) (fromIntegral ny) (fromIntegral nz)+            (realToFrac voxmax) (realToFrac level) nrowsPtr+  nrows <- peek nrowsPtr+  free nrowsPtr+  free voxelPtr+  return (result, fromIntegral nrows)++marchingCubes :: Voxel Double -> Double -> Bool -> IO (Matrix Double)+marchingCubes voxel level summary = do+  let (bds@(xbds, ybds, zbds), dims@(nx, ny, nz)) = snd voxel+  (ppCDouble, nrows) <- computeContour3d voxel level+  cpoints <- mapM (peekArray 3) =<< peekArray nrows ppCDouble+  let toDouble :: CDouble -> Double+      toDouble = realToFrac+      points = map (map toDouble) cpoints+  when summary $ do+    let tpoints = transpose points+        x = tpoints !! 0+        y = tpoints !! 1+        z = tpoints !! 2+        xm = minimum x - 1+        xM = maximum x - 1+        ym = minimum y - 1+        yM = maximum y - 1+        zm = minimum z - 1+        zM = maximum z - 1+        xmin = rescale xbds nx xm+        xmax = rescale xbds nx xM+        ymin = rescale ybds ny ym+        ymax = rescale ybds ny yM+        zmin = rescale zbds nz zm+        zmax = rescale zbds nz zM+    putStrLn "Prebounds:"+    print ((xm, ym, zm), (xM, yM, zM))+    putStrLn "Bounds:"+    print ((xmin, ymin, zmin), (xmax, ymax, zmax))+  return $ rescaleMatrix (fromLists points) bds dims++

@@ -0,0 +1,123 @@+{-# LANGUAGE FlexibleContexts #-}+{-# LANGUAGE TupleSections    #-}+module MarchingCubes.Mesh+  ( makeMesh+  , Mesh (..)+  ) where+import           Data.Bifunctor                 ( bimap )+import           Data.Foldable                  ( foldl' )+import           Data.IntMap.Strict             ( IntMap+                                                , elems+                                                , unionWith+                                                )+import qualified Data.IntMap.Strict            as IM+import           Data.List.Split                ( chunksOf )+import           Data.Map.Strict                ( insertLookupWithKey )+import qualified Data.Map.Strict               as M+import           Data.Matrix                    ( Matrix+                                                , nrows+                                                , toLists+                                                )+import           Data.Tuple.Extra               ( (***)+                                                , curry3+                                                , second+                                                , thd3+                                                )+import           Data.Vector.Unboxed            ( (!)+                                                , Unbox+                                                , Vector+                                                , fromList+                                                )+import           Linear                         ( V3(..)+                                                , (^+^)+                                                , (^-^)+                                                , cross+                                                , signorm+                                                )+import           MarchingCubes.MarchingCubes    ( Voxel+                                                , marchingCubes+                                                )+++data Mesh a = Mesh +  { +    _vertices :: Vector (a, a, a)+  , _faces :: [(Int, Int, Int)] +  , _normals :: Vector (a, a, a)+  } +  deriving Show++matrix2listOfTriples :: Matrix a -> [(a, a, a)]+matrix2listOfTriples mtrx =+  map (\row -> (row !! 0, row !! 1, row !! 2)) (toLists mtrx)++uniqueWithIndices :: (Unbox a, Ord a) => [a] -> (Vector a, Vector Int)+uniqueWithIndices list = (fromList *** fromList) $ go 0 M.empty list+ where+  go _ _ []       = ([], [])+  go i m (x : xs) = case insertLookupWithKey (curry3 thd3) x i m of+    (Nothing, m') -> bimap (x :) (i :) (go (i + 1) m' xs)+    (Just j , m') -> second (j :) (go i m' xs)++undupMesh+  :: (Unbox a, Ord a) => ([(a, a, a)], [[Int]]) -> (Vector (a, a, a), [[Int]])+undupMesh (vertices, faces) = (vertices', faces')+ where+  (vertices', idx) = uniqueWithIndices vertices+  faces'           = map (map (idx !)) faces++normal :: Floating a => (a, a, a) -> (a, a, a) -> (a, a, a) -> V3 a+normal v1 v2 v3 = signorm $ cross (v3' ^-^ v1') (v2' ^-^ v1') -- (NaN,NaN,NaN) if degenerate face+ where+  toV3 (x, y, z) = V3 x y z+  v1' = toV3 v1+  v2' = toV3 v2+  v3' = toV3 v3++faceNormals+  :: (Floating a, Unbox a) => Vector (a, a, a) -> [Int] -> IntMap (V3 a)+faceNormals vs face = IM.fromList nrmls+ where+  nrml  = normal (vs ! (face !! 0)) (vs ! (face !! 1)) (vs ! (face !! 2))+  nrmls = map (, nrml) face++normalsV3+  :: (Floating a, Unbox a) => (Vector (a, a, a), [[Int]]) -> IntMap (V3 a)+normalsV3 (vs, faces) = IM.map+  signorm+  (foldl' (unionWith (^+^)) IM.empty (map (faceNormals vs) faces))++normals+  :: (Floating a, Unbox a) => (Vector (a, a, a), [[Int]]) -> Vector (a, a, a)+normals mesh = fromList $ map fromV3 (elems $ normalsV3 mesh)+  where fromV3 (V3 x y z) = (x, y, z)++degenerateFace :: (Unbox a, Eq a) => Vector (a, a, a) -> [Int] -> Bool+degenerateFace vertices face = v1 == v2 || v1 == v3 || v2 == v3+ where+  v1 = vertices ! (face !! 0)+  v2 = vertices ! (face !! 1)+  v3 = vertices ! (face !! 2)++-- | Make mesh of isosurface.+makeMesh :: +     Voxel Double -- ^ voxel obtained with `makeVoxel`+  -> Double       -- ^ isovalue+  -> IO (Mesh Double)+makeMesh voxel level = do+  mtrx <- marchingCubes voxel level True+  let vertices            = matrix2listOfTriples mtrx+      faces               = chunksOf 3 [0 .. nrows mtrx - 1]+      (vertices', faces') = undupMesh (vertices, faces)+      faces''             = filter (not . degenerateFace vertices') faces'+      toTriplet ikj       = (ikj !! 0, ikj !! 2, ikj !! 1)+      faces'''            = map toTriplet faces''+      mesh                = (vertices', faces'')+  return $ Mesh +    {+      _vertices = vertices', _faces = faces''', _normals = normals mesh+    }++++