24小时热门版块排行榜    

查看: 3963  |  回复: 6
【奖励】 本帖被评价1次,作者zhangguangping增加金币 0.8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

[资源] 使用gfortan,openmpi以及BLAS, BLACS, Lapack, Scalapack编译SIESTA 已有1人参与

1.编译openmpi (http://www.open-mpi.org/software/ompi/v1.6/)
(1)
./configure --prefix=/home/x_guazh/software/openmpi-1.6.5 CC=gcc CXX=g++ F77=gfortran FC=gfortran
(3)
make all install
(4)
打开 ~/.bashrc 添加环境变量:
export PATH=/home/x_guazh/software/openmpi-1.6.5/bin:$PATH
export LD_LIBRARY_PATH= /home/x_guazh/software/openmpi-1.6.5/lib: $LD_LIBRARY_PATH
(5)
source ~/.bashrc

2.编译Lapack (http://www.netlib.org/lapack)
(1)
cp make.inc.example  make.inc
然后编辑Makefile,找到lib选项,编辑为:
lib: lapacklib blaslib lapackelib tmglib variants
这样就把所有的东西都编译出来,默认是不编译blas的
(2)
make
这样在目录下就会出来:
liblapack.a, blaslib.a, lapackelib.a, tmglib.a
3.编译BLACS
(1)
下载MPIBLACS: http://www.netlib.org/blacs/mpiblacs.tgz
下载补丁:http://www.netlib.org/blacs/mpiblacs-patch03.tgz
(2)
tar -zxvf mpiblacs.tgz
tar -zxvf mpiblacs-patch03.tgz
(3)
cd BLACS
cp BMAKES/Bmake.MPI-LINUX  ./Bmake.inc
(4)
vi Bmake.inc (参见:http://www.open-mpi.org/faq/?category=mpi-apps#blacs)

# Section 1:
BTOPdir = /home/x_guazh/software/BLACS
COMMLIB = MPI
# Ensure to use MPI for the communication layer
   COMMLIB = MPI
# The MPIINCdir macro is used to link in mpif.h and
# must contain the location of Open MPI's mpif.h.  
# The MPILIBdir and MPILIB macros are irrelevant
# and should be left empty.
MPIdir = /home/x_guazh/software/openmpi-1.6.5-gnu
MPILIBdir =
MPIINCdir = $(MPIdir)/include
MPILIB =
# Section 2:
SYSINC =
INTFACE = -Df77IsF2C
SENDIS =
BUFF =
TRANSCOMM = -DUseMpi2
WHATMPI =
SYSERRORS =

# Section 3:
# You may need to specify the full path to
# mpif77 / mpicc if they aren't already in
# your path.
   F77            = mpif77
   F77LOADFLAGS   =
   CC             = mpicc
   CCLOADFLAGS    =
(4) make mpi
编译完成后,在LIB下面出现:
blacsCinit_MPI-LINUX-0.a
blacsF77init_MPI-LINUX-0.a
blacs_MPI-LINUX-0.a

4.编译scalapack
参见:http://www.open-mpi.org/faq/?category=mpi-apps#scalapack
(1)
下载:http://www.netlib.org/scalapack/scalapack-2.0.2.tgz
(2)
tar zxvf scalapack-2.0.2.tgz
(3)
cd scalapack-2.0.2
(4)
cp SLmake.inc.example SLmake.inc
(5)
vi SLmake.inc
2.0.2版本的文件设置比较简单

BLASLIB       = /home/x_guazh/software/lapack-3.5.0/librefblas.a
LAPACKLIB     = /home/x_guazh/software/lapack-3.5.0/liblapack.a

5.编译siesta
把上面得到的.a文件全部拷到一个文件中,方便使用,比如~/software/lib
librefblas.a
liblapack.a
blacsCinit_MPI-LINUX-0.a
blacs_MPI-LINUX-0.a  
blacsF77init_MPI-LINUX-0.a
libscalapack.a

part of arch.make:

FFLAGS= -ftree-vectorize -O3
FFLAGS_DEBUG= -g

PATH_LIB=/home/x_guazh/software/lib

BLAS_LIBS=$(PATH_LIB)/librefblas.a
LAPACK_LIBS=$(PATH_LIB)/liblapack.a
BLACS_LIBS=$(PATH_LIB)/blacsCinit_MPI-LINUX-0.a $(PATH_LIB)/blacs_MPI-LINUX-0.a $(PATH_LIB)/blacsF77init_MPI-LINUX-0.a
SCALAPACK_LIBS=$(PATH_LIB)/libscalapack.a

#COMP_LIBS=dc_lapack.a # it it needed only if there is lack of LAPACK/BLAS
MPI_INTERFACE= libmpi_f90.a
MPI_INCLUDE=.
FPPFLAGS_MPI=-DMPI

[ Last edited by zhangguangping on 2014-1-17 at 11:43 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wiseguy2008

新虫 (初入文坛)


如果机器是64位,那中间编译过程有哪些需要修改
7楼2014-02-25 20:12:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

乐石||

木虫 (正式写手)



小木虫: 金币+0.5, 给个红包,谢谢回帖
顶一个!这大半夜的
2楼2014-01-16 08:03:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiongxiong5712

新虫 (小有名气)


先收藏了,guangping兄好勤奋啊,佩服佩服。另外借宝地想请教guangping兄和各位大侠们,用siesta做结构优化的时候有没有用过MD.RelaxCellOnly T这个设置先只优化晶胞参数?我现在在优化一个六方晶系的晶体,其他人一般的做法是:先设置一系列晶格常数,计算总能量,找到总能量最低的晶格常数,然后固定晶胞,只优化原子坐标。我看见siesta手册上有MD.RelaxCellOnly 这个参数,可以只优化晶胞参数,所以就直接拿这个参数来优化晶格常数和夹角,不知道这样做可不可以?
下面是晶胞和坐标一起优化的设置:
使用gfortan,openmpi以及BLAS, BLACS, Lapack, Scalapack编译SIESTA
结果是:
outcell: Cell vector modules (Ang)   :    7.942596    7.942622    5.742192
outcell: Cell angles (23,13,12) (deg):     90.0157     90.0139    120.0227

下面是只优化晶胞参数的设置:
使用gfortan,openmpi以及BLAS, BLACS, Lapack, Scalapack编译SIESTA-1
结果是:
outcell: Cell vector modules (Ang)   :    7.931178    7.931178    5.750268
outcell: Cell angles (23,13,12) (deg):     90.0000     90.0000    120.0000
请问一下这两种方法哪种是正确的,哪种不正确呢?还请赐教,多谢了。
3楼2014-01-16 10:05:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
3楼: Originally posted by xiongxiong5712 at 2014-01-16 03:05:45
先收藏了,guangping兄好勤奋啊,佩服佩服。另外借宝地想请教guangping兄和各位大侠们,用siesta做结构优化的时候有没有用过MD.RelaxCellOnly T这个设置先只优化晶胞参数?我现在在优化一个六方晶系的晶体,其他人一 ...

没有研究过这个选项。等待高人指点!
4楼2014-01-17 18:45:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 11:51 by xiangy672
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见