This problem is addressed by numexpr, which also handle multithreading.