[Help-gsl] fitting

[Petr Ent]

hi everyone,
is it possible to do following thing with non-linear least squares fitting...?

       double x_init[3] = { 1.0, 0.0, 0.0 };
       gsl_vector_view x = gsl_vector_view_array (x_init, p);
       T = gsl_multifit_fdfsolver_lmsder;
       s = gsl_multifit_fdfsolver_alloc (T, n, p);

for ( int i = 0; i < 100 000; i++ )
       
       gsl_multifit_fdfsolver_set (s, &f, &x.vector);

       do
         {
           iter++;
           status = gsl_multifit_fdfsolver_iterate (s);
     
           if (status)
             break;
     
           status = gsl_multifit_test_delta (s->dx, s->x,1e-4, 1e-4);
         }
       while (status == GSL_CONTINUE && iter < 500);
end of for cycle

or another way .... is it possible to allocate solver, create vector from 
x_init, set this to solver and let it do its job in loop (with same x_init 
vector)? does it somehow change the original x_init vector?

thank you in advance
best regards
Petr Ent