gistfile1.txt

@numba.jit(target='cpu', nopython=True)
def fast_simpson(i, j):
    a = 0
    b = 0
    c = 0
    
    M = i.shape[0]
    
    for idx in range(M):
        a += (i[idx] == 1) and (j[idx] == 1)
        b += (i[idx] == 1) and (j[idx] == 0)
        c += (i[idx] == 0) and (j[idx] == 1)
    
    return 1 - (a / float(min((a + b), (a +c))))