Repository URL to install this package:
|
Version:
3.15.0 ▾
|
"""Make unit test for addmul() in mont.c"""
from common import counter, make_main, split64
def make_test(t, a, k):
if k == -1:
k = 0xFFFFFFFFFFFFFFFF
assert(0 <= k < 0x10000000000000000)
# What we expect the function to compute
result = t + a*k
# Turn a[] and t[] into arrays of 64-bit words
a = split64(a)
t_in = split64(t)
result = split64(result)
# Computation does not depend on zero terms
result_len = max(len(result), 1 + len(a))
# Pad the output vector with as many padding zeroes as needed
for x in range(result_len - len(t_in)):
t_in.append("0")
for x in range(result_len - len(result)):
result.append("0")
test_nr = next(counter)
print("")
print("void test_%d() {" % test_nr)
#print ' printf("Test #%d\\n");' % test_nr
print(" const uint64_t a[] = {" + ", ".join(a) + "};")
print(" uint64_t t[] = {" + ", ".join(t_in) + ", 0xAAAAAAAAAAAAAAAAULL};")
print(" const uint64_t expected_t[] = {" + ", ".join(result) + "};")
print("")
print(" addmul(t, %d, a, %d, 0x%x);" % (result_len, len(a), k))
print(" assert(memcmp(t, expected_t, 8*%d) == 0);" % result_len)
print(" assert(t[%d] == 0xAAAAAAAAAAAAAAAAULL);" % result_len)
print("}")
print("")
print("#include <assert.h>")
print("#include <string.h>")
print("#include <stdint.h>")
print("#include <stdio.h>")
print("")
print("void addmul(uint64_t *t, size_t tw, const uint64_t *a, size_t aw, uint64_t k);")
make_test(0, 0, 0)
make_test(0, 1, 1)
make_test(0, 5, 5)
make_test(0, 0xFFFFFFFFFFFFFFFFFFF, -1)
make_test(0xFFFFFFFFFFFFFFFF, 1, 1)
make_test(32783243204234329232323, 9232922323, 39393938)
make_test(32783243204234329232323333333333333783839393,
92329223233877777777777777777777777838333, 0x1000000)
make_test(37893272389423987423987429837498237498237498274982374982734982374982734982743982374,
30309093333334930430493049304930940394039430303000009090909093434930493094039409340930493094309403940394039403940394039403940390493049304943,
0x1000000)
make_main()