Numba Jaccard

gistfile1.py

@numba.jit(target='cpu', nopython=True)
def _compute_overlaps(u, v):
    a = 0
    b = 0
    c = 0
    m = u.shape[0]
    for idx in range(m):
        a += u[idx] & v[idx]
        b += u[idx] & ~v[idx]
        c += ~u[idx] & v[idx]
    return a, b, c


@numba.jit(target='cpu', nopython=True)
def jaccard(u, v):
    a, b, c = _compute_overlaps(u, v)
    return 1 - (a / float(a + b + c))