Merge branch 'dev' into fixdev

Author: Fe-r-oz

src/Classical/weight_dist.jl

@@ -1019,6 +1019,17 @@ function minimum_words(C::AbstractLinearCode)
     end
 end
 
+
+# classical version of QDistRnd here
+
+
+
+
+
+
+
+
+
 #############################
     # Weight Enumerators
 #############################

src/CodingTheory.jl

@@ -362,7 +362,7 @@ include("Quantum/weight_dist.jl")
 # export weight_plot_CSS_X, weight_plot_CSS_Z, weight_plot_CSS, minimum_distance_X_Z,
 #     minimum_distance_X, minimum_distance_Z, is_pure, QDistRndCSS
 export minimum_distance_upper_bound!, random_information_set_minimum_distance_bound!,
-    QDistRand!
+    QDistRnd!
 
 #############################
 #  Quantum/product_codes.jl

src/Quantum/subsystem_code.jl

@@ -282,8 +282,8 @@ function SubsystemCode(S::CTMatrixTypes, L::CTMatrixTypes, G::CTMatrixTypes;
     args = _is_CSS_symplectic(S, signs, true)
     if args[1]
         # bare
-        X_logs = reduce(vcat, [log[1][:, 1:n] for log in logs])
-        Z_logs = reduce(vcat, [log[2][:, n + 1:end] for log in logs])
+        X_logs = reduce(vcat, [log[1][:, 1:n] for log in log_pairs])
+        Z_logs = reduce(vcat, [log[2][:, n + 1:end] for log in log_pairs])
         _, X_mat = rref(vcat(args[2], X_logs))
         anti = _remove_empty(X_mat, :rows) * transpose(Z_logs)
         u_bound_dx_bare, _ = _min_wt_row(X_mat[findall(!iszero(anti[i:i, :]) for i in axes(anti, 1)), :])
@@ -292,10 +292,10 @@ function SubsystemCode(S::CTMatrixTypes, L::CTMatrixTypes, G::CTMatrixTypes;
         u_bound_dz_bare, _ = _min_wt_row(Z_mat[findall(!iszero(anti[i:i, :]) for i in axes(anti, 1)), :])
 
         # dressed
-        _, X_mat = rref(vcat(X_mat, reduce(vcat, [log[1][:, 1:n] for log in gauge_ops])))
+        _, X_mat = rref(vcat(X_mat, reduce(vcat, [log[1][:, 1:n] for log in g_ops_pairs])))
         anti = _remove_empty(X_mat, :rows) * transpose(Z_logs)
         u_bound_dx_dressed, _ = _min_wt_row(X_mat[findall(!iszero(anti[i:i, :]) for i in axes(anti, 1)), :])
-        _, Z_mat = rref(vcat(Z_mat, reduce(vcat, [log[2][:, n + 1:end] for log in gauge_ops])))
+        _, Z_mat = rref(vcat(Z_mat, reduce(vcat, [log[2][:, n + 1:end] for log in g_ops_pairs])))
         anti = _remove_empty(Z_mat, :rows) * transpose(X_logs)
         u_bound_dz_dressed, _ =_min_wt_row(Z_mat[findall(!iszero(anti[i:i, :]) for i in axes(anti, 1)), :])
         return SubsystemCodeCSS(F, n, k, r, missing, missing, missing, missing, missing,