test_numpy_multithreading.py

#!/usr/bin/env python
import numpy
import sys
import timeit

try:
  import numpy.core._dotblas
  print 'FAST BLAS'
except ImportError:
  print 'slow blas'

print "version:", numpy.__version__
print "maxint:", sys.maxint
print

x = numpy.random.random((1000,1000)).astype(numpy.float32)

setup = "import numpy; x = numpy.random.random((1000,1000)).astype(numpy.float32)"
count = 5

t = timeit.Timer("numpy.dot(x, x.T)", setup=setup)
print "dot:", t.timeit(count)/count, "sec"

test_scipy_multithreading.py

#!/usr/bin/env python
import timeit

setup = "import numpy;\
         import scipy.linalg as linalg;\
         x = numpy.random.random((1000,1000));\
         z = numpy.dot(x, x.T)"
count = 5

t = timeit.Timer("linalg.cholesky(z, lower=True)", setup=setup)
print "cholesky:", t.timeit(count)/count, "sec"

t = timeit.Timer("linalg.svd(z)", setup=setup)
print "svd:", t.timeit(count)/count, "sec"