autojit_func = numba.autojit()(increment_counts_in_matrix_from_chain)
jit_func = numba.jit(argtypes=[numba.int64[:, ::1], numba.double[:, ::1]])(
    increment_counts_in_matrix_from_chain)
t = np.random.randint(0, 50, 500)
m1 = np.zeros((50, 50))
m2 = np.zeros((50, 50))
m3 = np.zeros((50, 50))
How can I speed up transition matrix creation in Numpy?