GSL矩阵计算之行列式2

#include <stdio.h>#include <gsl/gsl_cblas.h>#pragma comment(lib, "libgsl_d.lib")#pragma comment(lib, "libgslcblas_d.lib")intmain (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;  }