vlas-sokolov / pyspecnest

Nested sampling of pyspeckit models with PyMultiNest.
MIT License
2 stars 3 forks source link

Bayesian averaging floating point precision #7

Closed vlas-sokolov closed 6 years ago

vlas-sokolov commented 6 years ago

Extreme log-likelihood values get overflown when we attempt to sum their exponentials. However, using bigfloat package is an external library dependency that we can do without.

First course of action should be trying to calculate the e^ln(Z1) + e^ln(Z2) as a float128 type.