Exponentiel de grands nombres négatifs

Question

Comment peut-on calculer un nombre tel que np.exp (-28000) sur Python? La réponse est vers 5E-12161. On m'a dit qu'en raison du format à virgule flottante à double précision, je serais seulement capable de calculer un nombre> 1e-2048

Answer 1

Essayez le module décimal.

Decimal(math.exp(1))**-28000 

Answer 2

Essayez mpmath pour l'arithmétique en virgule flottante avec une précision arbitraire

Edit 1:

>>> import mpmath as mp 
>>> import numpy as np 
>>> a = np.matrix((0,0)) 
>>> print(a) 
[0.0 0.0] 
>>> b = mp.matrix(a.tolist()) 
>>> c = b.apply(mp.exp) 
>>> print(c) 
[1.0] 
[1.0] 

Answer 3

Vous pouvez définir une fonction pour calculer la « bigexponent » et l'appliquer à l'ensemble (le long un axe). Mais, notez que le tableau d'entrée doit être de type int.

# solution courtesy: http://stackoverflow.com/a/43084475/2956066 
In [97]: def bigexp(x): 
    ...:  return Decimal(math.exp(1))**x 


In [98]: np.apply_along_axis(bigexp, 1, arr)  

---

efficacité (par ordre décroissant)

# twice faster than applying along axis 0 
In [115]: %timeit np.apply_along_axis(bigexp, 1, a) 
1000 loops, best of 3: 272 µs per loop 

In [116]: %timeit np.apply_along_axis(bigexp, 0, a) 
1000 loops, best of 3: 465 µs per loop