Commit 2b2aee3f authored by Marco Ratto's avatar Marco Ratto
Browse files

fixed bug in sampling in the neighborhood of a critical point

parent 80479f2d
Pipeline #7072 passed with stages
in 136 minutes and 11 seconds
......@@ -203,7 +203,7 @@ if fload==0
for j=1:np
ub=xparam1(j+nshock)*(1+sign(xparam1(j+nshock))*neighborhood_width);
lb=xparam1(j+nshock)*(1-sign(xparam1(j+nshock))*neighborhood_width);
if bounds.ub(j+nshock)>=xparam1(j) && bounds.lb(j)<=xparam1(j+nshock)
if bounds.ub(j+nshock)>=xparam1(j+nshock) && bounds.lb(j+nshock)<=xparam1(j+nshock)
ub=min([bounds.ub(j+nshock) ub]);
lb=max([bounds.lb(j+nshock) lb]);
else
......
Supports Markdown
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