Deal with overflow in exp using numpy
You can use the bigfloat package. It supports arbitrary precision floating point operations. http://packages.python.org/bigfloat/ import bigfloat bigfloat.exp(5000,bigfloat.precision(100)) # -> BigFloat.exact(‘2.9676283840236670689662968052896e+2171′, precision=100) Are you using a function optimization framework? They usually implement value boundaries (using penalty terms). Try that. Are the relevant values really that extreme? In optimization it’s not uncommon to minimize log(f). (approximate log … Read more