X-Git-Url: https://git.libre-soc.org/?a=blobdiff_plain;f=src%2Fadd%2Funit_test_single.py;h=2b0d9e569e2a86d3382ef17579fc62cf63fa835e;hb=6bff1a997f3846872cf489c24b5c01426c4dc97c;hp=ca33efb141193bdccb3315c05b22f67d6773b3fd;hpb=4388bc894114ed31ad1f3bf671a02ce7768c82e8;p=ieee754fpu.git diff --git a/src/add/unit_test_single.py b/src/add/unit_test_single.py index ca33efb1..2b0d9e56 100644 --- a/src/add/unit_test_single.py +++ b/src/add/unit_test_single.py @@ -1,3 +1,6 @@ +from random import randint +from random import seed + import sys from sfpy import Float32 @@ -7,6 +10,9 @@ def get_mantissa(x): def get_exponent(x): return ((x & 0x7f800000) >> 23) - 127 +def set_exponent(x, e): + return (x & ~0x7f800000) | ((e+127) << 23) + def get_sign(x): return ((x & 0x80000000) >> 31) @@ -30,51 +36,112 @@ def match(x, y): (x == y) ) -def get_case(dut, a, b): +def get_rs_case(dut, a, b, mid): + in_a, in_b = dut.rs[0] + out_z = dut.res[0] + yield dut.ids.in_mid.eq(mid) + yield in_a.v.eq(a) + yield in_a.valid_i.eq(1) + yield + yield + yield + yield + a_ack = (yield in_a.ready_o) + assert a_ack == 0 + + yield in_a.valid_i.eq(0) + + yield in_b.v.eq(b) + yield in_b.valid_i.eq(1) + yield + yield + b_ack = (yield in_b.ready_o) + assert b_ack == 0 + + yield in_b.valid_i.eq(0) + + yield out_z.ready_i.eq(1) + + while True: + out_z_stb = (yield out_z.valid_o) + if not out_z_stb: + yield + continue + vout_z = yield out_z.v + #out_mid = yield dut.ids.out_mid + yield out_z.ready_i.eq(0) + yield + break + + return vout_z, mid + +def check_rs_case(dut, a, b, z, mid=None): + if mid is None: + mid = randint(0, 6) + mid = 0 + out_z, out_mid = yield from get_rs_case(dut, a, b, mid) + assert out_z == z, "Output z 0x%x not equal to expected 0x%x" % (out_z, z) + assert out_mid == mid, "Output mid 0x%x != expected 0x%x" % (out_mid, mid) + + +def get_case(dut, a, b, mid): + #yield dut.in_mid.eq(mid) yield dut.in_a.v.eq(a) - yield dut.in_a.stb.eq(1) + yield dut.in_a.valid_i_test.eq(1) + yield + yield yield yield - a_ack = (yield dut.in_a.ack) + a_ack = (yield dut.in_a.ready_o) assert a_ack == 0 + + yield dut.in_a.valid_i.eq(0) + yield dut.in_b.v.eq(b) - yield dut.in_b.stb.eq(1) - b_ack = (yield dut.in_b.ack) + yield dut.in_b.valid_i.eq(1) + yield + yield + b_ack = (yield dut.in_b.ready_o) assert b_ack == 0 + yield dut.in_b.valid_i.eq(0) + + yield dut.out_z.ready_i.eq(1) + while True: - yield - out_z_stb = (yield dut.out_z.stb) + out_z_stb = (yield dut.out_z.valid_o) if not out_z_stb: + yield continue - yield dut.in_a.stb.eq(0) - yield dut.in_b.stb.eq(0) - yield dut.out_z.ack.eq(1) - yield - yield dut.out_z.ack.eq(0) - yield + out_z = yield dut.out_z.v + #out_mid = yield dut.out_mid + yield dut.out_z.ready_i.eq(0) yield break - out_z = yield dut.out_z.v - return out_z + return out_z, mid # TODO: mid -def check_case(dut, a, b, z): - out_z = yield from get_case(dut, a, b) +def check_case(dut, a, b, z, mid=None): + if mid is None: + mid = randint(0, 6) + mid = 0 + out_z, out_mid = yield from get_case(dut, a, b, mid) assert out_z == z, "Output z 0x%x not equal to expected 0x%x" % (out_z, z) + assert out_mid == mid, "Output mid 0x%x != expected 0x%x" % (out_mid, mid) -def run_test(dut, stimulus_a, stimulus_b): +def run_test(dut, stimulus_a, stimulus_b, op, get_case_fn): expected_responses = [] actual_responses = [] for a, b in zip(stimulus_a, stimulus_b): + mid = randint(0, 6) + mid = 0 af = Float32.from_bits(a) bf = Float32.from_bits(b) - z = af + bf - expected_responses.append(z.get_bits()) - #print (af, bf, z) - actual = yield from get_case(dut, a, b) + z = op(af, bf) + expected_responses.append((z.get_bits(), mid)) + actual = yield from get_case_fn(dut, a, b, mid) actual_responses.append(actual) if len(actual_responses) < len(expected_responses): @@ -83,10 +150,15 @@ def run_test(dut, stimulus_a, stimulus_b): for expected, actual, a, b in zip(expected_responses, actual_responses, stimulus_a, stimulus_b): - passed = match(expected, actual) + passed = match(expected[0], actual[0]) + if expected[1] != actual[1]: # check mid + print ("MID failed", expected[1], actual[1]) + sys.exit(0) if not passed: + expected = expected[0] + actual = actual[0] print ("Fail ... expected:", hex(expected), "actual:", hex(actual)) print (hex(a)) @@ -111,3 +183,73 @@ def run_test(dut, stimulus_a, stimulus_b): sys.exit(0) +corner_cases = [0x80000000, 0x00000000, 0x7f800000, 0xff800000, + 0x7fc00000, 0xffc00000] + +def run_corner_cases(dut, count, op, get_case_fn): + #corner cases + from itertools import permutations + stimulus_a = [i[0] for i in permutations(corner_cases, 2)] + stimulus_b = [i[1] for i in permutations(corner_cases, 2)] + yield from run_test(dut, stimulus_a, stimulus_b, op, get_case_fn) + count += len(stimulus_a) + print (count, "vectors passed") + +def run_test_2(dut, stimulus_a, stimulus_b, op, get_case_fn): + yield from run_test(dut, stimulus_a, stimulus_b, op, get_case_fn) + yield from run_test(dut, stimulus_b, stimulus_a, op, get_case_fn) + +def run_cases(dut, count, op, fixed_num, num_entries, get_case_fn): + if isinstance(fixed_num, int): + stimulus_a = [fixed_num for i in range(num_entries)] + report = hex(fixed_num) + else: + stimulus_a = fixed_num + report = "random" + + stimulus_b = [randint(0, 1<<32) for i in range(num_entries)] + yield from run_test_2(dut, stimulus_a, stimulus_b, op, get_case_fn) + count += len(stimulus_a) + print (count, "vectors passed 2^32", report) + + # non-canonical NaNs. + stimulus_b = [set_exponent(randint(0, 1<<32), 128) \ + for i in range(num_entries)] + yield from run_test_2(dut, stimulus_a, stimulus_b, op, get_case_fn) + count += len(stimulus_a) + print (count, "vectors passed Non-Canonical NaN", report) + + # -127 + stimulus_b = [set_exponent(randint(0, 1<<32), -127) \ + for i in range(num_entries)] + yield from run_test_2(dut, stimulus_a, stimulus_b, op, get_case_fn) + count += len(stimulus_a) + print (count, "vectors passed exp=-127", report) + + # nearly zero + stimulus_b = [set_exponent(randint(0, 1<<32), -126) \ + for i in range(num_entries)] + yield from run_test_2(dut, stimulus_a, stimulus_b, op, get_case_fn) + count += len(stimulus_a) + print (count, "vectors passed exp=-126", report) + + # nearly inf + stimulus_b = [set_exponent(randint(0, 1<<32), 127) \ + for i in range(num_entries)] + yield from run_test_2(dut, stimulus_a, stimulus_b, op, get_case_fn) + count += len(stimulus_a) + print (count, "vectors passed exp=127", report) + + return count + +def run_edge_cases(dut, count, op, get_case_fn): + #edge cases + for testme in corner_cases: + count = yield from run_cases(dut, count, op, testme, 10, get_case_fn) + + for i in range(100000): + stimulus_a = [randint(0, 1<<32) for i in range(10)] + count = yield from run_cases(dut, count, op, stimulus_a, 10, + get_case_fn) + return count +