The project is composed of these parts:

• Making Logistic random sequences
  • create_input.py
  • make_output.sh
  • logistic.c
  • logistic32*.c
  • draw_picture.py
  • draw_bit_picture.py
  • analyze_out.sh
  • average_period.py
  • compute_period.sh
  • compute_pass.sh
  • average_pass.py
  • integrate_period.sh
  • draw_period_picture.py
  • how to use
  • formats
    • *.in
    • *.out
    • *.a
    • *.times

Create input data such as 32-3.600000-1.in in the ./data directory. Save a file mu_list.txt with format "3.600000\n" for further analysis.

Use every *.in file as input for logistic32.c, create the corresponding output *.out. And the times between a char's appearance with last time are saved in *.times. And save a file outlist.txt in that directory for drawing program to read.

According to the input mu and x, output the next PRINT_TIMES numbers in [0, 255].

Input format: (precise, mu, x)

32 3.6 0.4

For every substitution result, output only in [0, 256) in binary char to the stdout. Output the distance of the same char output to the stderr.

Draw pictures of the .out files in ./data directory.

The .out file contains only random bytes.

Use every .out file as input bytes to dieharder analization. Each report is in the .a file.

Use lines just contain an integer as input. Output is a float of the average distance of elements.

Use every .times file as input to average_period.py to compute its average period. Store the float into .p.

Use every .a file as input to average_pass.py to count its pass times. The output is in report.txt.

Count their PASSED, FAILED, WEAK.

Integrate all .p files into periods.txt. Output format:

32
3.600000
255.35253111111112

Then call draw_period_picture.py to create a curve of average curve.

Use periods.txt to draw 4 lines for each precision.

bash clear.sh
python3 create_input.py
bash make_output.sh
python3 draw_bit_picture.py
bash analyze_out.sh
bash compute_period.sh
bash compute_pass.sh
bash integrate_period.sh
python3 draw_period_picture.py

* these below are just samples *
cat *.out | dieharder -o -f example.input -t 10
dieharder -g 202 -f example.input -a

dieharder -g 201 -f example.bin -a

8 3.6 0.1

(unreadable bytes)

#=============================================================================#
#            dieharder version 3.31.1 Copyright 2003 Robert G. Brown          #
#=============================================================================#
   rng_name    |           filename             |rands/second|
 file_input_raw|               08-3.600000-2.out|  4.37e+07  |
#=============================================================================#
        test_name   |ntup| tsamples |psamples|  p-value |Assessment
#=============================================================================#
   diehard_bitstream|   0|   2097152|     100|0.00000000|  FAILED  

148
172
204
93
206
207