[Help-gsl] Implementing float gsl_multifit_nlin

[Gordan Bobic]

Hi,

I sort of bit this off in the hope that it's not more than I can chew, 
but it looks like the support for the majority of the required 
components for non-linear fitting are only implemented for gsl_vector 
and gsl_matrix types, rather than gsl_vector_float and gsl_matrix_float.

Is this really the case, or am I missing something?
I'm about a thousand lines into the patch, and more and more stuff keeps 
coming up. :-(

Where are:
compute_diag()
compute_gradient_direction()
compute_trial_step()
compute_rptdx()
scaled_enorm()
enorm()

defined?

I don't suppose anyone else would be interested in this, BTW (or if the 
maintainers are reading this, if they'd be interested in a patch)? My 
motivation us purely speed, as discussed earlier re: float vs. double 
performance. If there isn't much interest, I probably won't bother 
finishing this, as it seems to be heading in the epic direction.

If I do continue with it, then there's possibly APIs to discuss, e.g. 
overloading vs. new functions. At the moment I did everything with new 
functions with _float in the name, but maybe it would be better to 
overload the relevant functions instead.

Finally, has anybody thought about wrapping any of GSL in some classes 
for more sensible organization? Far fetched, I know, but it would likely 
make adding handling of another type considerably easier...

Thanks.

Gordan