| 查看: 3496 | 回复: 37 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
maomao1210金虫 (正式写手)
|
[交流]
【讨论】编译器，数学库安装。 已有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 ] |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有3人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
申请科研项目纯靠实力能成功吗？
已经有7人回复
-
极限拉扯，如何抉择
已经有14人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
zyz1981
至尊木虫 (文坛精英)
医学物理
- 应助: 1 (幼儿园)
- 贵宾: 1.976
- 金币: 16789.9
- 散金: 305
- 红花: 50
- 沙发: 17
- 帖子: 19504
- 在线: 761.6小时
- 虫号: 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 |
31楼2010-01-29 14:55:24
亲核取代
铁杆木虫 (知名作家)
- 应助: 13 (小学生)
- 贵宾: 0.761
- 金币: 13303.5
- 散金: 3309
- 红花: 17
- 沙发: 130
- 帖子: 9420
- 在线: 2971.6小时
- 虫号: 462242
- 注册: 2007-11-19
- 性别: GG
- 专业: 催化化学
2楼2010-01-28 11:48:04
maomao1210
金虫 (正式写手)
- 程序强帖: 5
- 应助: 2 (幼儿园)
- 金币: 1432.8
- 散金: 242
- 红花: 16
- 沙发: 1
- 帖子: 991
- 在线: 442.1小时
- 虫号: 253215
- 注册: 2006-05-20
- 性别: MM
- 专业: 考古理论
3楼2010-01-28 13:41:45
mao_yijian
铁杆木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 6747.2
- 红花: 1
- 帖子: 2683
- 在线: 132.3小时
- 虫号: 810764
- 注册: 2009-07-17
- 性别: GG
- 专业: 检验医学其他科学问题
4楼2010-01-28 13:46:58
10