'<float32 ** int> should return float32, not float64.'

def test_issue_2006(self):

    def foo(x, y):
        return np.power(x, y)
    pyfunc = foo
    cfunc = jit(nopython=True)(pyfunc)
    def check(x, y):
        got = cfunc(x, y)
        np.testing.assert_array_almost_equal(got, pyfunc(x, y))
        self.assertEqual(got.dtype, x.dtype)
    xs = [np.float32([1, 2, 3]), np.complex64([1j, 2, (3 - 3j)])]
    for x in xs:
        check(x, 3)
        check(x, np.uint64(3))
        check(x, np.int64([2, 2, 3]))
