def compare_implementations(size, coverage=1):
    NN = size * coverage
    ar = generate_positions(size, coverage)
    res = CS.computeIndices(ar[:, 0], ar[:, 1])
    resDP = CS.computeIndicesDP(ar[:, 0], ar[:, 1])
    assert (sum(res[:, 0:2] == resDP[:, 0:2]) == NN * 2)
def compare_implementations(size, coverage = 1):
    NN = size * coverage
    ar = generate_positions(size, coverage)
    res = CS.computeIndices(ar[:,0], ar[:,1])
    resDP = CS.computeIndicesDP(ar[:,0], ar[:,1])
    assert(sum(res[:,0:2] == resDP[:,0:2]) == NN*2)