li_ma_significance needs float64 for large numbers
Closed this issue · 0 comments
moralejo commented
noff = 1e7
nexcess = 1e4
print(li_ma_significance(np.float64(noff+nexcess), np.float64(noff / 0.2), 0.2))
print(li_ma_significance(np.float32(noff+nexcess), np.float32(noff / 0.2), 0.2))
2.8861902631323018
3.7806044
May be this is unavoidable, I don't know... or perhaps it can be fixed just by using float64 internally in the calculation?