int
gsl_multilarge_nlinear_init (const gsl_vector * x,
gsl_multilarge_nlinear_fdf * fdf,
gsl_multilarge_nlinear_workspace * w)
{
return gsl_multilarge_nlinear_winit(x, NULL, fdf, w);
}
static void
test_fdf(const gsl_multilarge_nlinear_type * T,
const gsl_multilarge_nlinear_parameters * params,
const double xtol, const double gtol, const double ftol,
const double epsrel, const double x0_scale,
test_fdf_problem *problem,
const double *wts)
{
gsl_multilarge_nlinear_fdf *fdf = problem->fdf;
const size_t n = fdf->n;
const size_t p = fdf->p;
const size_t max_iter = 2500;
gsl_vector *x0 = gsl_vector_alloc(p);
gsl_vector_view x0v = gsl_vector_view_array(problem->x0, p);
gsl_multilarge_nlinear_workspace *w =
gsl_multilarge_nlinear_alloc (T, params, n, p);
const char *pname = problem->name;
char buf[2048];
char sname[2048];
int status, info;
sprintf(buf, "%s/%s/solver=%s/scale=%s%s%s",
gsl_multilarge_nlinear_name(w),
params->trs->name,
params->solver->name,
params->scale->name,
problem->fdf->df ? "" : "/fdjac",
problem->fdf->fvv ? "" : "/fdfvv");
strcpy(sname, buf);
/* scale starting point x0 */
gsl_vector_memcpy(x0, &x0v.vector);
test_scale_x0(x0, x0_scale);
if (wts)
{
gsl_vector_const_view wv = gsl_vector_const_view_array(wts, n);
gsl_multilarge_nlinear_winit(x0, &wv.vector, fdf, w);
}
else
gsl_multilarge_nlinear_init(x0, fdf, w);
status = gsl_multilarge_nlinear_driver(max_iter, xtol, gtol, ftol,
NULL, NULL, &info, w);
gsl_test(status, "%s/%s did not converge, status=%s",
sname, pname, gsl_strerror(status));
/* check solution */
test_fdf_checksol(sname, pname, epsrel, w, problem);
if (wts == NULL)
{
/* test again with weighting matrix W = I */
gsl_vector *wv = gsl_vector_alloc(n);
sprintf(sname, "%s/weighted", buf);
gsl_vector_memcpy(x0, &x0v.vector);
test_scale_x0(x0, x0_scale);
gsl_vector_set_all(wv, 1.0);
gsl_multilarge_nlinear_winit(x0, wv, fdf, w);
status = gsl_multilarge_nlinear_driver(max_iter, xtol, gtol, ftol,
NULL, NULL, &info, w);
gsl_test(status, "%s/%s did not converge, status=%s",
sname, pname, gsl_strerror(status));
test_fdf_checksol(sname, pname, epsrel, w, problem);
gsl_vector_free(wv);
}
gsl_multilarge_nlinear_free(w);
gsl_vector_free(x0);
}
