Ñò
[ÐKc @ sb d Z d g Z d d k l Z d d k l Z d e f d „ ƒ YZ e d j o e ƒ GHn d S( sX
Machine arithmetics - determine the parameters of the
floating-point arithmetic system
t MachAriÿÿÿÿ( t any( t seterrc B s; e Z d Z e e e d „ d d „ Z d „ Z d „ Z RS( sò
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, substracted
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] usable floating value).
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
Same as `xmin`.
huge : float
Same as `xmax`.
precision : float
``- int(-log10(eps))``
resolution : float
``- 10**(-precision)``
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 C s d | S( s %24.16e( ( t v( ( s6 /usr/lib64/python2.6/site-packages/numpy/lib/machar.pyt a s s Python floating point numberc C s> t d d ƒ } z | i | | | | | ƒ Wd t | Xd S( s)
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
t undert ignoreN( R t _do_init( t selft
float_convt int_convt float_to_floatt float_to_strt titlet saverrstate( ( s6 /usr/lib64/python2.6/site-packages/numpy/lib/machar.pyt __init___ s
c . C s d } d } | d ƒ } | | } | | }
| } xd t | ƒ D]@ } | | } | | }
|
| } t | | |
j ƒ o Pq? q? Wt | | | i f ‚ | } xf t | ƒ D]B } | | } | | }
| |
| ƒ } t | d j ƒ o Pq¬ q¬ Wt | | | i f ‚ | } | | ƒ } d } | } xn t | ƒ D]J } | d } | | } | | }
|
| } t | | |
j ƒ o Pq3q3Wt | | | i f ‚ | | } | } xd t | ƒ D]@ } | | } | | }
|
| } t | | |
j ƒ o Pq´q´Wt | | | i f ‚ | | }
d } t |
| |
j ƒ o
d } n | | } | | }
| d j o! t |
| |
j ƒ o
d } n | d } | | } | } x t | ƒ D] } | | } q¨W| } x… t | ƒ D]a } | | }
t |
| |
j ƒ o Pn | | } | d } | d j o t d t ƒ ‚ qÏqÏWt | | | i f ‚ | } | } | d } | } xd t | ƒ D]@ } | | }
t |
| |
j ƒ o Pn | | } | d } quWt | | | i f ‚ | } d } | | }
| d j o% t |
| | |
j ƒ o
d } n d } d } | } | | } d }! x½ t | ƒ D]™ } | }" |" |" } | | } | | }
t | | |
j ƒ p t t | ƒ |" j ƒ o Pn |
| } t | | | j ƒ o Pn | d } | | } qFWt | | | i f ‚ | d j o | d }# | | }$ nD d }# | }% x& | |% j o |% | }% |# d }# q-W|% |% d }$ xÒ t | ƒ D]® } |" }& |" | }" |" | } |" | }
t | | |
j ƒ ol t t |" ƒ |& j ƒ oS | d } |
| } t | | |" j ƒ o$ t |
|" j ƒ o d }! |" }& Pq qnPqnWt | | | i f ‚ | }' |$ | | d j o% | d j o |$ |$ }$ |# d }# n |$ |' }( | |! } | d j o |( d }( n |( |' } | d j o | o |( d }( n | d
j o |( d }( n t | |" j ƒ o |( d }( n | | }) t |) | |) j ƒ o | | | }) n |) |& | | | }) |( |' d } x9 t | ƒ D]+ }* | d j o |) |) }) qs|) | }) qsW| | _ | | _ | | _ | | ƒ | _
| | ƒ | _ | | _ | | ƒ | _
| | ƒ | _ | | _ |# | _ |' | _ | |& ƒ | _ | |& ƒ | _ |( | _ | |) ƒ | _ | |) ƒ | _ | | _ | | _ | i
| _ | i | _ | i | _ d d k }+ t |+ i | | i
ƒ ƒ ƒ | _ | | | | | }, |, | i }- | |- ƒ | _ | |- ƒ | _! d S( Ni' s' Did not converge after %d tries with %si i iÿÿÿÿi i sA could not determine machine tolerance for 'negep', locals() -> %si
i (" t xrangeR t RuntimeErrort dtypet ranget localst abst ibetat itt negept epsnegt _str_epsnegt machept epst _str_epst ngrdt iexpt minexpt xmint _str_xmint maxexpt xmaxt _str_xmaxt irndR
t epsilont tinyt huget matht intt log10t precisiont
resolutiont _str_resolution(. R R R
R R R
t max_iterNt msgt onet twot zerot at _t tempt temp1t bt itempR t betaR t betahR& t tempaR t betaint iR R R R t kt zt tt nxrest yR t mxt izR! R R# R$ t jR* t tenR. ( ( s6 /usr/lib64/python2.6/site-packages/numpy/lib/machar.pyR r sr
$
(
0
0
* "
"c C s d | i S( Nsò 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)
---------------------------------------------------------------------
( t __dict__( R ( ( s6 /usr/lib64/python2.6/site-packages/numpy/lib/machar.pyt __str__E s
( t __name__t
__module__t __doc__t floatR+ R R RJ ( ( ( s6 /usr/lib64/python2.6/site-packages/numpy/lib/machar.pyR s O Ót __main__N( RM t __all__t numpy.core.fromnumericR t numpy.core.numericR t objectR RK ( ( ( s6 /usr/lib64/python2.6/site-packages/numpy/lib/machar.pyt s ÿ D