24小时热门版块排行榜    

查看: 2841  |  回复: 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 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyz1981

至尊木虫 (文坛精英)

医学物理

关于LAPACK

官方网站:
http://netlib.amss.ac.cn/lapack/index.html

用户手册:
http://netlib.amss.ac.cn/lapack/lug/index.html
http://www.netlib.org/lapack/lawn41/index.html


LAPACK (Linear Algebra PACKage) 是Oak Ridge 国家实验室、加州大学Davis 分校和Illinois 大学等联合开发的线性代数函数库,用于在不同高性能计算环境上高效求解数值线性代数问题[7] 。

LAPACK 采用标准Fortran 77 编写。LAPACK 支持实型和复型数据类型,完全支持单精度和双精度计算。LAPACK 可以在向量机, 高性能超标量工作站,和共享存储多处理机上高效运行,也可以在各种类型的单机(PC, 工作站, 大型机)上获得满意的结果。LAPACK 最新版本为lapack-3.0。

  

测试文件:Single Precision Example (sgeev):
program test_sgeev
c Purpose: Compute the eigenvalues and, optionally, the left and/or
c right eigenvectors for an n-by-n real nonsymmetric matrix.
c Usage: call sgeev(jobvl, jobvr, n, a, lda, wr, wi, vl, ldvl, vr, ldvr,
c work, lwork, info)

integer nout, n, lda, lwork, ldvr, ldvl
parameter (nout=6, n=3, lda=n, lwork=5*n, ldvl=n, ldvr=n)
integer i, j, info
real a(lda,n), vr(ldvr,n), vl(ldvl,n)
real work(lwork), wi(n), wr(n)
external sgeev
c
data a / 3, 7, -4, 2, 8, -4, -1, -1, 3 /
c
call sgeev('N', 'V', n, a, lda, wr, wi, vl, ldvl, vr, ldvr,
& work, lwork, info)
c
write (*,*)
print*,' *** the real part of the computed eigenvalues ***'
write(nout, 98) (wr(i), i=1,n)
c
write (*,*)
print*,' *** the imaginary part of the computed eigenvalues ***'
write(nout, 98) (wi(i), i=1,n)
c
write (*,*)
print*,' *** the right eigenvectors ***'
do 10 i = 1,ldvr
write(nout, 99) vr(i,1), vr(i,2), vr(i,3)
10 continue
98 format (1x,f10.4)
99 format (1x,3f10.4)
c
end


***** OUTPUT: *****

*** the real part of the computed eigenvalues ***
11.0000
1.0000
2.0000

*** the imaginary part of the computed eigenvalues ***
0.0000
0.0000
0.0000

*** the right eigenvectors ***
-0.2673 -0.7071 -0.3714
-0.8018 0.7071 0.5571
0.5345 0.0000 0.7428


*********************


DGGEV 广义特征值求解,测试程序
     implicit none
      integer     i,info
      real(8)       ::  aa(2,2),bb(2,2),x(2),vl(2,2),vr(2,2)
      real(8)       ::  WORK(16,2),A1(2), A2(2), B1(2)
      data bb/10., 0., 0., 10./
      data aa/20000, -10000, -10000, 10000/
      call DGGEV('N','V',2,AA,2,BB,2,A1,A2,B1,VL,2,VR,2,WORK,16,info)
      do i=1,2
            x(i)=A1(i)/B1(i)
      end do
      write(*,*)x
      end

解:2628   382


Windows环境编译方法
http://jsy.dyn.dhs.org/maths/for ... dows/#installLAPACK

转自:http://zhousicheng.googlepages.com/a_042
33楼2010-01-29 15:02:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 38 个回答

亲核取代

铁杆木虫 (知名作家)

-无关讨论,删除-

[ Last edited by nono2009 on 2010-1-29 at 10:11 ]
2楼2010-01-28 11:48:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)

感谢nono2009,余泽成领导的支持。呵呵。非常感谢。
3楼2010-01-28 13:41:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mao_yijian

铁杆木虫 (著名写手)

-无关讨论,删除-

[ Last edited by nono2009 on 2010-1-29 at 10:12 ]
坚持“四项基本原则”
4楼2010-01-28 13:46:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见