CLAPACK的安装与简单使用

 

CLAPACK的安装与简单使用

 

烤鱼片(@eii.dlmu)

cleverysm@163.com

 

 

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首先自从其主页上下载CLAPACK包http://www.netlib.org/clapack/clapack.tgz，解压。在其目录下的INSTALL目录中有make.inc的范本文件，比如在linux下就是make.inc.LINUX，将这个文件拷到CLAPACK目录下并改名为make.inc，并修改此文件中的相应参数，如果需要的话。CLAPACK需要F2CLIBS的libI77.a和libF77.a，一个tmglib和blas，这几个库都包含在了CLAPACK的安装包中，但是blas可以选用其他优化过的版本以得到更好的性能。在此我们使用CLAPACK中的blas。

 

首先编译F2CLIBS，用于将fortran转换为c语言，在CLAPACK目录下，make f2clib，在

CLAPACK/F2CLIBS下就会生成libI77.a和libF77.a。

编译tmglib， 在CLAPACK目录下，make tmglib，会生成CLAPACK/tmglib_LINUX.a。

编译blas， 在CLAPACK目录下，make blaslib，会生成CLAPACK/blas_LINUX.a。

最后是编译CLAPACK，make，会生成CLAPACK/lapack_LINUX.a

 

CLAPACK下的clapack.h就是所需要的头文件，除此之外还需要的一个头文件是F2CLIBS/f 2c .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，表示函数正常退出，否则表示出错

 

在此，我用了高教同济大学版线性代数第二版106页，2题(1)的数据。

 

A:

4,2,－1

3,－1,2

11,0,3

 

B:

2

10

8

 

代码：

#include <iostream>

using namespace std;

#include <F2CLIBS/f 2c .h>

//因为程序是C++，而CLAPACK是C语言写的，所以在此处用extern关键字

extern "C"

{

 

    #include <clapack.h>

}

int main(void) {

   

    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[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 EXIT_SUCCESS;

}

 

程序编译链接的时候，lapack_LINUX.a，tmglib_LINUX.a，blas_LINUX.a，libF77.a，libI77.a 这几个库文件都需要链接进去。

 

运算结果：

-14

56

54