zMachAr.j s
y1} z Python floating point numberc t d 5 | ||||| ddd dS # 1 swxY w Y dS )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_convint_convfloat_to_float float_to_strtitles r __init__zMachAr.__init__h s H
%
%
% U U MM*h e T T T U U U U U U U U U U U U U U U U U Us 8<<c
d}d} |d }||z } ||z
}
|} t | D ])} | | z } | |z }
|
| z
}t ||z
|
k r n*t || |j fz |}t | D ]/} ||z }| |z }
||
| z
}t |dk r n0t || |j fz |} || }d}|}t | D ].} |dz }||z }||z }
|
|z
}t ||z
|
k r n/t || |j fz || z }|} t | D ])} | | z } | |z }
|
| z
}t ||z
|
k r n*t || |j fz | |z }
d}t |
| z
|
k rd}| |z }||z }
|dk rt |
|z
|
k rd}|dz }||z }|} t | D ]}| |z } | }t | D ]M} || z
}
t |
|z
|
k r nH| |z } |dz
}|dk r t dt z Nt || |j fz | }| }| dz
}|} t | D ])} || z }
t |
|z
|
k r n$| |z } |dz }*t || |j fz | } d} || z }
|dk rt |
|z |z
|
k rd} d}d} |} || z } d}!t | D ]r} | }"|"|"z } | |z } | | z }
t | | z |
k s t t | |"k r nA|
|z }t ||z | k r n$|dz }| | z } st || |j fz |d k r |dz }#| | z }$n"d}#|}%| |%k r|%|z }%|#dz }#| |%k |%|%z dz
}$t | D ]} |"}&|"|z }"|"|z } |"| z }
t | | z |
k rZt t |" |&k r:| dz } |
|z }t ||z |"k rt |
|"k rd}!|"}& n nt || |j fz | }'|$| | z dz
k r|d k r
|$|$z }$|#dz }#|$|'z }(||!z }|dk r|(dz
}(|(|'z }|dk r|s|(dz
}(|d
k r|(dz
}(t | |"k r|(dz
}(||z
})t |)|z |)k 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 |,! || j | _" | | z | z | z | z }-|-| j" z }. ||. | _# ||. | _$ d S ) Ni' z'Did not converge after %d tries with %sr r zAcould not determine machine tolerance for 'negep', locals() -> %s
)%ranger RuntimeErrordtypelocalsabsibetaitnegepepsneg _str_epsnegmachepeps_str_epsngrdiexpminexpxmin _str_xminmaxexpxmax _str_xmaxirndr epsilontinyhugesmallest_normal_str_smallest_normalsmallest_subnormal_str_smallest_subnormalmathintlog10 precision
resolution_str_resolution)/r r r r r r max_iterNmsgonetwozeroa_temptemp1bitempr% betar& betahr5 tempar' betainir( r* r+ r- kztnxresyr. mxizr0 r/ r2 r3 jr; r= tenrA s/ r r zMachAr._do_initz s6 7jmmCiSy
y!! 5 5AAAs7D1HE53;$&''
sa^ 344 4 y!! 5 5AAAq5D HT!V$$E5A:
sa^ 344 4z% y!! 5 5AaBDAs7D1HE53;$&''
sa^ 344 4s
y!! 5 5AAAs7D1HE53;$&''
sa^ 344 45y tAv~ DDu} 199T%Z4 /009D Qt u AF
AA
y!! 5 5A7D48t#$$
DAAIEqyy" $ADJHH$N O O O sa^ 344 4
y!! 5 5A7D48t#$$
DAaZFF sa^ 344 4 Sy 199T#X^t 3449D
#Iy!!
5
5AA!A#AQ3D1Q3$;
3s1vv{#3#3
6ME5:?##
AAAAA sa^ 344 4 B;;q5DQBBDBr''%Zax r'' b1B y!! 5 5ADF
ACAq5DAEd?##
CFFTM(:(:
E v
uTzQ'' C NN E DE sa^ 344 4 Q??u{{bB!8Dfe| 199aZF VO A::a:aZF r66aZF qAv;; aZFV| tCx4 %f$DtDy~d*+ VOa q # #Azzd{d{ !4 5 5
$nV,, '<// !>#&&$ S))
"N4(( % d++ "N4(( % d++
x I I #y $0L$;$;!"0.1C"D"D '3|4F'G'G$ djj )A)A B B BCCCi#o#c)dn_-
(.44 +|J 7 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__S s U T]""r N) __name__
__module__ __qualname____doc__floatr> r r r_ r
r r r r sn T Tl #( % 4 45U U U U$W8 W8 W8r
#
#
#
#
#r __main__N) rc __all__ fromnumericr
_ufunc_configr _utilsr r r` printr
r r rk s * # # # # # #
O# O# O# O# O# O# O# O#d
z E&&((OOOOO r