Open pietroleonardi opened 1 year ago
Hmm, this is rather strange.
Can it be that you run out of live points?
According to the debug.log, you have a likelihood plateau, interestingly not at the beginning of the run but after a few iterations.
Maybe you want to check why all the points then return loglike=-5841.19.
That said, this error message is not very helpful, and it should be handled better. But I don't quite understand how this occurs.
Maybe putting a few prints in these lines could help debug:
"/home/pietro/anaconda3/envs/pandora_env/lib/python3.9/site-packages/ultranest/integrator.py", line 369
File "ultranest/mlfriends.pyx", line 855, in ultranest.mlfriends.MLFriends.compute_enlargement
File "ultranest/mlfriends.pyx", line 388, in ultranest.mlfriends.bounding_ellipsoid
File "ultranest/mlfriends.pyx", line 303, in ultranest.mlfriends.make_eigvals_positive
In particular: Is the list of live points empty? What is the highest level loop where this occurs? In integrator.py, perhaps the integration should stop when the live points cannot be replaced.
Description
I am using UltraNest to try to reproduce a injection retrieval example from Pandora, a open-source exomoon transit detection algorithm. However I am trying to fit the synthetic data with the model produced by the code Gefera, a generic code for light curves of two bodies mutually transiting a star.
My scrypt is extremely similar to the one found here, I just changed the likelihood function in order to create the model using gefera and modified the parameters.
After a few seconds of running, this is the error that I get:
What I Did
I start my scrypt by simulating my timing data, with the algorithm Pandora. Then I define the prior_transform function and the log_likelihood function:
These are the sampler and step.sampler:
I am attaching the debug.log