LAPACK(3)——线性方程组求解

解线性方程组A‘x=B使用函数LAPACKE_dgesv，函数原型如下，

lapack_int LAPACKE_dgesv( int matrix_order, lapack_int n, lapack_int nrhs,
double * a, lapack_int lda, lapack_int * ipiv,
double * b, lapack_int ldb );
/*

// 参数说明：
matrix_order 矩阵的顺序，两种取值，LAPACK_COL_MAJOR或是LAPACK_ROW_MAJOR。

　　　　　　　　　 　　如果是COL_MAJOR，表示a[0:N-1]是第一列，a[N:2N-1]是第二列，表明矩阵的元素存储顺序。

　　　　　　　　　　　LAPACK的优先选择是COL_MAJOR，当是ROW_MAJOR时，会先进行转置，然后进行计算，然后转置回来；
n 　　　　矩阵的size，系数矩阵应该是方阵，一个5x5的矩阵，则n=5；
nrhs 　　rhs的列数，LAPACK可以同时对一个系数矩阵，多个rhs进行求解，

　　　　　　　　　　　因为对系数矩阵只需要进行一次LU分解，多个rhs一起求解更方便；
a 　　　　系数矩阵；
lda 　　 leading dimension of a，lda>=max(1,n)；
ipiv 　　pivot的行交换记录，第i行被交换到第ipiv[i]行；
b 　　　 rhs，可以是多列。其元素存储顺序是
ldb 　　 leading dimension of b, ldb>=max(1,n)，

　　　　　　　　　　　本来看代码中应该是ldb>=max(1,nrhs)的，ldb=nrhs的时候报错了；
// 返回值：
info = 0, 程序正常运行结束。
*/

测试代码如下，参考 http://blog.sina.com.cn/s/blog_441ac3c9010000v0.html

#include < stdio.h >

// lapacke headers
#include " lapacke.h "
#include "