randominputPauli.c

   1 #include <stdio.h>
   2 #include <stdlib.h>
   3 #include <time.h>
   4
   5 // order of matrix elements is [row][column]!!!
   6
   7 int main( int argc, char *argv[])
   8 {
   9
  10   if(argc != 3) {
  11     printf("randominputPauli arguments: \"number of qubits\" \"number of T gates\"\n");
  12     exit(0);
  13   }
  14
  15   int N = atoi(argv[1]);              // number of qubits
  16   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)
  17
  18   printf("%d\n", N);
  19   printf("%d\n", T);
  20
  21   int r;
  22   srand((unsigned)time(NULL));
  23
  24
  25   int i, j;
  26
  27   for(i=0; i<N; i++) {
  28
  29     r = rand()%4; printf("%d\n",r); // omega
  30
  31     for(j=0; j<N; j++) {
  32       r = rand()%4;
  33       printf("%d %d %d %d\n", r==0, r==1, r==2, r==3);
  34     }
  35   }
  36
  37 }