Re: [Help-gsl] GSL Matrix from preexisting data

[Joakim Hove]

Mathias Wagner <address@hidden> writes:

> Hi,
>
> I would like to use GSL's LU Decomposition. The data are storred in A Blitz++ 
> (http://www.oonumerics.org/blitz/)  Array. I can get access the data in a 
> one-dimensional C-Array (via the .data() function)
> Is there a way to create a GSL Matrix from a one-dimensional C-Array ?

At least with the help of matrix views it is simple:

/----
| gsl_matrix_view A_view;
| gsl_matrix *A;
| 
| A_view = gsl_matrix_view_array(double * data, size_t N1, size_t N2);
| A = &A_view.matrix;
\----

Observe that the construction above does *not* allocate any new data.


HTH - Joakim


-- 
Joakim Hove
hove AT ift uib no
+47 (55 5)8 27 90
http://www.ift.uib.no/~hove/