for i in range(Ni):
    for j in range(Nj):
        for k in range(Nk):
            for l in range(Nl):
                fnew_v[j, k] += f[i, j, k, l] * b[i, l]
return fnew
Comparing Python accelerators (Cython,Numba,f2py) to Numpy einsum