It would be useful if NumPy had a modular exponentiation function. Besides performance improvements over a**b%m, it also allows computations where a**b would not fit in the datatype. NumPy also seems ...