fixed phases of random Pauli input generators

[strong_simulation_stabilizer_rank.git] / randominputcommutingHermitianPauli2.c

diff --git a/randominputcommutingHermitianPauli2.c b/randominputcommutingHermitianPauli2.c
index 1c79fbb58acbbf70b6b86572c8aa3d4df7e69ee9..596d32759b0ad910830e237244c6246d41e30094 100644 (file)
--- a/randominputcommutingHermitianPauli2.c
+++ b/randominputcommutingHermitianPauli2.c
@@ -63,7 +63,7 @@ int main( int argc, char *argv[])
          Pauli[j][qubita] = 3; // X->Y
        else if(Pauli[j][qubita] == 3) {
          Pauli[j][qubita] = 2; // Y->-X
-         omega[j] = (2*omega[j])%2;
+         omega[j] = (2+omega[j])%4;
        }
       }
     } else { // CNOT
@@ -89,11 +89,11 @@ int main( int argc, char *argv[])
        } else if((Pauli[j][qubita] == 3) && (Pauli[j][qubitb] == 2)) {// YX ->-ZY
          Pauli[j][qubita] == 1;
          Pauli[j][qubitb] == 3;
-         omega[j] = -omega[j];
+         omega[j] = (2+omega[j])%4;
        } else if((Pauli[j][qubita] == 1) && (Pauli[j][qubitb] == 3)) {// ZY ->-YX
          Pauli[j][qubita] == 3;
          Pauli[j][qubitb] == 2;
-         omega[j] = -omega[j];
+         omega[j] = (2+omega[j])%4;
        }
       }
     }