additional cleanup

Author: amicciche

src/LDPC_functions.jl

@@ -6,60 +6,6 @@ function stab_from_cxcz(Cx, Cz)
     return stab
 end
 
-# TODO delte this - these codes are not good at all
-# function getGoodLDPC(n=1)
-#     # Absolute paths to Cx and Cz npz files:
-#     if n==1
-#         Cx = npzread("/Users/micciche/Research/QuantumInfo23/JuliaProjects/codes_for_hardware_test/1_ra1_rb2_X_rankX120_rankZ179_minWtX2_minWtZ2.npz");
-#         Cz = npzread("/Users/micciche/Research/QuantumInfo23/JuliaProjects/codes_for_hardware_test/1_ra1_rb2_Z_rankX120_rankZ179_minWtX2_minWtZ2.npz");
-#         stab = stab_from_cxcz(Cx,Cz);
-#     elseif n==2
-#         Cx = npzread("/Users/micciche/Research/QuantumInfo23/JuliaProjects/codes_for_hardware_test/1_ra2_rb2_X_rankX226_rankZ120_minWtX2_minWtZ2.npz");
-#         Cz = npzread("/Users/micciche/Research/QuantumInfo23/JuliaProjects/codes_for_hardware_test/1_ra2_rb2_Z_rankX226_rankZ120_minWtX2_minWtZ2.npz");
-#         stab = stab_from_cxcz(Cx,Cz);
-#     elseif n==3
-#         Cx = npzread("/Users/micciche/Research/QuantumInfo23/JuliaProjects/codes_for_hardware_test/3_ra1_rb2_X_rankX120_rankZ179_minWtX2_minWtZ2.npz");
-#         Cz = npzread("/Users/micciche/Research/QuantumInfo23/JuliaProjects/codes_for_hardware_test/3_ra1_rb2_Z_rankX120_rankZ179_minWtX2_minWtZ2.npz");
-#         stab = stab_from_cxcz(Cx,Cz);
-#     end
-#     return stab, Cx, Cz
-# end
-
-function real_LDPC_numbers()
-    println("First one")
-    stab = getGoodLDPC(1);
-    data_qubits = size(stab)[2]
-    scirc, anc_bits, _ = naive_syndrome_circuit(stab)
-    total_qubits = anc_bits+data_qubits
-
-    println("Naive syndrome:", CircuitCompilation2xn.comp_numbers(scirc, total_qubits))
-    cat, scirc, anc_bits, _ = shor_syndrome_circuit(stab)
-    total_qubits = anc_bits+data_qubits
-    println("Shor syndrome:", CircuitCompilation2xn.comp_numbers(scirc, total_qubits))
-
-    println("\nSecond one")
-    stab = getGoodLDPC(2);
-    data_qubits = size(stab)[2]
-    scirc, anc_bits, _ = naive_syndrome_circuit(stab)
-    total_qubits = anc_bits+data_qubits
-
-    println("Naive syndrome:", CircuitCompilation2xn.comp_numbers(scirc, total_qubits))
-    cat, scirc, anc_bits, _ = shor_syndrome_circuit(stab)
-    total_qubits = anc_bits+data_qubits
-    println("Shor syndrome:", CircuitCompilation2xn.comp_numbers(scirc, total_qubits))
-
-    println("\nThird one")
-    stab = getGoodLDPC(3);
-    data_qubits = size(stab)[2]
-    scirc, anc_bits, _ = naive_syndrome_circuit(stab)
-    total_qubits = anc_bits+data_qubits
-
-    println("Naive syndrome:", CircuitCompilation2xn.comp_numbers(scirc, total_qubits))
-    cat, scirc, anc_bits, _ = shor_syndrome_circuit(stab)
-    total_qubits = anc_bits+data_qubits
-    println("Shor syndrome:", CircuitCompilation2xn.comp_numbers(scirc, total_qubits))
-end
-
 function generate_LDPC_matrix(n, k, w_r)
     function generate_row(n,w_r)
         onePart = ones(Bool, w_r)

src/plotting_CC2xn.jl

@@ -6,7 +6,7 @@ using QuantumClifford
 using QuantumClifford.ECC: Steane7, Shor9, naive_syndrome_circuit, shor_syndrome_circuit, parity_checks, code_s, code_n, code_k
 using QuantumClifford.ECC: naive_encoding_circuit, Cleve8, AbstractECC, Perfect5, CSS, PyBeliefPropOSDecoder
 using CairoMakie
-
+using LaTeXStrings
 using QuantumClifford.ECC: AbstractSyndromeDecoder, TableDecoder, evaluate_decoder, Toric, PyMatchingDecoder
 
 function my_plot_both_synd(code::AbstractECC, decoder::AbstractSyndromeDecoder, p_shift=0.0001, p_wait=1-exp(-14.5/28_000); name=string(typeof(code)))
@@ -269,7 +269,7 @@ function my_plot_both_synd(code::AbstractECC, decoder::AbstractSyndromeDecoder,
     lines!(f_shor_z, [0,0], [0,0], label="Shor-syndrome specialized comp", color=:orange)
 
     scatter!(f_shor_z, [0,0], [0,0], label="Constant near-term noise parameters \n(shift error, decoherence, gate noise)", color=:gray, marker=:utriangle)
-    scatter!(f_shor_z, [0,0], [0,0], label="Near-term noise parameters are \nmultiplied by 10*p_mem", color=:gray, marker=:circle)
+    scatter!(f_shor_z, [0,0], [0,0], label="Near-term noise parameters are \nmultiplied by 10  ×  p_mem", color=:gray, marker=:circle)
 
     xlims!(ax, high=-0.5, low=-4.0)
     ylims!(ax, high=-0.5, low=-4.0)
@@ -421,7 +421,7 @@ function plot_for_paper_figure(code::AbstractECC, decoder::AbstractSyndromeDecod
     lines!(f_shor_x, [0,0], [0,0], label="Shor-syndrome specialized comp", color=:orange)
 
     scatter!(f_shor_x, [0,0], [0,0], label="Constant near-term noise parameters \n(shift error, decoherence, gate noise)", color=:gray, marker=:utriangle)
-    scatter!(f_shor_x, [0,0], [0,0], label="Near-term noise parameters are \nmultiplied by 10*p_mem", color=:gray, marker=:circle)
+    scatter!(f_shor_x, [0,0], [0,0], label="Near-term noise parameters are \nmultiplied by 10 × p_mem", color=:gray, marker=:circle)
 
     xlims!(ax, high=-0.5, low=-4.0)
     ylims!(ax, high=-0.5, low=-4.0)
@@ -430,15 +430,15 @@ function plot_for_paper_figure(code::AbstractECC, decoder::AbstractSyndromeDecod
     return f 
 end
 
-f_Steane = my_plot_both_synd(Steane7(), TableDecoder(Steane7()))
+#f_Steane = my_plot_both_synd(Steane7(), TableDecoder(Steane7()))
 #f_Shor = my_plot_both_synd(Shor9(), TableDecoder(Shor9()))
 # f_Cleve = my_plot_both_synd(Cleve8(), TableDecoder(Cleve8()))
 # f_P5 = my_plot_both_synd(Perfect5(), TableDecoder(Perfect5()))
 
 # f_x_Gottesman3 = my_plot_both_synd(Gottesman(3), TableDecoder(Gottesman(3)))
 
-# f_t3 = my_plot_both_synd(Toric(3, 3), PyMatchingDecoder(Toric(3, 3)), name="Toric 3x3")
-f_t3 = plot_for_paper_figure(Toric(3, 3), PyMatchingDecoder(Toric(3, 3)), name="Toric 3x3")
+f_t3_both = my_plot_both_synd(Toric(3, 3), PyMatchingDecoder(Toric(3, 3)), name="Toric 3x3")
+#f_t3 = plot_for_paper_figure(Toric(3, 3), PyMatchingDecoder(Toric(3, 3)), name="Toric 3x3")
 
 #f_t4 = my_plot_both_synd(Toric(4, 4), PyMatchingDecoder(Toric(4, 4)), name="Toric 4x4")
 #f_t6 = my_plot_both_synd(Toric(6, 6), PyMatchingDecoder(Toric(6, 6)), name="Toric 6x6")