Extend power operation

This is something that can be done with numpy.array but not with galois.FieldArray at the moment.

import numpy as np
import galois

2**np.arange(10)
array([  1,   2,   4,   8,  16,  32,  64, 128, 256, 512])

GF = galois.GF(2**9)
2**GF.Range(0, 10)
# TypeError: Operation 'power' requires the first operand to be a GF(2^9) array, not <class 'int'>.

GF(2)**GF.Range(0, 10)
# TypeError: Operation 'power' requires the second operand to be an integer array, not <class 'numpy.ndarray over GF(2^9)'>.

I think some kind of casting can be done here, similarly to numpy. What do you think?

[iyanmv]

Ups, you are right… I always was thinking in prime finite fields only, but the moment you try to do these things with extended fields, things get ugly.

[iyanmv]

I think I will close this because I think the only scenario where it would make a bit of sense is with prime finite fields. But then, choosing which field should be use to cast the arrays is completely arbitrary (first, highest order?).