I probably have a case of the stupids today...but just in case I don't,
here's a question. There seems to be a whole section of
unreachable code in brent.c (see below)
since the first if-else construct always returns, so the next if-else
construct is never executed. Is this intentional?
Take care,
Andrew
gsl-1.10:brent.c (beginning at line 184):
if (u < z)
{
*x_lower = u;
*f_lower = f_u;
return GSL_SUCCESS;
}
else
{
*x_upper = u;
*f_upper = f_u;
return GSL_SUCCESS;
}
if (f_u <= f_w || w == z)
{
state->v = w;
state->f_v = f_w;
state->w = u;
state->f_w = f_u;
return GSL_SUCCESS;
}
else if (f_u <= f_v || v == z || v == w)
{
state->v = u;
state->f_v = f_u;
return GSL_SUCCESS;
}
-----------
_______________________________________________
Help-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/help-gsl
