randominputPauli.c

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