From 515e5dd387898880a3477d5d02ddfac55ba2c5f1 Mon Sep 17 00:00:00 2001
From: Lucas K <lkocia@s3miclassical.com>
Date: Wed, 17 Mar 2021 16:15:51 -0700
Subject: [PATCH] fixed phases of random Pauli input generators

---
 randominputPauli.c                    | 8 ++++++--
 randominputcommutingHermitianPauli2.c | 6 +++---
 2 files changed, 9 insertions(+), 5 deletions(-)

diff --git a/randominputPauli.c b/randominputPauli.c
index da17abe..ee919d6 100644
--- a/randominputPauli.c
+++ b/randominputPauli.c
@@ -1,6 +1,7 @@
 #include <stdio.h>
 #include <stdlib.h>
 #include <time.h>
+#include <sys/auxv.h>
 
 // order of matrix elements is [row][column]!!!
 
@@ -19,14 +20,17 @@ int main( int argc, char *argv[])
   printf("%d\n", T);
 
   int r;
-  srand((unsigned)time(NULL));
+  unsigned int *seed;
+
+  seed = (unsigned int *)getauxval(AT_RANDOM);
+  srand(*seed);
 
 
   int i, j;
   
   for(i=0; i<N; i++) {
     
-    r = rand()%4; printf("%d\n",r); // omega
+    r = (rand()%2)*2; printf("%d\n",r); // omega
     
     for(j=0; j<N; j++) {
       r = rand()%4;
diff --git a/randominputcommutingHermitianPauli2.c b/randominputcommutingHermitianPauli2.c
index 1c79fbb..596d327 100644
--- 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;
 	}
       }
     }
-- 
2.30.2