In [1]: import numpy as np

In [2]: import numexpr as ne

In [3]: a = np.arange(1e6)

In [4]: b = np.arange(1e6)

In [5]: ne.evaluate("a + 1")
Out[5]: 
array([  1.00000000e+00,   2.00000000e+00,   3.00000000e+00, ...,
         9.99998000e+05,   9.99999000e+05,   1.00000000e+06])

In [6]: ne.evaluate('a*b-4.1*a > 2.5*b')
Out[6]: array([False, False, False, ...,  True,  True,  True], dtype=bool)

In [7]: ne.evaluate("sin(a) + arcsinh(a/b)")
Out[7]: 
array([        nan,  1.72284457,  1.79067101, ...,  1.09567006,
        0.17523598, -0.09597844])