X-Git-Url: https://s3miclassical.com/gitweb/?p=strong_simulation_gauss_sum_rank.git;a=blobdiff_plain;f=randommultipleinputPaulis.c;fp=randommultipleinputPaulis.c;h=d49feec0fae6b7cb2a1e94bd90a2fc0204f4ca80;hp=0000000000000000000000000000000000000000;hb=41acfb8cfb426c10f7a16e8cb04337012b25e3fb;hpb=13d8b4a4c252c341d36d3d90f66ed28600777224

diff --git a/randommultipleinputPaulis.c b/randommultipleinputPaulis.c
new file mode 100644
index 0000000..d49feec
--- /dev/null
+++ b/randommultipleinputPaulis.c
@@ -0,0 +1,38 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <time.h>
+
+// order of matrix elements is [row][column]!!!
+
+int main( int argc, char *argv[])
+{
+
+  if(argc != 4) {
+    printf("randommultipleinputPaulis arguments: \"number of qubits\" \"number of T gates\" \"number of Paulis\"\n");
+    exit(0);
+  }
+  
+  int N = atoi(argv[1]);              // number of qubits
+  int T = atoi(argv[2]);              // number of T gate magic states (set to the first 'K' of the 'N' qubits -- the rest are set to the '0' computational basis state)
+  int numPaulis = atoi(argv[3]);
+
+  printf("%d\n", N);
+  //printf("%d\n", T);
+
+  int r;
+  srand((unsigned)time(NULL));
+
+
+  int i, j;
+  
+  for(i=0; i<numPaulis; i++) {
+    
+    //r = rand()%4; printf("%d\n",r); // omega
+    
+    for(j=0; j<N; j++) {
+      r = rand()%4;
+      printf("%d %d %d %d\n", r==0, r==1, r==2, r==3);
+    }
+  }
+
+}