summaryrefslogtreecommitdiff
path: root/gsl-1.9/doc/examples/vectorview.c
diff options
context:
space:
mode:
Diffstat (limited to 'gsl-1.9/doc/examples/vectorview.c')
-rw-r--r--gsl-1.9/doc/examples/vectorview.c30
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;
+}