Hi Gideon, > I'm working on a project that involves both numerical integration of a system > of ODEs and the solution of nonlinear systems of equations, and I was > wondering if there was a reason that the interfaces were not better unified. > In particular, the ODE routines take, as arguments, double arrays, while the > nonlinear solvers use the gsl_vector data structure. It's easy enough to map > from one to the other, but I was wondering, why not have more consistency?
>From what I have gathered by lurking on the mailing list and digging through documentation, for logic which could be isolated from the gsl_vector classes, double* arrays were preferred so that folks could copy'n'paste the code without needing all of the GSL. It is likely the case that the nonlinear solvers require enough GSL machinery that they could not simply be ripped out. - Rhys
