Ë
|
Kñúh. ã ó` d Z dgZddlmZ ddlmZ G d d« Zedk( r e e« « yy)z
|
Machine arithmetic - determine the parameters of the
|
floating-point arithmetic system
|
|
Author: Pearu Peterson, September 2003
|
|
ÚMachAré )Úerrstate)Úanyc ó0 e Zd ZdZeeed dfdZd Zd Zy)r aw
|
Diagnosing machine parameters.
|
|
Attributes
|
----------
|
ibeta : int
|
Radix in which numbers are represented.
|
it : int
|
Number of base-`ibeta` digits in the floating point mantissa M.
|
machep : int
|
Exponent of the smallest (most negative) power of `ibeta` that,
|
added to 1.0, gives something different from 1.0
|
eps : float
|
Floating-point number ``beta**machep`` (floating point precision)
|
negep : int
|
Exponent of the smallest power of `ibeta` that, subtracted
|
from 1.0, gives something different from 1.0.
|
epsneg : float
|
Floating-point number ``beta**negep``.
|
iexp : int
|
Number of bits in the exponent (including its sign and bias).
|
minexp : int
|
Smallest (most negative) power of `ibeta` consistent with there
|
being no leading zeros in the mantissa.
|
xmin : float
|
Floating-point number ``beta**minexp`` (the smallest [in
|
magnitude] positive floating point number with full precision).
|
maxexp : int
|
Smallest (positive) power of `ibeta` that causes overflow.
|
xmax : float
|
``(1-epsneg) * beta**maxexp`` (the largest [in magnitude]
|
usable floating value).
|
irnd : int
|
In ``range(6)``, information on what kind of rounding is done
|
in addition, and on how underflow is handled.
|
ngrd : int
|
Number of 'guard digits' used when truncating the product
|
of two mantissas to fit the representation.
|
epsilon : float
|
Same as `eps`.
|
tiny : float
|
An alias for `smallest_normal`, kept for backwards compatibility.
|
huge : float
|
Same as `xmax`.
|
precision : float
|
``- int(-log10(eps))``
|
resolution : float
|
``- 10**(-precision)``
|
smallest_normal : float
|
The smallest positive floating point number with 1 as leading bit in
|
the mantissa following IEEE-754. Same as `xmin`.
|
smallest_subnormal : float
|
The smallest positive floating point number with 0 as leading bit in
|
the mantissa following IEEE-754.
|
|
Parameters
|
----------
|
float_conv : function, optional
|
Function that converts an integer or integer array to a float
|
or float array. Default is `float`.
|
int_conv : function, optional
|
Function that converts a float or float array to an integer or
|
integer array. Default is `int`.
|
float_to_float : function, optional
|
Function that converts a float array to float. Default is `float`.
|
Note that this does not seem to do anything useful in the current
|
implementation.
|
float_to_str : function, optional
|
Function that converts a single float to a string. Default is
|
``lambda v:'%24.16e' %v``.
|
title : str, optional
|
Title that is printed in the string representation of `MachAr`.
|
|
See Also
|
--------
|
finfo : Machine limits for floating point types.
|
iinfo : Machine limits for integer types.
|
|
References
|
----------
|
.. [1] Press, Teukolsky, Vetterling and Flannery,
|
"Numerical Recipes in C++," 2nd ed,
|
Cambridge University Press, 2002, p. 31.
|
|
c ó
|
| dS )Nz24.16e© )Úvs úFH:\Change_password\venv_build\Lib\site-packages\numpy/_core/_machar.pyú<lambda>zMachAr.<lambda>i s
|
¨1¨V¨* ó z Python floating point numberc óp t d¬« 5 | j |||||« ddd« y# 1 sw Y yxY w)a!
|
|
float_conv - convert integer to float (array)
|
int_conv - convert float (array) to integer
|
float_to_float - convert float array to float
|
float_to_str - convert array float to str
|
title - description of used floating point numbers
|
|
Úignore)ÚunderN)r Ú_do_init)ÚselfÚ
|
float_convÚint_convÚfloat_to_floatÚ float_to_strÚtitles r
|
Ú__init__zMachAr.__init__g s9 ô HÔ %ñ UØ MM* h°À ÈeÔ T÷ U÷ Uñ Uús ,¬5c ó² d}d} |d« }||z } ||z
|
}
|
|} t |« D ]$ } | | z } | |z } | | z
|
}t ||z
|
|
|
k7 « s$ n t | |j fz « |}t |« D ]' } ||z }| |z } || | z
|
« }t |dk7 « s' n t || |j fz « |} ||« }d}|}t |« D ]) } |dz }||z }||z } | |z
|
}t ||z
|
|
|
k7 « s) n t || |j fz « || z }|} t |« D ]$ } | | z } | |z } | | z
|
}t ||z
|
|
|
k7 « s$ n t || |j fz « | |z } d}t | | z
|
|
|
k7 « rd}| |z }||z } |dk( rt | |z
|
|
|
k7 « rd}|dz }||z }|} t |« D ] }| |z } | }t |« D ]? } || z
|
} t | |z
|
|
|
k7 « r nA| |z } |dz
|
}|dk s+t dt « z « t || |j fz « | }| }| dz
|
}|} t |« D ]$ } || z } t | |z
|
|
|
k7 « r n&| |z } |dz }& t || |j fz « | } d} || z } |dk( rt | |z |z
|
|
|
k7 « rd} d}d} |} || z } d}!t |« D ]_ } | }"|"|"z } | |z } | | z } t | | z |
|
k( « st t | « |"k\ « r n>| |z }t ||z | k( « r n&|dz }| | z } a t || |j fz « |d k7 r |dz }#| | z }$n!d}#|}%| |%k\ r|%|z }%|#dz }#| |%k\ r|%|%z dz
|
}$t |« D ]m } |"}&|"|z }"|"|z } |"| z } t | | z |
|
k7 « rHt t |"« |&k « r1| dz } | |z }t ||z |"k( « sXt | |"k7 « sgd}!|"}& n n t || |j fz « | }'|$| | z dz
|
k r|d k7 r
|
|$|$z }$|#dz }#|$|'z }(||!z }|dk\ r|(dz
|
}(|(|'z }|dk( r|s|(dz
|
}(|d
|
kD r|(dz
|
}(t | |"k7 « r|(dz
|
}(||z
|
})t |)|z |)k7 « r|||z z
|
})|)|&|z |z |z z })|(|'z dz }t |« D ] }*|dk( r|)|)z })|)|z }) t |&||z z « }+|| _ || _ || _ ||« | _ ||« | _
|
|| _ || « | _ || « | _ | | _ |#| _ |'| _ ||&« | _ ||&« | _ |(| _ ||)« | _ ||)« | _ || _ || _ | j | _ | j" | _ | j( | _ | j" | _ || j" « | _ ||+« | _ ||+« | _ dd l },tA |,jC || j « « « | _" | | z | z | z | z }-|-| jD z }. ||.« | _# ||.« | _$ y ) Ni' z'Did not converge after %d tries with %sr é éÿÿÿÿé é zAcould not determine machine tolerance for 'negep', locals() -> %sé
|
é )%Úranger Ú RuntimeErrorÚdtypeÚlocalsÚabsÚibetaÚitÚnegepÚepsnegÚ _str_epsnegÚmachepÚepsÚ_str_epsÚngrdÚiexpÚminexpÚxminÚ _str_xminÚmaxexpÚxmaxÚ _str_xmaxÚirndr ÚepsilonÚtinyÚhugeÚsmallest_normalÚ_str_smallest_normalÚsmallest_subnormalÚ_str_smallest_subnormalÚmathÚintÚlog10Ú precisionÚ
|
resolutionÚ_str_resolution)/r r r r r r Ú max_iterNÚmsgÚoneÚtwoÚzeroÚaÚ_ÚtempÚtemp1ÚbÚitempr$ Úbetar% Úbetahr4 Útempar& ÚbetainÚir' r) r* r, ÚkÚzÚtÚnxresÚyr- ÚmxÚizr/ r. r1 r2 Újr: r< Útenr@ s/ r
|
r zMachAr._do_inity sÝ Ø Ø7ÙmØCiØSyð ÜyÓ!ò 5AØAAØs7DØ1HEÜ53; $Ñ&Õ'Ùð 5ô s a¨¯© ^Ñ 3Ó4Ð 4Ø ÜyÓ!ò 5AØAAØq5DÙ T AXÓ&EÜ5A:Ùð 5ô s a¨¯© ^Ñ 3Ó4Ð 4ØÙ%Ó ð Ø ÜyÓ!ò 5AØaBØDAØs7DØ1HEÜ53; $Ñ&Õ'Ùð 5ô s a¨¯© ^Ñ 3Ó4Ð 4às
|
Ø ÜyÓ!ò 5AØAAØs7DØ1HEÜ53; $Ñ&Õ'Ùð 5ô s a¨¯© ^Ñ 3Ó4Ð 4Ø5yØÜ tax4Ñ Ô ØDØDØu}Ø 19T E\¨TÑ 1Ô2ØDð QØtØ Üu ò AØF
|
Að à ÜyÓ!ò 5AØ7DÜ4#: Ñ%Ô&ÙØDAØAIEàqyÜ"ð $AÜDJÃHñ$Nó Oð Oð 5ô s a¨¯© ^Ñ 3Ó4Ð 4ØØð Ø äyÓ!ò 5AØ7DÜ4#: Ñ%Ô&ÙØDAØaZFð 5ô s a¨¯© ^Ñ 3Ó4Ð 4Øð ØSyØ 19T CZ¨#Ñ -°Ñ 5Ô6ØDð Ø Ø Ø #IØÜyÓ!ò 5AØAØAAØCAØq5DÜ1q5D=Ô!¤S¬¨Q«°1©Ô%5ÙØ6MEÜ54< 1Ñ$Ô%ÙØAAØAAð 5ô s a¨¯© ^Ñ 3Ó4Ð 4Ø B;Øq5DØQBàDØBØr'Ø%ZØaxð r'ð b1Bô yÓ!ò 5AØDØF
|
AØCAØq5DÜAEd?Ô#¬¬C°«F°T©MÔ(:ØEØ v Üut| qÑ(Õ)¬c°$¸!±)nØ EØ DÙáð 5ô s a¨¯© ^Ñ 3Ó4Ð 4Øð Q?u¨{ØbBØ!8DØfØe|Ø 19ØaZFØ VOØ A:aØaZFØ r6ØaZFÜ qAv;ØaZFØV|Ü tczTÑ!Ô "Ø Ñ&DØtd{ TÑ)¨DÑ0Ñ1Ø VOaÑ Üqò #AØzØd{àd{ð #ô ! ¨°© Ñ!4Ó 5Ðà
|
ØØ
|
Ù$ VÓ, Ù'¨Ó/ÔØ Ù! #Ó&Ù$ SÓ) Ø Ø Ø Ù" 4Ó( Ù% dÓ+Ø Ù" 4Ó( Ù% dÓ+Ø à
|
àxx ØII ØII Ø #yyÔ Ù$0°·±Ó$;Ô!Ù"0Ð1CÓ"DÔ Ù'3Ð4FÓ'GÔ$ãÜ djj© ¸¿¹Ó)AÓ BÐ BÓCØCi#o Ñ# cÑ)Ødnn_Ñ-
|
Ù(¨Ó4Ù +¨JÓ 7Õ r c ó$ d}|| j z S )NaC Machine parameters for %(title)s
|
---------------------------------------------------------------------
|
ibeta=%(ibeta)s it=%(it)s iexp=%(iexp)s ngrd=%(ngrd)s irnd=%(irnd)s
|
machep=%(machep)s eps=%(_str_eps)s (beta**machep == epsilon)
|
negep =%(negep)s epsneg=%(_str_epsneg)s (beta**epsneg)
|
minexp=%(minexp)s xmin=%(_str_xmin)s (beta**minexp == tiny)
|
maxexp=%(maxexp)s xmax=%(_str_xmax)s ((1-epsneg)*beta**maxexp == huge)
|
smallest_normal=%(smallest_normal)s smallest_subnormal=%(smallest_subnormal)s
|
---------------------------------------------------------------------
|
)Ú__dict__)r Úfmts r
|
Ú__str__zMachAr.__str__R s ð Uð ð T]]Ñ"Ð"r N) Ú__name__Ú
|
__module__Ú __qualname__Ú__doc__Úfloatr= r r r^ r r r
|
r r s+ ñTðl #(°#Ø %Ù 5Ø5óUò$W8ór #r Ú__main__N) rb Ú__all__Ú _ufunc_configr Ú fromnumericr r r_ Úprintr r r
|
ú<module>ri s> ðñð *å #Ý ÷
|
O#ñ O#ðd
|
zÒÙ &(
Oð r
|
