Commit 4d9dbccf authored by Jakob Stierhof's avatar Jakob Stierhof

Bug in reduce width initialization

Hopefully this works, not sure what I was thinking...
parent 44514cf9
Pipeline #6492 passed with stages
in 42 seconds
......@@ -990,24 +990,22 @@ private define emceeInitUniformPick (init, engine) %{{{
variable file = engine.leader.inFile;
variable par = __parameters(engine.fit.object);
variable numParameter = length(par.value);
variable want_t, max_t,t; % calculate where we want to move a smaller
% box, and how much we can shift it at most
variable minb, maxb;
if (file.mode & EMCEE_FILE_RANGE)
(feed, ) = file.read(engine, engine.totalNumberWalkers);
else
feed = @par;
want_t = par.value-init.width*0.5*(par.max-par.min)+par.min;
max_t = (1-init.width)*(par.max-par.min)+par.min;
t = _max(_min(want_t, max_t), 0.0);
minb = _max(par.min, par.value-init.width*0.5*(par.max-par.min));
maxb = _min(par.max, par.value+init.width*0.5*(par.max-par.min));
(p,c) = __emcee_params(par, feed, engine.totalNumberWalkers);
k = 0;
_for i (0, length(c)-1) {
_for j (0, c[i]-1)
engine.walkers[j+k] = _max(p[i].min, _min(p[i].max, rand_uniform(numParameter)*init.width*(p[i].max-p[i].min)+t[i]));
engine.walkers[j+k] = _max(p[i].min, _min(p[i].max, rand_uniform(numParameter)*(maxb-minb)+minb));
k += c[i];
}
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment