Временная тема

[capzap]

unsigned int bit_count(unsigned int n) {
  unsigned int c;
  for (c=0;n>0;n=n&(n-1)) c++;
  return c;
}

unsigned int rol(unsigned int x, unsigned int n){
        return (x << n)|(x >> (32-n));
}

void rdm()
{
        int             i, counter;
        unsigned bit;
        int             num_0s, num_1s, v, bitsRead;
        unsigned char x = 42;
        unsigned char mask;
        unsigned int T_PLC_MS = 7*24*3600*1000;
        unsigned int tn;
        unsigned int tc;

        float m;
        float rx;
        float rdm;
        float last;
        float E = exp(1.0);

        if ( ((epsilon = (BitSequence *) calloc(tp.n, sizeof(BitSequence))) == NULL) ) {
                printf("Insufficient memory available.\n");
                exit(1);
        }
        counter = 1;
        last = 42;
        for ( v=0; v < tp.numOfBitStreams; v++ ) {
                mask = 0;
                num_0s = 0;
                num_1s = 0;
                bitsRead = 0;
                for ( i=0; i < tp.n; i++ ) {
//                      if (mask == 0) {
                                mask = 0x800000;
                                tn = clock() * 1000.0/ CLOCKS_PER_SEC;
                                tc = bit_count(tn);
                                tn |= (tn&(1 << (2))) << (31-2);
                                tn |= (tn&(1 << (5))) << (30-5);
                                tn |= (tn&(1 << (4))) << (29-4);
                                tn |= (tn&(1 << (1))) << (28-1);
                                tn |= (tn&(1 << (0))) << (27-0);
                                tn |= (tn&(1 << (7))) << (26-7);
                                tn |= (tn&(1 << (6))) << (25-6);
                                tn |= (tn&(1 << (3))) << (24-3);
                                tn = rol(tn, bit_count(tn)) | 0x80000001;
                                tn = tn % 71474513 + (tc + 77);
                                rdm = ((float)tn) / 10000000.0 * (E - (last < 0.0 ? 0 : (last > 1.0?1.0:last) ));
                                rdm = rdm - floor(rdm);
                                last = rdm;
//                              T_PLC_MS+=10;
//                      }
                        if ( rdm < 0.5 ) {//(x & mask) == 0 ) {
                                bit = 0;
                                num_0s++;
                        }
                        else {
                                bit = 1;
                                num_1s++;
                        }
                        mask /= 2;
                        bitsRead++;
                        epsilon[i] = bit;
                }
                fprintf(freqfp, "\t\tBITSREAD = %d 0s = %d 1s = %d\n", bitsRead, num_0s, num_1s); fflush(freqfp);
                nist_test_suite();
                }
        free(epsilon);
}

Что на этот раз?

Пересчет времени в секунды выполняется по формуле: clock() / CLOCKS_PER_SEC

следовательно Вы хотите получить миллисекунды, для работы в цикле, да еще на Вашем мощном компе у большинства итераций это будет постоянное число и это можно было бы не принимать в расчет, если бы следовать рекомендациям

in order to use rdm more then once within one sps cycle it need to be called with different seed numbers LAST

оно к сожалению у Вас тоже постоянное, а не значение предыдущей итерации как это делал я