really small values

[Dee Asbury]

In multiplying a value of xe^-325 with ye^-4, Python is returning zero. How
do I get it to give me back my tiny value?

Thanks!
Dee

"The most exciting phrase to hear in science, the one that
heralds new discoveries, is not 'Eureka!' but 'That's funny...'"

- Isaac Asimov

[mensanator]

On Jul 17, 4:13?pm, "Dee Asbury" <[e-mail]:

In multiplying a value of xe^-325 with ye^-4, Python is returning zero. How
do I get it to give me back my tiny value?

Use the right tool for the right job.

>>import gmpy
help(gmpy.mpf)

Help on built-in function mpf in module gmpy:

mpf(...)
mpf(n): builds an mpf object with a numeric value n (n may be any
Python number, or an mpz, mpq, or mpf object) and a
default
precision (in bits) depending on the nature of n
mpf(n,bits=0): as above, but with the specified number of bits (0
means to use default precision, as above)
mpf(s,bits=0,base=10): builds an mpf object from a string s made
up of
digits in the given base, possibly with fraction-part
(with
period as a separator) and/or exponent-part (with exponent
marker 'e' for base<=10, else '@'). If base=256, s must be
a gmpy.mpf portable binary representation as built by the
function gmpy.fbinary (and the .binary method of mpf
objects).
The resulting mpf object is built with a default precision
(in
bits) if bits is 0 or absent, else with the specified
number
of bits.

>>a = gmpy.mpf('3e-325')
a

mpf('3.e-325')

>>b = gmpy.mpf('2e-4')
b

mpf('2.e-4')

>>a*b

mpf('6.e-329')

>
Thanks!
Dee

"The most exciting phrase to hear in science, the one that
heralds new discoveries, is not 'Eureka!' but 'That's funny...'"

- Isaac Asimov

[Zentrader]

On Jul 17, 2:13 pm, "Dee Asbury" <[e-mail]:

In multiplying a value of xe^-325 with ye^-4, Python is returning zero. How
do I get it to give me back my tiny value?

Thanks!
Dee

Also, Python's decimal class allows theoretically unlimited
precision. I have extremely limited knowledge here. It gives the
following for 2**-325. I have no idea if the answer is correct.
You'll have to see if gmpy or decimal works better for this. One
piece of advice is to use whichever exclusively. If you use a float
and then covert to either one, the result will be corrupted.
import decimal
decimal.getcontext().prec = 375 ## set precision at 375
print "2**-325 =", decimal.Decimal(str(2**-325))

2**-325 = 1.46302386084E-98

[mensanator]

On Jul 19, 5:11?pm, Zentrader <[e-mail]:

On Jul 17, 2:13 pm, "Dee Asbury" <[e-mail]:

In multiplying a value of xe^-325 with ye^-4, Python is returning zero. How
do I get it to give me back my tiny value?

Thanks!
Dee

Also, Python's decimal class allows theoretically unlimited
precision. I have extremely limited knowledge here. It gives the
following for 2**-325. I have no idea if the answer is correct.

Looks ok, you should get about 0.3 as many decimal digits as bits.

You'll have to see if gmpy or decimal works better for this. One
piece of advice is to use whichever exclusively. If you use a float
and then covert to either one, the result will be corrupted.
import decimal
decimal.getcontext().prec = 375 ## set precision at 375
print "2**-325 =", decimal.Decimal(str(2**-325))

2**-325 = 1.46302386084E-98

Interestingly, using str() causes you to lose precision.

>>a = 2**-325
a

1.463023860841312e-098

>>b = str(a)
b

'1.46302386084e-098'

Setting the precision to 375 didn't help because you
corrupted the value before you converted it to decimal.

Although the 375 will allow you to do this:

>>a = decimal.Decimal('1e-325')
b = decimal.Decimal('1e-4')
a

Decimal("1E-325")

>>b

Decimal("0.0001")

>>a+b

Decimal("0.000100000000000000000000000000000000000 00000000000000000000000000000000000000000000000000 00000000000000000000000000000000000000000000000000 00000000000000000000000000000000000000000000000000 00000000000000000000000000000000000000000000000000 00000000000000000000000000000000000000000000000000 000000000000000000000000000000000001")
You can also do this in gmpy, although you set precision in bits, not
decimal digits.