improve tests

Author: Aniketsy

Lib/test/test_complex.py

@@ -8,7 +8,7 @@
 )
 
 from random import random
-from math import isnan, copysign, ulp
+from math import isnan, copysign
 import operator
 
 INF = float("inf")
@@ -517,19 +517,13 @@ def test_pow_with_small_integer_exponents(self):
                         self.assertNotEqual(r1.imag, 0.0)
                     self.assertTrue(r2.real == 0.0 or r2.imag == 0.0)
 
-    @support.requires_IEEE_754
     def test_pow_small_negative_integer_exponents(self):
-        z = complex(float.fromhex('0x1.47e9c711723f5p+81'),
-                    float.fromhex('0x1.38afd1168e49fp+85'))
-        expected = complex(float.fromhex('0x0.4000000000000p-1022'),
-                           float.fromhex('0x0.3ffffffffffffp-1022'))
-        for exponent in (-12, -12.0, complex(-12.0, 0.0)):
-            with self.subTest(exponent=exponent):
-                result = z ** exponent
-                self.assertLessEqual(abs(result.real - expected.real),
-                                     4 * ulp(expected.real))
-                self.assertLessEqual(abs(result.imag - expected.imag),
-                                     4 * ulp(expected.imag))
+        for z, expected in [(complex(2**86, 0), complex(2**-1032, 0)),
+                            (complex(0, 2**86), complex(2**-1032, 0)),
+                            (complex(2**85, 2**85), complex(-2**-1026, 0))]:
+            for exponent in (-12, -12.0, complex(-12)):
+                with self.subTest(z=z, exponent=exponent):
+                    self.assertEqual(z**exponent, expected)
 
     def test_boolcontext(self):
         for i in range(100):

Objects/complexobject.c

@@ -368,23 +368,21 @@ c_powu(Py_complex x, long n)
 static Py_complex
 c_powi(Py_complex x, long n)
 {
-     if (n > 0)
-        return c_powu(x,n);
-    else if (n == 0)
-         return (Py_complex){1., 0.};
-
+    if (n > 0) {
+        return c_powu(x, n);
+    }
+    if (n == 0) {
+        return (Py_complex){1., 0.};
+    }
     double m = fabs(x.real) > fabs(x.imag) ? x.real : x.imag;
-
     if (m && isfinite(m)) {
         int e;
-
         frexp(m, &e);
-
         if (-e*n > 800) {
             x = (Py_complex){ldexp(x.real, -e), ldexp(x.imag, -e)};
             x = _Py_rc_quot(1.0, c_powu(x, -n));
-            x.real = ldexp(x.real, (int)(e * n));
-            x.imag = ldexp(x.imag, (int)(e * n));
+            x.real = ldexp(x.real, (int)(e*n));
+            x.imag = ldexp(x.imag, (int)(e*n));
             return x;
         }
     }