Add python implementation

Author: antonwolfy

dpnp/dpnp_iface_mathematical.py

@@ -73,6 +73,7 @@
     DPNPRound,
     DPNPSinc,
     DPNPUnaryFunc,
+    DPNPUnaryTwoOutputsFunc,
     acceptance_fn_gcd_lcm,
     acceptance_fn_negative,
     acceptance_fn_positive,
@@ -112,6 +113,7 @@
     "fmax",
     "fmin",
     "fmod",
+    "frexp",
     "gcd",
     "gradient",
     "heaviside",
@@ -2342,6 +2344,77 @@ def ediff1d(ary, to_end=None, to_begin=None):
     mkl_impl_fn="_fmod",
 )
 
+
+_FREXP_DOCSTRING = """
+Decompose each element :math:`x_i` of the input array `x` into the mantissa and
+the twos exponent.
+
+For full documentation refer to :obj:`numpy.frexp`.
+
+Parameters
+----------
+x : {dpnp.ndarray, usm_ndarray}
+    Array of numbers to be decomposed, expected to have a real-valued
+    floating-point data type.
+out1 : {None, dpnp.ndarray, usm_ndarray}, optional
+    Output array for the mantissa to populate. Array must have the same shape
+    as `x` and the expected data type.
+
+    Default: ``None``.
+out2 : {None, dpnp.ndarray, usm_ndarray}, optional
+    Output array for the exponent to populate. Array must have the same shape
+    as `x` and the expected data type.
+
+    Default: ``None``.
+order : {None, "C", "F", "A", "K"}, optional
+    Memory layout of the newly output array, if parameter `out` is ``None``.
+
+    Default: ``"K"``.
+
+Returns
+-------
+mantissa : dpnp.ndarray
+    Floating values between -1 and 1.
+exponent : dpnp.ndarray
+    Integer exponents of 2.
+
+Limitations
+-----------
+Parameters `where` and `subok` are supported with their default values.
+Keyword argument `kwargs` is currently unsupported.
+Otherwise ``NotImplementedError`` exception will be raised.
+
+See Also
+--------
+:obj:`dpnp.ldexp` : Compute :math:`y = x1 * 2^{x2}`, inverse to
+                    :obj:`dpnp.frexp`.
+
+Notes
+-----
+Complex dtypes are not supported, they will raise a ``TypeError``.
+
+Examples
+--------
+>>> import dpnp as np
+>>> x = np.arange(9)
+>>> y1, y2 = np.frexp(x)
+>>> y1
+array([0.   , 0.5  , 0.5  , 0.75 , 0.5  , 0.625, 0.75 , 0.875, 0.5  ])
+>>> y2
+array([0, 1, 2, 2, 3, 3, 3, 3, 4], dtype=int32)
+>>> y1 * 2**y2
+array([0., 1., 2., 3., 4., 5., 6., 7., 8.])
+
+"""
+
+frexp = DPNPUnaryTwoOutputsFunc(
+    "_frexp",
+    ufi._frexp_result_type,
+    ufi._frexp,
+    _FREXP_DOCSTRING,
+)
+
+
 _GCD_DOCSTRING = r"""
 Returns the greatest common divisor of :math:`\abs{x1}` and :math:`\abs{x2}`.
 

pyproject.toml

@@ -154,7 +154,7 @@ max-returns = 8
 
 [tool.pylint.format]
 max-line-length = 80
-max-module-lines = 5000
+max-module-lines = 8000
 
 [tool.pylint.imports]
 allow-wildcard-with-all = true