diff options
Diffstat (limited to 'gsl-1.9/doc/examples/vectorview.c')
-rw-r--r-- | gsl-1.9/doc/examples/vectorview.c | 30 |
1 files changed, 30 insertions, 0 deletions
diff --git a/gsl-1.9/doc/examples/vectorview.c b/gsl-1.9/doc/examples/vectorview.c new file mode 100644 index 0000000..9225c28 --- /dev/null +++ b/gsl-1.9/doc/examples/vectorview.c @@ -0,0 +1,30 @@ +#include <math.h> +#include <stdio.h> +#include <gsl/gsl_matrix.h> +#include <gsl/gsl_blas.h> + +int +main (void) +{ + size_t i,j; + + gsl_matrix *m = gsl_matrix_alloc (10, 10); + + for (i = 0; i < 10; i++) + for (j = 0; j < 10; j++) + gsl_matrix_set (m, i, j, sin (i) + cos (j)); + + for (j = 0; j < 10; j++) + { + gsl_vector_view column = gsl_matrix_column (m, j); + double d; + + d = gsl_blas_dnrm2 (&column.vector); + + printf ("matrix column %d, norm = %g\n", j, d); + } + + gsl_matrix_free (m); + + return 0; +} |