Add hyperbolic tests

Author: juntyr

quaddtype/release_tracker.md

@@ -50,12 +50,12 @@
 | arctan        | ✅    | ❌ _Need: basic tests + edge cases (NaN/inf/0/asymptotes)_           |
 | arctan2       | ✅    | ❌ _Need: basic tests + edge cases (NaN/inf/0/quadrant coverage)_    |
 | hypot         |       |                                                                      |
-| sinh          |       |                                                                      |
-| cosh          |       |                                                                      |
-| tanh          |       |                                                                      |
-| arcsinh       |       |                                                                      |
-| arccosh       |       |                                                                      |
-| arctanh       |       |                                                                      |
+| sinh          | ✅    | ✅                                                                   |
+| cosh          | ✅    | ✅                                                                   |
+| tanh          | ✅    | ✅                                                                   |
+| arcsinh       | ✅    | ✅                                                                   |
+| arccosh       | ✅    | ✅                                                                   |
+| arctanh       | ✅    | ✅                                                                   |
 | degrees       |       |                                                                      |
 | radians       |       |                                                                      |
 | deg2rad       |       |                                                                      |

quaddtype/tests/test_quaddtype.py

@@ -100,9 +100,9 @@ def test_array_minmax(op, a, b):
     np.testing.assert_array_equal(quad_res.astype(float), float_res)
 
     # Check sign for zero results
-    if float_result == 0.0:
-        assert np.signbit(float_result) == np.signbit(
-            quad_result), f"Zero sign mismatch for {op}({a}, {b})"
+    if float_res == 0.0:
+        assert np.signbit(float_res) == np.signbit(
+            quad_res), f"Zero sign mismatch for {op}({a}, {b})"
 
 
 @pytest.mark.parametrize("op", ["amin", "amax", "nanmin", "nanmax"])
@@ -119,9 +119,9 @@ def test_array_aminmax(op, a, b):
     np.testing.assert_array_equal(np.array(quad_res).astype(float), float_res)
 
     # Check sign for zero results
-    if float_result == 0.0:
-        assert np.signbit(float_result) == np.signbit(
-            quad_result), f"Zero sign mismatch for {op}({a}, {b})"
+    if float_res == 0.0:
+        assert np.signbit(float_res) == np.signbit(
+            quad_res), f"Zero sign mismatch for {op}({a}, {b})"
 
 
 @pytest.mark.parametrize("op", ["negative", "positive", "absolute", "sign", "signbit", "isfinite", "isinf", "isnan", "sqrt", "square", "reciprocal"])
@@ -478,3 +478,58 @@ def test_mod(a, b, backend, op):
         numpy_negative = numpy_result < 0
 
         assert result_negative == numpy_negative, f"Sign mismatch for {a} % {b}: quad={result_negative}, numpy={numpy_negative}"
+
+
+@pytest.mark.parametrize("op", ["sinh", "cosh", "tanh", "arcsinh", "arccosh", "arctanh"])
+@pytest.mark.parametrize("val", [
+    # Basic cases
+    "0.0", "-0.0", "1.0", "-1.0", "2.0", "-2.0",
+    # Small values
+    "1e-10", "-1e-10", "1e-15", "-1e-15",
+    # Values near one
+    "0.9", "-0.9", "0.9999", "-0.9999",
+    "1.1", "-1.1", "1.0001", "-1.0001",
+    # Medium values
+    "10.0", "-10.0", "20.0", "-20.0",
+    # Large values
+    "100.0", "200.0", "700.0", "1000.0",
+    "-100.0", "-200.0", "-700.0", "-1000.0",
+    # Fractional values
+    "0.5", "-0.5", "1.5", "-1.5", "2.5", "-2.5",
+    # Special values
+    "inf", "-inf", "nan", "-nan"
+])
+def test_hyperbolic_functions(op, val):
+    """Comprehensive test for hyperbolic functions: sinh, cosh, tanh, arcsinh, arccosh, arctanh"""
+    op_func = getattr(np, op)
+
+    quad_val = QuadPrecision(val)
+    float_val = float(val)
+
+    quad_result = op_func(quad_val)
+    float_result = op_func(float_val)
+
+    # Handle NaN cases
+    if np.isnan(float_result):
+        assert np.isnan(
+            float(quad_result)), f"Expected NaN for {op}({val}), got {float(quad_result)}"
+        return
+
+    # Handle infinity cases
+    if np.isinf(float_result):
+        assert np.isinf(
+            float(quad_result)), f"Expected inf for {op}({val}), got {float(quad_result)}"
+        assert np.sign(float_result) == np.sign(
+            float(quad_result)), f"Infinity sign mismatch for {op}({val})"
+        return
+
+    # For finite non-zero results
+    # Use relative tolerance for exponential functions due to their rapid growth
+    rtol = 1e-14 if abs(float_result) < 1e100 else 1e-10
+    np.testing.assert_allclose(float(quad_result), float_result, rtol=rtol, atol=1e-15,
+                               err_msg=f"Value mismatch for {op}({val})")
+
+    # Check sign for zero results
+    if float_result == 0.0:
+        assert np.signbit(float_result) == np.signbit(
+            quad_result), f"Zero sign mismatch for {op}({val})"