numba-cuda 0.22.0__cp312-cp312-manylinux_2_24_aarch64.manylinux_2_28_aarch64.whl

numba_cuda/numba/cuda/tests/cudapy/test_math.py

@@ -0,0 +1,842 @@
+# SPDX-FileCopyrightText: Copyright (c) 2025 NVIDIA CORPORATION & AFFILIATES. All rights reserved.
+# SPDX-License-Identifier: BSD-2-Clause
+
+import sys
+import numpy as np
+import pytest
+from numba.cuda.testing import (
+    skip_unless_cc_53,
+    unittest,
+    CUDATestCase,
+    skip_on_cudasim,
+)
+from numba.cuda.np import numpy_support
+from numba import cuda, vectorize
+from numba.cuda import float32, float64, int32, void, int64
+import math
+
+
+def math_acos(A, B):
+    i = cuda.grid(1)
+    B[i] = math.acos(A[i])
+
+
+def math_asin(A, B):
+    i = cuda.grid(1)
+    B[i] = math.asin(A[i])
+
+
+def math_atan(A, B):
+    i = cuda.grid(1)
+    B[i] = math.atan(A[i])
+
+
+def math_acosh(A, B):
+    i = cuda.grid(1)
+    B[i] = math.acosh(A[i])
+
+
+def math_asinh(A, B):
+    i = cuda.grid(1)
+    B[i] = math.asinh(A[i])
+
+
+def math_atanh(A, B):
+    i = cuda.grid(1)
+    B[i] = math.atanh(A[i])
+
+
+def math_cos(A, B):
+    i = cuda.grid(1)
+    B[i] = math.cos(A[i])
+
+
+def math_sin(A, B):
+    i = cuda.grid(1)
+    B[i] = math.sin(A[i])
+
+
+def math_tan(A, B):
+    i = cuda.grid(1)
+    B[i] = math.tan(A[i])
+
+
+def math_cosh(A, B):
+    i = cuda.grid(1)
+    B[i] = math.cosh(A[i])
+
+
+def math_sinh(A, B):
+    i = cuda.grid(1)
+    B[i] = math.sinh(A[i])
+
+
+def math_tanh(A, B):
+    i = cuda.grid(1)
+    B[i] = math.tanh(A[i])
+
+
+def math_atan2(A, B, C):
+    i = cuda.grid(1)
+    C[i] = math.atan2(A[i], B[i])
+
+
+def math_exp(A, B):
+    i = cuda.grid(1)
+    B[i] = math.exp(A[i])
+
+
+def math_exp2(A, B):
+    i = cuda.grid(1)
+    B[i] = math.exp2(A[i])
+
+
+def math_erf(A, B):
+    i = cuda.grid(1)
+    B[i] = math.erf(A[i])
+
+
+def math_erfc(A, B):
+    i = cuda.grid(1)
+    B[i] = math.erfc(A[i])
+
+
+def math_expm1(A, B):
+    i = cuda.grid(1)
+    B[i] = math.expm1(A[i])
+
+
+def math_fabs(A, B):
+    i = cuda.grid(1)
+    B[i] = math.fabs(A[i])
+
+
+def math_gamma(A, B):
+    i = cuda.grid(1)
+    B[i] = math.gamma(A[i])
+
+
+def math_lgamma(A, B):
+    i = cuda.grid(1)
+    B[i] = math.lgamma(A[i])
+
+
+def math_log(A, B):
+    i = cuda.grid(1)
+    B[i] = math.log(A[i])
+
+
+def math_log2(A, B):
+    i = cuda.grid(1)
+    B[i] = math.log2(A[i])
+
+
+def math_log10(A, B):
+    i = cuda.grid(1)
+    B[i] = math.log10(A[i])
+
+
+def math_log1p(A, B):
+    i = cuda.grid(1)
+    B[i] = math.log1p(A[i])
+
+
+def math_remainder(A, B, C):
+    i = cuda.grid(1)
+    C[i] = math.remainder(A[i], B[i])
+
+
+def math_nextafter(A, B, C):
+    i = cuda.grid(1)
+    C[i] = math.nextafter(A[i], B[i])
+
+
+def math_sqrt(A, B):
+    i = cuda.grid(1)
+    B[i] = math.sqrt(A[i])
+
+
+def math_hypot(A, B, C):
+    i = cuda.grid(1)
+    C[i] = math.hypot(A[i], B[i])
+
+
+def math_pow(A, B, C):
+    i = cuda.grid(1)
+    C[i] = math.pow(A[i], B[i])
+
+
+def math_ceil(A, B):
+    i = cuda.grid(1)
+    B[i] = math.ceil(A[i])
+
+
+def math_floor(A, B):
+    i = cuda.grid(1)
+    B[i] = math.floor(A[i])
+
+
+def math_copysign(A, B, C):
+    i = cuda.grid(1)
+    C[i] = math.copysign(A[i], B[i])
+
+
+def math_fmod(A, B, C):
+    i = cuda.grid(1)
+    C[i] = math.fmod(A[i], B[i])
+
+
+def math_modf(A, B, C):
+    i = cuda.grid(1)
+    B[i], C[i] = math.modf(A[i])
+
+
+def math_isnan(A, B):
+    i = cuda.grid(1)
+    B[i] = math.isnan(A[i])
+
+
+def math_isinf(A, B):
+    i = cuda.grid(1)
+    B[i] = math.isinf(A[i])
+
+
+def math_isfinite(A, B):
+    i = cuda.grid(1)
+    B[i] = math.isfinite(A[i])
+
+
+def math_degrees(A, B):
+    i = cuda.grid(1)
+    B[i] = math.degrees(A[i])
+
+
+def math_radians(A, B):
+    i = cuda.grid(1)
+    B[i] = math.radians(A[i])
+
+
+def math_trunc(A, B):
+    i = cuda.grid(1)
+    B[i] = math.trunc(A[i])
+
+
+def math_pow_binop(A, B, C):
+    i = cuda.grid(1)
+    C[i] = A[i] ** B[i]
+
+
+def math_mod_binop(A, B, C):
+    i = cuda.grid(1)
+    C[i] = A[i] % B[i]
+
+
+class TestCudaMath(CUDATestCase):
+    def unary_template_float16(self, func, npfunc, start=0, stop=1):
+        self.unary_template(func, npfunc, np.float16, np.float16, start, stop)
+
+    def unary_template_float32(self, func, npfunc, start=0, stop=1):
+        self.unary_template(func, npfunc, np.float32, np.float32, start, stop)
+
+    def unary_template_float64(self, func, npfunc, start=0, stop=1):
+        self.unary_template(func, npfunc, np.float64, np.float64, start, stop)
+
+    def unary_template_int64(self, func, npfunc, start=0, stop=50):
+        self.unary_template(func, npfunc, np.int64, np.float64, start, stop)
+
+    def unary_template_uint64(self, func, npfunc, start=0, stop=50):
+        self.unary_template(func, npfunc, np.uint64, np.float64, start, stop)
+
+    def unary_template(self, func, npfunc, npdtype, nprestype, start, stop):
+        nelem = 50
+        A = np.linspace(start, stop, nelem, dtype=npdtype)
+        B = np.empty_like(A, dtype=nprestype)
+        arytype = numpy_support.from_dtype(npdtype)[::1]
+        restype = numpy_support.from_dtype(nprestype)[::1]
+        cfunc = cuda.jit((arytype, restype))(func)
+        cfunc[1, nelem](A, B)
+
+        # When this test was originally written it used
+        # assertTrue(np.allclose(...), which has different default tolerance
+        # values to assert_allclose. The tolerance values here are chosen as
+        # the tightest under which the tests will pass.
+        if npdtype == np.float64:
+            rtol = 1e-13
+        elif npdtype == np.float32:
+            rtol = 1e-6
+        else:
+            rtol = 1e-3
+        np.testing.assert_allclose(npfunc(A), B, rtol=rtol)
+
+    def unary_bool_special_values(self, func, npfunc, npdtype, npmtype):
+        fi = np.finfo(npdtype)
+        denorm = fi.tiny / 4
+        A = np.array(
+            [0.0, denorm, fi.tiny, 0.5, 1.0, fi.max, np.inf, np.nan],
+            dtype=npdtype,
+        )
+        B = np.empty_like(A, dtype=np.int32)
+        cfunc = cuda.jit((npmtype[::1], int32[::1]))(func)
+
+        cfunc[1, A.size](A, B)
+        np.testing.assert_array_equal(B, npfunc(A))
+
+        cfunc[1, A.size](-A, B)
+        np.testing.assert_array_equal(B, npfunc(-A))
+
+    def unary_bool_special_values_float32(self, func, npfunc):
+        self.unary_bool_special_values(func, npfunc, np.float32, float32)
+
+    def unary_bool_special_values_float64(self, func, npfunc):
+        self.unary_bool_special_values(func, npfunc, np.float64, float64)
+
+    def unary_bool_template_float32(self, func, npfunc, start=0, stop=1):
+        self.unary_template(func, npfunc, np.float32, np.float32, start, stop)
+
+    def unary_bool_template_float64(self, func, npfunc, start=0, stop=1):
+        self.unary_template(func, npfunc, np.float64, np.float64, start, stop)
+
+    def unary_bool_template_int32(self, func, npfunc, start=0, stop=49):
+        self.unary_template(func, npfunc, np.int32, np.int32, start, stop)
+
+    def unary_bool_template_int64(self, func, npfunc, start=0, stop=49):
+        self.unary_template(func, npfunc, np.int64, np.int64, start, stop)
+
+    def unary_bool_template(self, func, npfunc, npdtype, npmtype, start, stop):
+        nelem = 50
+        A = np.linspace(start, stop, nelem).astype(npdtype)
+        B = np.empty(A.shape, dtype=np.int32)
+        iarytype = npmtype[::1]
+        oarytype = int32[::1]
+        cfunc = cuda.jit((iarytype, oarytype))(func)
+        cfunc[1, nelem](A, B)
+        np.testing.assert_allclose(npfunc(A), B)
+
+    def binary_template_float32(self, func, npfunc, start=0, stop=1):
+        self.binary_template(func, npfunc, np.float32, np.float32, start, stop)
+
+    def binary_template_float64(self, func, npfunc, start=0, stop=1):
+        self.binary_template(func, npfunc, np.float64, np.float64, start, stop)
+
+    def binary_template_int64(self, func, npfunc, start=0, stop=50):
+        self.binary_template(func, npfunc, np.int64, np.float64, start, stop)
+
+    def binary_template_uint64(self, func, npfunc, start=0, stop=50):
+        self.binary_template(func, npfunc, np.uint64, np.float64, start, stop)
+
+    def binary_template(self, func, npfunc, npdtype, nprestype, start, stop):
+        nelem = 50
+        A = np.linspace(start, stop, nelem).astype(npdtype)
+        B = np.empty_like(A).astype(nprestype)
+        arytype = numpy_support.from_dtype(npdtype)[::1]
+        restype = numpy_support.from_dtype(nprestype)[::1]
+        cfunc = cuda.jit((arytype, arytype, restype))(func)
+        cfunc[1, nelem](A, A, B)
+        np.testing.assert_allclose(npfunc(A, A), B)
+
+    # ---------------------------------------------------------------------------
+    # test_math_acos
+
+    def test_math_acos(self):
+        self.unary_template_float32(math_acos, np.arccos)
+        self.unary_template_float64(math_acos, np.arccos)
+        # For integers we can only test with zero, since <=-1 and >=1 result in
+        # invalid values.
+        self.unary_template_int64(math_acos, np.arccos, start=0, stop=0)
+        self.unary_template_uint64(math_acos, np.arccos, start=0, stop=0)
+
+    # ---------------------------------------------------------------------------
+    # test_math_asin
+
+    def test_math_asin(self):
+        self.unary_template_float32(math_asin, np.arcsin)
+        self.unary_template_float64(math_asin, np.arcsin)
+        # For integers we can only test with zero, since <=-1 and >=1 result in
+        # invalid values.
+        self.unary_template_int64(math_asin, np.arcsin, start=0, stop=0)
+        self.unary_template_uint64(math_asin, np.arcsin, start=0, stop=0)
+
+    # ---------------------------------------------------------------------------
+    # test_math_atan
+
+    def test_math_atan(self):
+        self.unary_template_float32(math_atan, np.arctan)
+        self.unary_template_float64(math_atan, np.arctan)
+        self.unary_template_int64(math_atan, np.arctan)
+        self.unary_template_uint64(math_atan, np.arctan)
+
+    # ---------------------------------------------------------------------------
+    # test_math_acosh
+
+    def test_math_acosh(self):
+        self.unary_template_float32(math_acosh, np.arccosh, start=1, stop=2)
+        self.unary_template_float64(math_acosh, np.arccosh, start=1, stop=2)
+        self.unary_template_int64(math_acosh, np.arccosh, start=1, stop=2)
+        self.unary_template_uint64(math_acosh, np.arccosh, start=1, stop=2)
+
+    # ---------------------------------------------------------------------------
+    # test_math_asinh
+
+    def test_math_asinh(self):
+        self.unary_template_float32(math_asinh, np.arcsinh)
+        self.unary_template_float64(math_asinh, np.arcsinh)
+        self.unary_template_int64(math_asinh, np.arcsinh)
+        self.unary_template_uint64(math_asinh, np.arcsinh)
+
+    # ---------------------------------------------------------------------------
+    # test_math_atanh
+
+    def test_math_atanh(self):
+        self.unary_template_float32(math_atanh, np.arctanh, start=0, stop=0.9)
+        self.unary_template_float64(math_atanh, np.arctanh, start=0, stop=0.9)
+        self.unary_template_int64(math_atanh, np.arctanh, start=0, stop=0.9)
+        self.unary_template_uint64(math_atanh, np.arctanh, start=0, stop=0.9)
+
+    # ---------------------------------------------------------------------------
+    # test_math_cos
+
+    def test_math_cos(self):
+        self.unary_template_float32(math_cos, np.cos)
+        self.unary_template_float64(math_cos, np.cos)
+        self.unary_template_int64(math_cos, np.cos)
+        self.unary_template_uint64(math_cos, np.cos)
+
+    @skip_unless_cc_53
+    def test_math_fp16(self):
+        self.unary_template_float16(math_sin, np.sin)
+        self.unary_template_float16(math_cos, np.cos)
+        self.unary_template_float16(math_exp, np.exp)
+        self.unary_template_float16(math_log, np.log, start=1)
+        self.unary_template_float16(math_log2, np.log2, start=1)
+        self.unary_template_float16(math_log10, np.log10, start=1)
+        self.unary_template_float16(math_fabs, np.fabs, start=-1)
+        self.unary_template_float16(math_sqrt, np.sqrt)
+        self.unary_template_float16(math_ceil, np.ceil)
+        self.unary_template_float16(math_floor, np.floor)
+        if sys.version_info >= (3, 11):
+            self.unary_template_float16(math_exp2, np.exp2)
+
+    @skip_on_cudasim("numpy does not support trunc for float16")
+    @skip_unless_cc_53
+    def test_math_fp16_trunc(self):
+        self.unary_template_float16(math_trunc, np.trunc)
+
+    # ---------------------------------------------------------------------------
+    # test_math_sin
+
+    def test_math_sin(self):
+        self.unary_template_float32(math_sin, np.sin)
+        self.unary_template_float64(math_sin, np.sin)
+        self.unary_template_int64(math_sin, np.sin)
+        self.unary_template_uint64(math_sin, np.sin)
+
+    # ---------------------------------------------------------------------------
+    # test_math_tan
+
+    def test_math_tan(self):
+        self.unary_template_float32(math_tan, np.tan)
+        self.unary_template_float64(math_tan, np.tan)
+        self.unary_template_int64(math_tan, np.tan)
+        self.unary_template_uint64(math_tan, np.tan)
+
+    # ---------------------------------------------------------------------------
+    # test_math_cosh
+
+    def test_math_cosh(self):
+        self.unary_template_float32(math_cosh, np.cosh)
+        self.unary_template_float64(math_cosh, np.cosh)
+        self.unary_template_int64(math_cosh, np.cosh)
+        self.unary_template_uint64(math_cosh, np.cosh)
+
+    # ---------------------------------------------------------------------------
+    # test_math_sinh
+
+    def test_math_sinh(self):
+        self.unary_template_float32(math_sinh, np.sinh)
+        self.unary_template_float64(math_sinh, np.sinh)
+        self.unary_template_int64(math_sinh, np.sinh)
+        self.unary_template_uint64(math_sinh, np.sinh)
+
+    # ---------------------------------------------------------------------------
+    # test_math_tanh
+
+    def test_math_tanh(self):
+        self.unary_template_float32(math_tanh, np.tanh)
+        self.unary_template_float64(math_tanh, np.tanh)
+        self.unary_template_int64(math_tanh, np.tanh)
+        self.unary_template_uint64(math_tanh, np.tanh)
+
+    # ---------------------------------------------------------------------------
+    # test_math_atan2
+
+    def test_math_atan2(self):
+        self.binary_template_float32(math_atan2, np.arctan2)
+        self.binary_template_float64(math_atan2, np.arctan2)
+        self.binary_template_int64(math_atan2, np.arctan2)
+        self.binary_template_uint64(math_atan2, np.arctan2)
+
+    # ---------------------------------------------------------------------------
+    # test_math_erf
+
+    def test_math_erf(self):
+        @vectorize
+        def ufunc(x):
+            return math.erf(x)
+
+        self.unary_template_float32(math_erf, ufunc)
+        self.unary_template_float64(math_erf, ufunc)
+        self.unary_template_int64(math_erf, ufunc)
+        self.unary_template_uint64(math_erf, ufunc)
+
+    # ---------------------------------------------------------------------------
+    # test_math_erfc
+
+    def test_math_erfc(self):
+        @vectorize
+        def ufunc(x):
+            return math.erfc(x)
+
+        self.unary_template_float32(math_erfc, ufunc)
+        self.unary_template_float64(math_erfc, ufunc)
+        self.unary_template_int64(math_erfc, ufunc)
+        self.unary_template_uint64(math_erfc, ufunc)
+
+    # ---------------------------------------------------------------------------
+    # test_math_exp
+
+    def test_math_exp(self):
+        self.unary_template_float32(math_exp, np.exp)
+        self.unary_template_float64(math_exp, np.exp)
+        self.unary_template_int64(math_exp, np.exp)
+        self.unary_template_uint64(math_exp, np.exp)
+
+    # ---------------------------------------------------------------------------
+    # test_math_exp2
+
+    @unittest.skipUnless(sys.version_info >= (3, 11), "Python 3.11+ required")
+    def test_math_exp2(self):
+        self.unary_template_float32(math_exp2, np.exp2)
+        self.unary_template_float64(math_exp2, np.exp2)
+        self.unary_template_int64(math_exp2, np.exp2)
+        self.unary_template_uint64(math_exp2, np.exp2)
+
+    # ---------------------------------------------------------------------------
+    # test_math_expm1
+
+    def test_math_expm1(self):
+        self.unary_template_float32(math_expm1, np.expm1)
+        self.unary_template_float64(math_expm1, np.expm1)
+        self.unary_template_int64(math_expm1, np.expm1)
+        self.unary_template_uint64(math_expm1, np.expm1)
+
+    # ---------------------------------------------------------------------------
+    # test_math_fabs
+
+    def test_math_fabs(self):
+        self.unary_template_float32(math_fabs, np.fabs, start=-1)
+        self.unary_template_float64(math_fabs, np.fabs, start=-1)
+        self.unary_template_int64(math_fabs, np.fabs, start=-1)
+        self.unary_template_uint64(math_fabs, np.fabs, start=0)
+
+    # ---------------------------------------------------------------------------
+    # test_math_gamma
+
+    def test_math_gamma(self):
+        @vectorize
+        def ufunc(x):
+            return math.gamma(x)
+
+        self.unary_template_float32(math_gamma, ufunc, start=0.1)
+        self.unary_template_float64(math_gamma, ufunc, start=0.1)
+        self.unary_template_int64(math_gamma, ufunc, start=1)
+        self.unary_template_uint64(math_gamma, ufunc, start=1)
+
+    # ---------------------------------------------------------------------------
+    # test_math_lgamma
+
+    def test_math_lgamma(self):
+        @vectorize
+        def ufunc(x):
+            return math.lgamma(x)
+
+        self.unary_template_float32(math_lgamma, ufunc, start=0.1)
+        self.unary_template_float64(math_lgamma, ufunc, start=0.1)
+        self.unary_template_int64(math_lgamma, ufunc, start=1)
+        self.unary_template_uint64(math_lgamma, ufunc, start=1)
+
+    # ---------------------------------------------------------------------------
+    # test_math_log
+
+    def test_math_log(self):
+        self.unary_template_float32(math_log, np.log, start=1)
+        self.unary_template_float64(math_log, np.log, start=1)
+        self.unary_template_int64(math_log, np.log, start=1)
+        self.unary_template_uint64(math_log, np.log, start=1)
+
+    # ---------------------------------------------------------------------------
+    # test_math_log2
+
+    def test_math_log2(self):
+        self.unary_template_float32(math_log2, np.log2, start=1)
+        self.unary_template_float64(math_log2, np.log2, start=1)
+        self.unary_template_int64(math_log2, np.log2, start=1)
+        self.unary_template_uint64(math_log2, np.log2, start=1)
+
+    # ---------------------------------------------------------------------------
+    # test_math_log10
+
+    def test_math_log10(self):
+        self.unary_template_float32(math_log10, np.log10, start=1)
+        self.unary_template_float64(math_log10, np.log10, start=1)
+        self.unary_template_int64(math_log10, np.log10, start=1)
+        self.unary_template_uint64(math_log10, np.log10, start=1)
+
+    # ---------------------------------------------------------------------------
+    # test_math_log1p
+
+    def test_math_log1p(self):
+        self.unary_template_float32(math_log1p, np.log1p)
+        self.unary_template_float64(math_log1p, np.log1p)
+        self.unary_template_int64(math_log1p, np.log1p)
+        self.unary_template_uint64(math_log1p, np.log1p)
+
+    # ---------------------------------------------------------------------------
+    # test_math_remainder
+
+    def test_math_remainder(self):
+        self.binary_template_float32(math_remainder, np.remainder, start=1e-11)
+        self.binary_template_float64(math_remainder, np.remainder, start=1e-11)
+        self.binary_template_int64(math_remainder, np.remainder, start=1)
+        self.binary_template_uint64(math_remainder, np.remainder, start=1)
+
+    @skip_on_cudasim("math.remainder(0, 0) raises a ValueError on CUDASim")
+    def test_math_remainder_0_0(self):
+        @cuda.jit(void(float64[::1], int64, int64))
+        def test_0_0(r, x, y):
+            r[0] = math.remainder(x, y)
+
+        r = np.zeros(1, np.float64)
+        test_0_0[1, 1](r, 0, 0)
+        self.assertTrue(np.isnan(r[0]))
+
+    # ---------------------------------------------------------------------------
+    # test_math_nextafter
+
+    def test_math_nextafter(self):
+        self.binary_template_float32(math_nextafter, np.nextafter, start=1e-11)
+        self.binary_template_float64(math_remainder, np.remainder, start=1e-11)
+        self.binary_template_int64(math_remainder, np.remainder, start=1)
+        self.binary_template_uint64(math_remainder, np.remainder, start=1)
+
+    # ---------------------------------------------------------------------------
+    # test_math_sqrt
+
+    def test_math_sqrt(self):
+        self.unary_template_float32(math_sqrt, np.sqrt)
+        self.unary_template_float64(math_sqrt, np.sqrt)
+        self.unary_template_int64(math_sqrt, np.sqrt)
+        self.unary_template_uint64(math_sqrt, np.sqrt)
+
+    # ---------------------------------------------------------------------------
+    # test_math_hypot
+
+    def test_math_hypot(self):
+        self.binary_template_float32(math_hypot, np.hypot)
+        self.binary_template_float64(math_hypot, np.hypot)
+        self.binary_template_int64(math_hypot, np.hypot)
+        self.binary_template_uint64(math_hypot, np.hypot)
+
+    # ---------------------------------------------------------------------------
+    # test_math_pow
+
+    def pow_template_int32(self, npdtype):
+        nelem = 50
+        A = np.linspace(0, 25, nelem).astype(npdtype)
+        B = np.arange(nelem, dtype=np.int32)
+        C = np.empty_like(A)
+        arytype = numpy_support.from_dtype(npdtype)[::1]
+        cfunc = cuda.jit((arytype, int32[::1], arytype))(math_pow)
+        cfunc[1, nelem](A, B, C)
+
+        # NumPy casting rules result in a float64 output always, which doesn't
+        # match the overflow to inf of math.pow and libdevice.powi for large
+        # values of float32, so we compute the reference result with math.pow.
+        Cref = np.empty_like(A)
+        for i in range(len(A)):
+            Cref[i] = math.pow(A[i], B[i])
+
+        np.testing.assert_allclose(np.power(A, B).astype(npdtype), C, rtol=1e-6)
+
+    def test_math_pow(self):
+        self.binary_template_float32(math_pow, np.power)
+        self.binary_template_float64(math_pow, np.power)
+        with pytest.warns(RuntimeWarning, match="overflow encountered in cast"):
+            self.pow_template_int32(np.float32)
+        self.pow_template_int32(np.float64)
+
+    # ---------------------------------------------------------------------------
+    # test_math_pow_binop
+
+    def test_math_pow_binop(self):
+        self.binary_template_float32(math_pow_binop, np.power)
+        self.binary_template_float64(math_pow_binop, np.power)
+
+    # ---------------------------------------------------------------------------
+    # test_math_ceil
+
+    def test_math_ceil(self):
+        self.unary_template_float32(math_ceil, np.ceil)
+        self.unary_template_float64(math_ceil, np.ceil)
+        self.unary_template_int64(math_ceil, np.ceil)
+        self.unary_template_uint64(math_ceil, np.ceil)
+
+    # ---------------------------------------------------------------------------
+    # test_math_floor
+
+    def test_math_floor(self):
+        self.unary_template_float32(math_floor, np.floor)
+        self.unary_template_float64(math_floor, np.floor)
+        self.unary_template_int64(math_floor, np.floor)
+        self.unary_template_uint64(math_floor, np.floor)
+
+    # ---------------------------------------------------------------------------
+    # test_math_trunc
+    #
+    # Note that math.trunc() is only supported on NumPy float64s, and not
+    # other float types or int types. See NumPy Issue #13375:
+    #
+    # - https://github.com/numpy/numpy/issues/13375 - "Add methods from the
+    #   builtin float types to the numpy floating point types"
+
+    def test_math_trunc(self):
+        self.unary_template_float64(math_trunc, np.trunc)
+
+    @skip_on_cudasim("trunc only supported on NumPy float64")
+    def test_math_trunc_non_float64(self):
+        self.unary_template_float32(math_trunc, np.trunc)
+        self.unary_template_int64(math_trunc, np.trunc)
+        self.unary_template_uint64(math_trunc, np.trunc)
+
+    # ---------------------------------------------------------------------------
+    # test_math_copysign
+
+    def test_math_copysign(self):
+        self.binary_template_float32(math_copysign, np.copysign, start=-1)
+        self.binary_template_float64(math_copysign, np.copysign, start=-1)
+
+    # ---------------------------------------------------------------------------
+    # test_math_modf
+
+    def test_math_modf(self):
+        def modf_template_nan(dtype, arytype):
+            A = np.array([np.nan], dtype=dtype)
+            B = np.zeros_like(A)
+            C = np.zeros_like(A)
+            cfunc = cuda.jit((arytype, arytype, arytype))(math_modf)
+            cfunc[1, len(A)](A, B, C)
+            self.assertTrue(np.isnan(B))
+            self.assertTrue(np.isnan(C))
+
+        def modf_template_compare(A, dtype, arytype):
+            A = A.astype(dtype)
+            B = np.zeros_like(A)
+            C = np.zeros_like(A)
+            cfunc = cuda.jit((arytype, arytype, arytype))(math_modf)
+            cfunc[1, len(A)](A, B, C)
+            D, E = np.modf(A)
+            self.assertTrue(np.array_equal(B, D))
+            self.assertTrue(np.array_equal(C, E))
+
+        nelem = 50
+        # 32 bit float
+        with self.subTest("float32 modf on simple float"):
+            modf_template_compare(
+                np.linspace(0, 10, nelem), dtype=np.float32, arytype=float32[:]
+            )
+        with self.subTest("float32 modf on +- infinity"):
+            modf_template_compare(
+                np.array([np.inf, -np.inf]),
+                dtype=np.float32,
+                arytype=float32[:],
+            )
+        with self.subTest("float32 modf on nan"):
+            modf_template_nan(dtype=np.float32, arytype=float32[:])
+
+        # 64 bit float
+        with self.subTest("float64 modf on simple float"):
+            modf_template_compare(
+                np.linspace(0, 10, nelem), dtype=np.float64, arytype=float64[:]
+            )
+        with self.subTest("float64 modf on +- infinity"):
+            modf_template_compare(
+                np.array([np.inf, -np.inf]),
+                dtype=np.float64,
+                arytype=float64[:],
+            )
+        with self.subTest("float64 modf on nan"):
+            modf_template_nan(dtype=np.float64, arytype=float64[:])
+
+    # ---------------------------------------------------------------------------
+    # test_math_fmod
+
+    def test_math_fmod(self):
+        self.binary_template_float32(math_fmod, np.fmod, start=1)
+        self.binary_template_float64(math_fmod, np.fmod, start=1)
+
+    # ---------------------------------------------------------------------------
+    # test_math_mod_binop
+
+    def test_math_mod_binop(self):
+        self.binary_template_float32(math_mod_binop, np.fmod, start=1)
+        self.binary_template_float64(math_mod_binop, np.fmod, start=1)
+
+    # ---------------------------------------------------------------------------
+    # test_math_isnan
+
+    def test_math_isnan(self):
+        self.unary_bool_template_float32(math_isnan, np.isnan)
+        self.unary_bool_template_float64(math_isnan, np.isnan)
+        self.unary_bool_template_int32(math_isnan, np.isnan)
+        self.unary_bool_template_int64(math_isnan, np.isnan)
+        self.unary_bool_special_values_float32(math_isnan, np.isnan)
+        self.unary_bool_special_values_float64(math_isnan, np.isnan)
+
+    # ---------------------------------------------------------------------------
+    # test_math_isinf
+
+    def test_math_isinf(self):
+        self.unary_bool_template_float32(math_isinf, np.isinf)
+        self.unary_bool_template_float64(math_isinf, np.isinf)
+        self.unary_bool_template_int32(math_isinf, np.isinf)
+        self.unary_bool_template_int64(math_isinf, np.isinf)
+        self.unary_bool_special_values_float32(math_isinf, np.isinf)
+        self.unary_bool_special_values_float64(math_isinf, np.isinf)
+
+    # ---------------------------------------------------------------------------
+    # test_math_isfinite
+
+    def test_math_isfinite(self):
+        self.unary_bool_template_float32(math_isfinite, np.isfinite)
+        self.unary_bool_template_float64(math_isfinite, np.isfinite)
+        self.unary_bool_template_int32(math_isfinite, np.isfinite)
+        self.unary_bool_template_int64(math_isfinite, np.isfinite)
+        self.unary_bool_special_values_float32(math_isfinite, np.isfinite)
+        self.unary_bool_special_values_float64(math_isfinite, np.isfinite)
+
+    # ---------------------------------------------------------------------------
+    # test_math_degrees
+
+    def test_math_degrees(self):
+        self.unary_bool_template_float32(math_degrees, np.degrees)
+        self.unary_bool_template_float64(math_degrees, np.degrees)
+
+    # ---------------------------------------------------------------------------
+    # test_math_radians
+
+    def test_math_radians(self):
+        self.unary_bool_template_float32(math_radians, np.radians)
+        self.unary_bool_template_float64(math_radians, np.radians)
+
+
+if __name__ == "__main__":
+    unittest.main()