1 files changed, 33 insertions, 0 deletions
diff --git a/gsl-1.9/doc/examples/cblas.c b/gsl-1.9/doc/examples/cblas.c
new file mode 100644
index 0000000..cb6b31a
--- /dev/null
+++ b/gsl-1.9/doc/examples/cblas.c
@@ -0,0 +1,33 @@
+#include <stdio.h>
+#include <gsl/gsl_cblas.h>
+
+int
+main (void)
+{
+  int lda = 3;
+
+  float A[] = { 0.11, 0.12, 0.13,
+                0.21, 0.22, 0.23 };
+
+  int ldb = 2;
+  
+  float B[] = { 1011, 1012,
+                1021, 1022,
+                1031, 1032 };
+
+  int ldc = 2;
+
+  float C[] = { 0.00, 0.00,
+                0.00, 0.00 };
+
+  /* Compute C = A B */
+
+  cblas_sgemm (CblasRowMajor, 
+               CblasNoTrans, CblasNoTrans, 2, 2, 3,
+               1.0, A, lda, B, ldb, 0.0, C, ldc);
+
+  printf ("[ %g, %g\n", C[0], C[1]);
+  printf ("  %g, %g ]\n", C[2], C[3]);
+
+  return 0;  
+}