hgp canonical codeword fixed

Author: mkangquantum

doc/intro.ipynb

@@ -331,14 +331,6 @@
     "print('# layers of entangling gates: ', depth)"
    ]
   },
-  {
-   "cell_type": "markdown",
-   "id": "27121110",
-   "metadata": {},
-   "source": [
-    "For the BPC code, logical operators are written in the form in Eq. 30 of arXiv:2411.03302"
-   ]
-  },
   {
    "cell_type": "code",
    "execution_count": 14,
@@ -349,14 +341,15 @@
      "name": "stdout",
      "output_type": "stream",
      "text": [
-      "Z0 :  [45 46 48 49 51 52 54 55 57 58] , X0 :  [ 0  1  3  4  6  7  9 10 12 13]\n",
-      "Z1 :  [46 47 49 50 52 53 55 56 58 59] , X1 :  [ 1  2  4  5  7  8 10 11 13 14]\n",
-      "Z2 :  [60 61 63 64 66 67 69 70 72 73] , X2 :  [15 16 18 19 21 22 24 25 27 28]\n",
-      "Z3 :  [61 62 64 65 67 68 70 71 73 74] , X3 :  [16 17 19 20 22 23 25 26 28 29]\n",
-      "Z4 :  [ 0  3  6  9 12 15 18 21 24 27] , X4 :  [45 48 51 54 57 60 63 66 69 72]\n",
-      "Z5 :  [ 1  4  7 10 13 16 19 22 25 28] , X5 :  [46 49 52 55 58 61 64 67 70 73]\n",
-      "Z6 :  [15 18 21 24 27 30 33 36 39 42] , X6 :  [60 63 66 69 72 75 78 81 84 87]\n",
-      "Z7 :  [16 19 22 25 28 31 34 37 40 43] , X7 :  [61 64 67 70 73 76 79 82 85 88]\n"
+      "Z0 :  [ 2  5  8 11 14 17 20 23 26 29] , X0 :  [15 16 18 19 21 22 24 25 27 28]\n",
+      "Z1 :  [ 0  3  6  9 12 30 33 36 39 42] , X1 :  [15 17 18 20 21 23 24 26 27 29]\n",
+      "Z2 :  [ 1  4  7 10 13 31 34 37 40 43] , X2 :  [30 31 33 34 36 37 39 40 42 43]\n",
+      "Z3 :  [ 0  3  6  9 12 15 18 21 24 27] , X3 :  [30 32 33 35 36 38 39 41 42 44]\n",
+      "Z4 :  [ 1  2  3  4  5  7  8 12 13 14 15 18 19 20 24 31 32 33 36 37 38 39 40 41\n",
+      " 45 46] , X4 :  [15 17 20 22 24 25 32 37 39 42 45 60]\n",
+      "Z5 :  [ 1  3  7 15 16 18 21 24 25 30 31 33 34 39 45 47] , X5 :  [16 18 21 23 25 26 30 31 34 36 37 38 39 42 46 61]\n",
+      "Z6 :  [ 0  1  2  3  7  8 15 21 24 25 26 30 33 34 35 39 60 61] , X6 :  [ 0  5 10 12 15 17 20 22 24 25 30 35 40 42 45 75]\n",
+      "Z7 :  [ 3  4  6  7 12 16 18 22 30 31 33 36 39 40 60 62] , X7 :  [ 0  3  4  7  9 10 11 12 16 18 21 23 25 26 30 33 34 37 39 40 41 42 46 76]\n"
      ]
     }
    ],
@@ -400,7 +393,7 @@
      "name": "stderr",
      "output_type": "stream",
      "text": [
-      "100%|██████████| 100/100 [00:10<00:00,  9.45it/s]"
+      "100%|██████████| 100/100 [00:10<00:00,  9.50it/s]"
      ]
     },
     {
@@ -455,7 +448,7 @@
      "name": "stderr",
      "output_type": "stream",
      "text": [
-      "100%|██████████| 100/100 [00:09<00:00, 10.29it/s]\n"
+      "100%|██████████| 100/100 [00:09<00:00, 10.09it/s]\n"
      ]
     }
    ],

src/quits/qldpc_code.py

@@ -620,7 +620,8 @@ def __init__(self, p1, p2, lift_size, factor):
 
         self.hz = np.concatenate((h2, h1T), axis=1)
         self.hx = np.concatenate((h1, h2T), axis=1)
-        self.lz, self.lx = self.get_logicals()    # logical operators in the "canonical form"
+        self.lz, self.lx = compute_lz_and_lx(self.hx, self.hz)
+        # self.lz, self.lx = self.get_logicals()    # logical operators in the "canonical form"
 
     def get_block_mat(self, power):
         gen_mat = self.get_circulant_mat(self.factor, 1)
@@ -633,6 +634,7 @@ def get_block_mat(self, power):
         mat = mat * mat_placeholder * self.factor
         return mat, mat_placeholder
 
+    # WRONG; SHOULD BE FIXED LATER
     def get_logicals(self):
         '''
         :return: Logical operators of the code as a list of tuples (logical_z, logical_x)