# test floating point floor divide and modulus
# it has some tricky corner cases

def test(x, y):
    div, mod = divmod(x, y)
    print('%.8f %.8f %.8f %.8f' % (x // y, x % y, div, mod))
    print(div == x // y, mod == x % y, abs(div * y + mod - x) < 1e-15)

test(1.23456, 0.7)
test(-1.23456, 0.7)
test(1.23456, -0.7)
test(-1.23456, -0.7)

a = 1.23456
b = 0.7
test(a, b)
test(a, -b)
test(-a, b)
test(-a, -b)

for i in range(25):
    x = (i - 12.5) / 6
    for j in range(25):
        y = (j - 12.5) / 6
        test(x, y)
