https://github.com/JuliaLang/julia/blob/master/test/perf/per...

я так и думал - от NumPy там только фактически import numpy да создание массивов. вот примерчик:

def mandel(z):
    n = 0
    c = z
    for n in xrange(0,79):
        if abs(z) > 2:
            n -= 1
            break
        z = z**2 + c
    return n + 1

def mandelperf():
    r1 = numpy.arange(-2.0, 0.5, 0.1)
    r2 = numpy.arange(-1.0, 1.0, 0.1)
    M = numpy.zeros((len(r1)*len(r2)))
    count = 0
    for r in r1:
        for i in r2:
            M[count] = mandel(complex(r,i))
            count += 1
    return M

Единственная тестовая функция, действительно задействующая NumPy для вычислений:

def randmatmul(n):
    A = matrix(numpy.random.rand(n,n))
    B = matrix(numpy.random.rand(n,n))
    return A*B

смотрим в таблицу - так и есть, для этого теста разница практически незаметная.