In branch ysosilfit I am trying to generate an example that is modeled after pgQuasarSample.py. This example is called ysoSilFit.py
The very short run with the initial guess announced in line 66 does not work properly. In particular, after succesfully calling opt.optimise in line 67, the procedure breaks on lnprob = opt.sampler.lnprobability, as this somehow gets populated entirely with -inf values. This has all sorts of weird consequences downstream of course.
I tried varying nsamples in line 67, nwalkers in line 31 (30 is the minimum), and also the starting positions of the walkers in line 54, by varying the division factor (currently set to 100). None of it gave any different result. It was quite a lot of work to figure out where the problem was in the first place, as the error message was complaining about something else.
In branch ysosilfit I am trying to generate an example that is modeled after pgQuasarSample.py. This example is called ysoSilFit.py The very short run with the initial guess announced in line 66 does not work properly. In particular, after succesfully calling opt.optimise in line 67, the procedure breaks on lnprob = opt.sampler.lnprobability, as this somehow gets populated entirely with -inf values. This has all sorts of weird consequences downstream of course. I tried varying nsamples in line 67, nwalkers in line 31 (30 is the minimum), and also the starting positions of the walkers in line 54, by varying the division factor (currently set to 100). None of it gave any different result. It was quite a lot of work to figure out where the problem was in the first place, as the error message was complaining about something else.