RIOT/tests/sys/puf_sram/tests/example_test.py

#!/usr/bin/env python3

# Copyright (c) 2018 Kevin Weiss, for HAW Hamburg  <kevin.weiss@haw-hamburg.de>
# Copyright (C) 2018 Peter Kietzmann, for HAW Hamburg <peter.kietzmann@haw-hamburg.de>
#
# This file is subject to the terms and conditions of the GNU Lesser
# General Public License v2.1. See the file LICENSE in the top level
# directory for more details.

import argparse
import puf_sram_if
import numpy

DEFAULT_POWER_CYCLES = 500
DEFAULT_OFF_TIME = 1
DEFAULT_BAUDRATE = 115200
DEFAULT_PORT = '/dev/ttyUSB0'
DEFAULT_INFO = True


def min_erntropy(all_meas):
    p1 = numpy.zeros(len(all_meas[0]))
    # number of ones for each bit
    for i in range(0, len(all_meas[0])):
        tmp = list(map(lambda x: int(x[i]), all_meas))
        p1[i] = numpy.count_nonzero(tmp)

    # probability of ones
    p1 = numpy.divide(p1, float(len(all_meas)))
    # probability of zeros
    p0 = 1 - p1

    p0_1_max = numpy.maximum(p1, p0)
    log2_p0_1_max = numpy.log2(p0_1_max)
    H_min = numpy.sum(-log2_p0_1_max)
    H_min_rel = 100 * H_min/len(p0_1_max)

    return [H_min, H_min_rel]


def main_func():
    p = argparse.ArgumentParser()
    p.add_argument("-n", "--number", type=int, default=DEFAULT_POWER_CYCLES,
                   help="Number of iterations, default: %s" % DEFAULT_POWER_CYCLES)
    p.add_argument("-t", "--off_time", type=int, default=DEFAULT_OFF_TIME,
                   help="Off time, default: %s [s]" % DEFAULT_OFF_TIME)
    p.add_argument("-p", "--port", type=str, default=DEFAULT_PORT,
                   help="Serial port, default: %s" % DEFAULT_PORT)
    p.add_argument("-b", "--baudrate", type=int, default=DEFAULT_BAUDRATE,
                   help="Baudrate of the serial port, default: %d" % DEFAULT_BAUDRATE)
    p.add_argument("-d", "--disable_output", default=DEFAULT_INFO, action='store_false',
                   help="Disable verbose output")
    args = p.parse_args()

    puf_sram = puf_sram_if.PufSram(port=args.port, baud=args.baudrate)
    seeds = puf_sram.get_seed_list(n=args.number, off_time=args.off_time,
                                   allow_print=args.disable_output)
    seeds = [format(x, '0>32b') for x in seeds]
    H_min, H_min_rel = min_erntropy(seeds)

    print("Number of seeds: %i       " % len(seeds))
    print("Seed length    : %i Bit   " % len(seeds[0]))
    print("Abs. Entropy   : %02.02f Bit   " % H_min)
    print("Rel. Entropy   : %02.02f perc. " % H_min_rel)


if __name__ == "__main__":
    main_func()