Loading...

help-gsl@gnu.org

[Prev] Thread [Next]  |  [Prev] Date [Next]

Re: [Help-gsl] GSL Multidimensional Minimizer accepting a matrix instead of vector ? Rhys Ulerich Tue Feb 07 12:00:24 2012

> I am trying to write a minimizer for a function that takes a data
> matrix (N x p) and a weight vector (P-dimensional) and a class table
> vector Y (N -dim). The Multidimensional Minimizer in GSL only allows a
> vector data type for the functions  gsl_multimin_function and
> gsl_multimin_function_fdf. The two-dimensional paraboloid example
> shows how to pass a single vector
> of data (x,y). In my case I have a matrix. Is it possible to do such
> minimization with gsl ?

Can you pack your matrix, weight vector, and vector into a contiguous
region of memory and then use a gsl_vector_view
(http://www.gnu.org/software/gsl/manual/html_node/Vector-views.html)
to view that memory as a single "vector" to accomplish what you want?

- Rhys