Problemi con arrotondamenti galleggianti in Python

Question

Sto riscontrando un problema con np.round, np.around dove non sta arrotondando correttamente. Non riesco a comprendere il codice, perché quando lo faccio impostare manualmente il valore (al contrario di utilizzare i miei dati), le opere di ritorno, ma qui è l'output:

In [177]: a 
Out[177]: 0.0099999998 

In [178]: np.round(a,2) 
Out[178]: 0.0099999998 


In [179]: np.round(a,1) 
Out[179]: 0.0 

Che cosa mi manca? Il dtype di a è float32, devo cambiarlo?

Answer 1

Prova a creare np.float32(0.01) e vedrai la tua risposta. Stai ottenendo la precisione che puoi già.

>>> import numpy as np 
>>> x = 0.01 
>>> epsilon = 0.01 - np.float32(0.01) 
>>> for n in np.arange(x - 10*epsilon, x + 10*epsilon, epsilon): 
...  print(repr(np.float32(n))) 
...  
0.0099999979 
0.0099999979 
0.0099999979 
0.0099999988 
0.0099999988 
0.0099999988 
0.0099999988 
0.0099999998 
0.0099999998 
0.0099999998 
0.0099999998 
0.0099999998 
0.010000001 
0.010000001 
0.010000001 
0.010000001 
0.010000002 
0.010000002 
0.010000002 
0.010000002 

Answer 2

Nota sembra che ci sia un problema con la funzione di pitone round e numpy.float64 tipi. Vedere il seguente esempio:

In [88]: round(np.float64(16.259766999999947), 4) 
Out[88]: 16.259799999999998 

L'unico modo ho potuto risolvere questo problema è quello di convertire il numpy.float64 ad un galleggiante prima di utilizzare la funzione di turno come di seguito:

In [89]: round(float(np.float64(16.259766999999947)), 4) 
Out[89]: 16.2598