Created
October 5, 2012 21:27
-
-
Save osdf/3842524 to your computer and use it in GitHub Desktop.
Testing numpy and scipy setups
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#!/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 | |
x = numpy.random.random((1000,1000)) | |
setup = "import numpy; x = numpy.random.random((1000,1000))" | |
count = 5 | |
t = timeit.Timer("numpy.dot(x, x.T)", setup=setup) | |
print "dot:", t.timeit(count)/count, "sec" |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#!/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" |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Note that in recent numpy versions (>1.10) it will always report SLOW BLAS due to the missing
numpy.core._dotblas
.Here's a more updated way of checking if it's linked properly:
http://stackoverflow.com/questions/21671040/link-atlas-mkl-to-an-installed-numpy/21673585#21673585