htswanalysis/MACS/lib/gsl/gsl-1.11/doc/examples/vectorview.c

   1 #include <math.h>
   2 #include <stdio.h>
   3 #include <gsl/gsl_matrix.h>
   4 #include <gsl/gsl_blas.h>
   5
   6 int
   7 main (void)
   8 {
   9   size_t i,j;
  10
  11   gsl_matrix *m = gsl_matrix_alloc (10, 10);
  12
  13   for (i = 0; i < 10; i++)
  14     for (j = 0; j < 10; j++)
  15       gsl_matrix_set (m, i, j, sin (i) + cos (j));
  16
  17   for (j = 0; j < 10; j++)
  18     {
  19       gsl_vector_view column = gsl_matrix_column (m, j);
  20       double d;
  21
  22       d = gsl_blas_dnrm2 (&column.vector);
  23
  24       printf ("matrix column %d, norm = %g\n", j, d);
  25     }
  26
  27   gsl_matrix_free (m);
  28
  29   return 0;
  30 }