CLAPACK for Windows [转]

下载和教程（光看这两个大概就会安装和是使用了）

　　clapack ： http://icl.cs.utk.edu/lapack-for-windows/clapack/index.html

　　lapack ： http://icl.cs.utk.edu/lapack-for-windows/lapack/#libraries_mingw

　　（把下载的文件放在VS的合适目录下即可）

----

　　

CLAPACK是LAPACK的C语言接口。LAPACK的全称是Linear Algebra PACKage，是非常著名的线性代数库。LAPACK是用Fortran写的，为了方便C/C++程序的使用，就有了LAPACK的C接口库CLAPACK。LAPACK的主页是http://www.netlib.org/lapack/，CLAPACK的主页在http://www.netlib.org/clapack/。

CLAPACK有Linux和Window版。安装CLAPACK-Windows首先从其主页上下载CLAPACK-Windows包，解压，然后用VC（我用的是VS2005)打开clapack.dsw，编绎既可。CLAPACK下的 clapack.h 是所需要的头文件，除此之外还需要的一个头文件是F2CLIBS/f2c.h。

 

现在通过使用CLAPACK中的一个函数sgesv_解线性方程组来学习一下使用的方法。

 

包括此函数在内的所有函数可以在CLAPACK/SRC下找到源代码，并在代码中有函数参数的说明信息。sgesv_的代码文件就是sgesv.c。 

int sgesv_(integer *n, integer *nrhs, real *a, integer *lda, integer *ipiv, real *b, integer *ldb, integer *info)

sgesv_的功能是使用LU分解法解线性方程组AX=B，其中A是一个n*n的方阵。

integer *n, 方程的个数，也既A的行数和列数

integer *nrhs, B的列数

real *a, 存储矩阵A数据的一维数组，在fortran中，数组是列主序存储，在此a中的二维数据也必须是列主序

integer *lda, 等于n

integer *ipiv, 一个输出数据数组，数组大小是n，具体什么功能不太明白，但是似乎不影响最后结果，谁明白请告诉我

real *b,存储矩阵B数据的一维数组，在fortran中，数组是列主序存储，在此b中的二维数据也必须是列主序

integer *ldb, 等于n

integer *info,输出参数，如果返回此参数为0，表示函数正常退出，否则表示出错

 

以下是代码：

// mytest.cpp : Defines the entry point for the console application.
//

#include "stdafx.h"

#include <iostream>

using namespace std;

#include <F2CLIBS/f2c.h>

//因为程序是C++，而CLAPACK是C语言写的，所以在此处用extern关键字
extern "C"
{
#include <clapack.h> 
}


int _tmain(int argc, _TCHAR* argv[])
{
integer M=3 ;
integer N=1;
real a[9]={4,3,11,2,-1,0,-1,2,3};
real b[3]={2,10,8};
integer lda;
integer ldb;
integer INFO;

lda=M;
ldb=M; 
integer *ipiv;
ipiv = (integer *)new integer[M];

sgesv_(&M, &N, a, &lda,ipiv, b, &ldb, &INFO);

if(INFO==0)
{
   for(int i=0; i<M; i++)
   {      
    cout<<b[i]<<endl;
   } 
}
else
{
   cout<<"Failed."<<endl;
}  

return 0;
}

 

编译链接时要将clapack.lib加进链接库里。

 

[转]http://hi.baidu.com/xuyungui/blog/item/c43ef9026e8cff713912bbd3.html