24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

maomao1210

金虫 (正式写手)

程序强帖: 5
应助: 2 (幼儿园)
金币: 1431.3
散金: 242
红花: 16
沙发: 1
帖子: 991
在线: 441.5小时
虫号: 253215
注册: 2006-05-20
性别: MM
专业: 考古理论

[交流] 【讨论】编译器，数学库安装。已有1人参与

经常遇到编译器的安装，以及各种数学库的安装，还有一些常用软件的安装，最基本的能力。呵呵
如果能有个高手把这个总结一下的话，岂不是很好？让大家都来学习一下。
编译器：pgi, ifort,等等
数学库： blas, lapack, atlas, mkl 等等
常用软件就很多了，分子轨道理论方面的，密度泛函方面的，呵呵
希望能有高手给总结一下。让我们都来学习操作操作。
可惜我是一个linux盲，所以只能在这里呼吁了。呵呵。
为了此个项目，我决定拿出我全部的金币做为资助。
pgi 的安装过程，100 金币
ifort 的安装过程，100金币
mpich 的安装过程，100金币
blas+lapack的安装过程，160金币
atlas的安装过程，100金币
mkl的安装过程，100金币
共660金币。
呵呵。总结后，大家都可以来学习一下。
金币就这么多了，所以常用软件的安装就没有金币贡献了，呵呵，等赚了钱再继续来求贤，呵呵。

[ Last edited by maomao1210 on 2010-1-28 at 10:45 ]

回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

1楼 2010-01-28 10:43:18

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zyz1981

至尊木虫 (文坛精英)

医学物理

应助: 1 (幼儿园)
贵宾: 1.976
金币: 16749.1
散金: 305
红花: 50
沙发: 17
帖子: 19496
在线: 761.5小时
虫号: 828345
注册: 2009-08-15
性别: GG
专业: 肿瘤物理治疗

CLAPACK的安装与简单使用

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/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，表示函数正常退出，否则表示出错

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

A:
4,2,－1
3,－1,2
11,0,3

B:
2
10
8

代码：
#include
usingnamespace std;
#include
//因为程序是C++，而CLAPACK是C语言写的，所以在此处用extern关键字
extern"C"
{

#include
}
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    {
         cout<<    }
}
else
{
   cout<<"Failed."< }

return EXIT_SUCCESS;
}

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

运算结果：
-14
56
54

转自：http://blog.csdn.net/daiyuchao/archive/2008/01/04/2026162.aspx