Fix cupy subnormals

Author: crusaderky

tests/test_funcs.py

@@ -235,6 +235,11 @@ def test_hypothesis(  # type: ignore[no-any-explicit,no-any-decorated]
             data.draw(npst.arrays(dtype=dtype, shape=(), elements=elements))
         )
         float_fill_value = float(fill_value)
+        if library is Backend.CUPY and dtype is np.float32:
+            # Avoid data-dependent dtype promotion when encountering subnormals
+            # close to the max float32 value
+            float_fill_value = float(np.clip(float_fill_value, -1e38, 1e38))
+
         arrays = tuple(
             xp.asarray(
                 data.draw(npst.arrays(dtype=dtype, shape=shape, elements=elements))