o

���g	e�
@sj
dZd
dgZ
ddlZddlmZ
ddlmZ
dd	lmZ
dd
lm	Z

ddlmZmZm
Z

ddlmZmZmZmZ
d
d�Zdd�ZGdd�d�Ze
je
je
je
je
je
jiZdZe
jee
j de�!d�
d�e
jee
j"de�!d�
d�e
j#ee
j$de�!d�
d�e
j%ee
j&de�!d�
d�iZ'iZ(dd�Z)iZ*dd �Z+d!d"�Z,d#d$�Z-ed%�
Gd&d
�d
��
Z.ed%�
Gd'd�d��
Z/dS)(zJMachine limits for Float32 and Float64 and (long double) if available...

�finfo�iinfo�N�)
�
set_module�
)
�MachAr)
�numeric)
�numerictypes)�array�inf�NaN)�log10�exp2�	nextafter�isnanc

C�|jd
kr|�
�}d|_|S)zfix rank-0 --> rank-1r)
r)�ndim�copy�shape�
�
a�r�?/usr/local/lib/python3.10/dist-packages/numpy/core/getlimits.py�_fr0�



rc

Cr)zfix rank > 0 --> rank-0rr)�sizerrrrrr�_fr1rrc@sNeZ
dZd
Zdd�
dd�Zedd��
Zedd	��
Zd
d�Zdd
�Z	dd�Z
dS)�
MachArLikez$ Object to simulate MachAr instance N)
�smallest_subnormalc	Ks4
t|
|_
|
|_|j
d
|_|s t|�d�
|�d�
|jd�|_n||_|�|�
|_|_|�|�
|_	|�|�
|_
|_|�|�
|_|�|�
|_
|_|j
d|�
|_|j�|�

tt|j�
�
|_|�|�d�
|j�
|_|�|j�
|_|�|j	�
|_|�|j�
|_|�|j
�
|_|�|j�
|_|�|j�
|_dS)N�titlerr�
�dtype�itype�
)�_MACHAR_PARAMS�params�ftyperr�_smallest_subnormal�_float_to_float�epsilon�eps�epsneg�xmax�huge�xmin�smallest_normal�tiny�ibeta�__dict__�update�intr
�	precision�_float_conv�
resolution�
_float_to_str�_str_eps�_str_epsneg�	_str_xmin�	_str_xmax�_str_resolution�_str_smallest_normal)	�selfr&r*r+r-r0r1r�kwargsrrr�__init__"s2





�









�




zMachArLike.__init__c
Cs6|j}
|�
d
�
|
krtjd�|j
�
tdd�
|�|
�
S)a 
Return the value for the smallest subnormal.

        Returns
        -------
        smallest_subnormal : float
            value for the smallest subnormal.

        Warns
        -----
        UserWarning
            If the calculated value for the smallest subnormal is zero.
        rz8The value of the smallest subnormal for {} type is zero.r�
�
stacklevel)r'r&�warnings�warn�format�UserWarningr(�r?�valuerrrr>s



�
zMachArLike.smallest_subnormalc

Cs|�|j
�
S)
z;Return the string representation of the smallest subnormal.)r8r�
r?rrr�_str_smallest_subnormalVsz"MachArLike._str_smallest_subnormalcCst|�
|
�
�
S)
z�Converts float to float.

        Parameters
        ----------
        value : float
            value to be converted.
        )rr6rHrrrr([�zMachArLike._float_to_floatcCst|
g
|j
�S)
z�Converts float to conv.

        Parameters
        ----------
        value : float
            value to be converted.
        )r
r&rHrrrr6erLzMachArLike._float_convcCs|jd
t
t|
�
d|j�S)zConverts float to str.

        Parameters
        ----------
        value : float
            value to be converted.
        �fmtr)r%r
rr&rHrrrr8oszMachArLike._float_to_str)�__name__�
__module__�__qualname__�__doc__rA�propertyrrKr(r6r8rrrrr s
�





rz(numpy {} precision floating point numberz%24.16e�double)r"rMrz%15.7e�singlez%szlong doublez%12.5e�halfcCs|t|
<dS�
N)
�_KNOWN_TYPES)�machar�bytepatrrr�_register_type�s
rZcCs�tj
}t|d
ddddddddt|d
�
�
t|d�
�
|d	�
|d
�
d�}
t|
d�
|
td<tj}t|d
ddddddddt|d
�
�
t|d�
�
|ddd�
t|d�
�
d�}t|d�
|td<tj}d}d}t|dddddddddd|d|||d �
|d�}t|d!�
|td"<tj}t|d#�
�
}	t|d$�
�
}
t	j
d%d&�
�
|d'�
|	|
|d �
}Wd�
n1s�w



Y
t|d(d#d$d)d*d+dddt|d(�
�
|	||
d�}t|d,�
|td<t|d-�
�
}
t|d$�
�
}t	j
d%d&�
�
|d'�
|
||d �
}Wd�
n1s�w



Y
t|d.d-d$d)d/d+dddt|d.�
�
|
||d�}t|d0�
|td1<t|t�
|d�
|d2�}t
}|td3d��
}t|d4d5ddd6ddddt|d4�
�
t|d5�
�
|||d7�}t|d8�
t|d9�
|td:<dS);Ni����i����i��r#�rri�g?)
�machep�negep�minexp�maxexp�it�iexpr1�irnd�ngrdr*r+r-r0sf�i��i��i�������g��?s��̽� g�<gi���i���i���i�4�g�<g�?�s���������@i����i����ignore)
�allri����i@�p�s����������������i����i�����?s
����������Pr gi����i�����i)r]r^r_r`rarbr1rcrdr*r+r-r0rs������Y<��������s��������������Y<�dd)�ntypes�float16rrrZ�	_float_ma�float32�float64�
longdoubler�errstaterrr)�f16�
float16_ma�f32�
float32_ma�f64�
epsneg_f64�tiny_f64�
float64_ma�ld�epsneg_f128�	tiny_f128�	huge_f128�float128_ma�
epsneg_f80�tiny_f80�huge_f80�
float80_ma�huge_dd�smallest_normal_dd�smallest_subnormal_dd�float_dd_marrr�_register_known_types�s
















�


















�

















�




�













�
�

�













�


















�
�
�r�c
Cs�t�
|�
}
|
d
urtt|�
�
�
|d�
|d�
�d�
��}d
}|tjkr,t�
|d
d��
}|d
ur5t�
|�
}|d
urQt	|�
dkrQdd�t�
�D�
}|�
|d
d��
}|d
u
rW|Stjd	|�d
|�d�t
dd
�
t|�
S)a� Get MachAr instance or MachAr-like instance

    Get parameters for floating point type, by first trying signatures of
    various known floating point types, then, if none match, attempting to
    identify parameters by analysis.

    Parameters
    ----------
    ftype : class
        Numpy floating point type class (e.g. ``np.float64``)

    Returns
    -------
    ma_like : instance of :class:`MachAr` or :class:`MachArLike`
        Object giving floating point parameters for `ftype`.

    Warns
    -----
    UserWarning
        If the binary signature of the float type is not in the dictionary of
        known float types.
    Ng�g$@�
<r#r[c
Ss*i|]\}
}t|
�
dkr|
d
d�|�qS)r[Nr#)
�len)�.0�
k�
vrrr�
<dictcomp>_
s*z_get_machar.<locals>.<dictcomp>z
Signature z for zz does not match any known type: falling back to type probe function.
This warnings indicates broken support for the dtype!rrB)r$�get�
ValueError�repr�newbyteorder�tobytesrurzrWr��itemsrDrErG�_discovered_machar)r&r%�key�ma_like�_ktrrr�_get_machar3
s&












�r�c

s>t��
t
�f
d
d��
f
dd��f
dd���
fdd��
d�S)z� Create MachAr instance with found information on float types

    TODO: MachAr should be retired completely ideally.  We currently only
          ever use it system with broken longdouble (valgrind, WSL).
    c

st|g
��SrV)
r
�
r��
r&rr�<lambda>s
sz$_discovered_machar.<locals>.<lambda>c

st|�
�d
�
�
dS)Nr"r)r�astyper�)
r%rrr�t
sc

stt
|�
d
��S)Nr�r
rr�r�rrr�u
sc

s�
d
tt
|�
d��S)NrMrr�r��r&r%rrr�v
sr)r$rr�rr�rr�l
s






�r��numpyc@sLeZ
dZd
ZiZdd�Zdd�Zdd�Zdd	�Ze	d
d��
Z
e	dd
��
ZdS)r
a

    finfo(dtype)

    Machine limits for floating point types.

    Attributes
    ----------
    bits : int
        The number of bits occupied by the type.
    dtype : dtype
        Returns the dtype for which `finfo` returns information. For complex
        input, the returned dtype is the associated ``float*`` dtype for its
        real and complex components.
    eps : float
        The difference between 1.0 and the next smallest representable float
        larger than 1.0. For example, for 64-bit binary floats in the IEEE-754
        standard, ``eps = 2**-52``, approximately 2.22e-16.
    epsneg : float
        The difference between 1.0 and the next smallest representable float
        less than 1.0. For example, for 64-bit binary floats in the IEEE-754
        standard, ``epsneg = 2**-53``, approximately 1.11e-16.
    iexp : int
        The number of bits in the exponent portion of the floating point
        representation.
    machep : int
        The exponent that yields `eps`.
    max : floating point number of the appropriate type
        The largest representable number.
    maxexp : int
        The smallest positive power of the base (2) that causes overflow.
    min : floating point number of the appropriate type
        The smallest representable number, typically ``-max``.
    minexp : int
        The most negative power of the base (2) consistent with there
        being no leading 0's in the mantissa.
    negep : int
        The exponent that yields `epsneg`.
    nexp : int
        The number of bits in the exponent including its sign and bias.
    nmant : int
        The number of bits in the mantissa.
    precision : int
        The approximate number of decimal digits to which this kind of
        float is precise.
    resolution : floating point number of the appropriate type
        The approximate decimal resolution of this type, i.e.,
        ``10**-precision``.
    tiny : float
        An alias for `smallest_normal`, kept for backwards compatibility.
    smallest_normal : float
        The smallest positive floating point number with 1 as leading bit in
        the mantissa following IEEE-754 (see Notes).
    smallest_subnormal : float
        The smallest positive floating point number with 0 as leading bit in
        the mantissa following IEEE-754.

    Parameters
    ----------
    dtype : float, dtype, or instance
        Kind of floating point or complex floating point
        data-type about which to get information.

    See Also
    --------
    iinfo : The equivalent for integer data types.
    spacing : The distance between a value and the nearest adjacent number
    nextafter : The next floating point value after x1 towards x2

    Notes
    -----
    For developers of NumPy: do not instantiate this at the module level.
    The initial calculation of these parameters is expensive and negatively
    impacts import times.  These objects are cached, so calling ``finfo()``
    repeatedly inside your functions is not a problem.

    Note that ``smallest_normal`` is not actually the smallest positive
    representable value in a NumPy floating point type. As in the IEEE-754
    standard [1]_, NumPy floating point types make use of subnormal numbers to
    fill the gap between 0 and ``smallest_normal``. However, subnormal numbers
    may have significantly reduced precision [2]_.

    This function can also be used for complex data types as well. If used,
    the output will be the same as the corresponding real float type
    (e.g. numpy.finfo(numpy.csingle) is the same as numpy.finfo(numpy.single)).
    However, the output is true for the real and imaginary components.

    References
    ----------
    .. [1] IEEE Standard for Floating-Point Arithmetic, IEEE Std 754-2008,
           pp.1-70, 2008, http://www.doi.org/10.1109/IEEESTD.2008.4610935
    .. [2] Wikipedia, "Denormal Numbers",
           https://en.wikipedia.org/wiki/Denormal_number

    Examples
    --------
    >>> np.finfo(np.float64).dtype
    dtype('float64')
    >>> np.finfo(np.complex64).dtype
    dtype('float32')

    cCsh
z|j�
|
�
}|du
r|WSWn	ty


Yn
w|
dur%tjd
tdd�
zt�|
�
}
Wnty<


t�t|
�
�
}
Yn
w|j�
|
�
}|du
rI|S|
g
}t�	|
�
}||
u
r\|�
|�

|}
t|
tj�sht
d|
�
�
|j�
|
�
}|du
rt|St|
tj�s�t|
}||
u
r�|�
|�

|}
|j�
|
d�}|du
r�|D]}||j|<q�|St�|�
�|
�
}|D]}||j|<q�|S)Nzifinfo() dtype cannot be None. This behavior will raise an error in the future. (Deprecated in NumPy 1.25)rrBzdata type %r not inexact)�_finfo_cacher��	TypeErrorrDrE�DeprecationWarningrr!�type�
obj2sctype�append�
issubclass�inexactr��floating�_convert_to_float�object�__new__�_init)�clsr!�obj�dtypes�newdtype�dtrrrr��
s\



�
�

�

�

























z
finfo.__new__cCs�t�
|
�
|_
t|
�
}d
D]}t||t||��
qdD]}t||t||�jd�
q|j
jd|_|jjd|_	|j	|_
|jjd|_|j|_
|j|_||_|j��|_|j��|_|j��|_|j��|_|j��|_|j��|_|j��|_|S)N)r5rbr`r_r^r])r7r+rrrg)rr!r��setattr�getattr�flat�itemsize�bitsr-�max�minr*rb�nexpra�nmant�_macharr;�strip�	_str_tinyr<�_str_maxr:r9r=r>rK)r?r!rX�wordrrrr�s*



















zfinfo._initc
Csd
}
|
|jS)Na+Machine parameters for %(dtype)s
---------------------------------------------------------------
precision = %(precision)3s   resolution = %(_str_resolution)s
machep = %(machep)6s   eps =        %(_str_eps)s
negep =  %(negep)6s   epsneg =     %(_str_epsneg)s
minexp = %(minexp)6s   tiny =       %(_str_tiny)s
maxexp = %(maxexp)6s   max =        %(_str_max)s
nexp =   %(nexp)6s   min =        -max
smallest_normal = %(_str_smallest_normal)s   smallest_subnormal = %(_str_smallest_subnormal)s
---------------------------------------------------------------
)
r2�r?rMrrr�__str__6s�

z
finfo.__str__c
Cs"|jj
}
|j��}|
|d
<d|S)N�klasszZ%(klass)s(resolution=%(resolution)s, min=-%(_str_max)s, max=%(_str_max)s, dtype=%(dtype)s))�	__class__rNr2r)r?�
c�
drrr�__repr__Fs





�zfinfo.__repr__c

Cs0t|j
jjd
�
rtjdtdd�
|j
jjd
S)a7
Return the value for the smallest normal.

        Returns
        -------
        smallest_normal : float
            Value for the smallest normal.

        Warns
        -----
        UserWarning
            If the calculated value for the smallest normal is requested for
            double-double.
        rz;The value of smallest normal is undefined for double doublerrB)rr�r/r�rDrErGrJrrrr/Ms


�zfinfo.smallest_normalc


Cs|jS)
aQ
Return the value for tiny, alias of smallest_normal.

        Returns
        -------
        tiny : float
            Value for the smallest normal, alias of smallest_normal.

        Warns
        -----
        UserWarning
            If the calculated value for the smallest normal is requested for
            double-double.
        )
r/rJrrrr0dsz
finfo.tinyN)rNrOrPrQr�r�r�r�r�rRr/r0rrrrr
z
sf8


c@sHeZ
dZd
ZiZiZdd�Zedd��
Zedd��
Z	dd	�Z
d
d�ZdS)
ra�
    iinfo(type)

    Machine limits for integer types.

    Attributes
    ----------
    bits : int
        The number of bits occupied by the type.
    dtype : dtype
        Returns the dtype for which `iinfo` returns information.
    min : int
        The smallest integer expressible by the type.
    max : int
        The largest integer expressible by the type.

    Parameters
    ----------
    int_type : integer type, dtype, or instance
        The kind of integer data type to get information about.

    See Also
    --------
    finfo : The equivalent for floating point data types.

    Examples
    --------
    With types:

    >>> ii16 = np.iinfo(np.int16)
    >>> ii16.min
    -32768
    >>> ii16.max
    32767
    >>> ii32 = np.iinfo(np.int32)
    >>> ii32.min
    -2147483648
    >>> ii32.max
    2147483647

    With instances:

    >>> ii32 = np.iinfo(np.int32(10))
    >>> ii32.min
    -2147483648
    >>> ii32.max
    2147483647

    cCs|zt�
|
�
|_
Wnty


t�
t|
�
�
|_
Yn
w|j
j|_|j
jd
|_d|j|jf|_|jdv
r<td|jf
�
�
dS)Nrgz%s%d�iuzInvalid integer data type %r.)	rr!r�r��kindr�r�r�r�)r?�int_typerrrrA�s



�





�ziinfo.__init__c
CsV|jd
krdSz	t
j|j}
W|
Sty*


td|jd>�
}
|
t
j|j<Y|
Sw)zMinimum value of given dtype.�
urr)r�r�	_min_valsr��KeyErrorr4r��r?�valrrrr��s


�


�z	iinfo.minc
Cshz	tj
|j}
W|
Sty3


|jd
krtd|j>d�
}
ntd|jd>d�
}
|
tj
|j<Y|
Sw)zMaximum value of given dtype.r�r)r�	_max_valsr�r�r�r4r�r�rrrr��s
�




�z	iinfo.maxc
Csd
}
|
|j|j
|jd�S)zString representation.z�Machine parameters for %(dtype)s
---------------------------------------------------------------
min = %(min)s
max = %(max)s
---------------------------------------------------------------
�r!r�r�r�r�rrrr��s�z
iinfo.__str__c

Csd
|jj
|j|j|jfS)Nz%s(min=%s, max=%s, dtype=%s))r�rNr�r�r!rJrrrr��s

�ziinfo.__repr__N)rNrOrPrQr�r�rArRr�r�r�r�rrrrrvs2




)0rQ�__all__rD�_utilsrr�r�rr	rur
rr�umathr
rrrrrr�csinglerT�complex_�float_�
clongfloat�	longfloatr��
_title_fmtrS�dict�int64rF�int32rz�longlongrU�int16r$rWrZrwr�r�r�r
rrrrr�<module>
sd




[

�


�


�


�


��
9
|