[Help-gsl] define various gsl_vector from column of a gsl_matrix

[Abraham Zamudio]

Hi People ,

I try define a set of "m" gsl_vector from columns of a gsl_matrix (size =
nxm)


My code starting with  :

#include <stdio.h>
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_vector.h>

int main ()
{
int i,j;
int n,m;

/*
El tamaño ideal de una matriz de 10000x10000
*/
printf("Ingresa las componentes\n");
printf("n=");
scanf("%d",&n);
printf("m=");
scanf("%d",&m);


//    printf("Creando gsl_matrix");
gsl_matrix * A = gsl_matrix_alloc(n,m);
//    printf("Creada  gsl_matrix");
gsl_vector * v = gsl_vector_alloc (m);

for (i = 0; i < n; i++)
    for (j = 0; j < m; j++)
        gsl_matrix_set (A, i, j, 0.23 + 100*i + j);

FILE * f = fopen("Matrix_A.txt","wb");
gsl_matrix_fwrite (f, A);
fclose (f);

for (i = 0; i < n; i++)
{
    for (j = 0; j < m; j++)
        gsl_vector_set(v,j,gsl_matrix_get(A, i, j));
        printf ("v(%d) = %lf\n",  j, gsl_vector_get (v, j));
}

gsl_matrix_free (A);
gsl_vector_free (v);

return 0;
     }


-- 
Abraham Zamudio Ch.