24小时热门版块排行榜    

查看: 2520  |  回复: 6

lei0819

金虫 (正式写手)

xcpu


[求助] Abinit 并行运算和串行运算的差别

大家好,我是刚深入学习使用Abinit,我在自己的电脑上面串行编译了Abinit的程序,之后又在学校集群上面装了同一个版本并行编译后的Abinit。
我在自己电脑上面试验计算ZrO2例子的时候,发现能正常输出结果和分析结果,用了大概1个多小时。
但是同一个输入文件我上传到集群上,用12个核,并行了将近一天,还是没有结束,感觉里面应该有问题,不得已杀死任务,请问这是怎么回事?

并行提交任务指令:
bsub -q normal -n 12 -o %J.log -e %J.err -i zroo.files mpijob abinit

下面是计算的输入文件:

# ZrO2-cubic
#
# Computation of the band structure.
# First, a SCF density computation, then a non-SCF band structure calculation.

ndtset 2

#Dataset 1 : usual self-consistent calculation

kptopt1 1          # Option for the automatic generation of k points,
                   # taking into account the symmetry
nshiftk1 4
shiftk1  0.5 0.5 0.5  # These shifts will be the same for all grids
         0.5 0.0 0.0
         0.0 0.5 0.0
         0.0 0.0 0.5
ngkpt1  8 8 8  
prtden1  1         # Print the density, for use by dataset 2
toldfe1  1.0d-20    # This value is way too large for most realistic studies of materials
iscf1 5

#Dataset 2 : the band structure
iscf2    -2
getden2  -1
kptopt2  -5
nband2   30
ndivk2   10 12 18 8 8       # 10, 12 and 18,8,8 divisions of the 5 segments, delimited
                            # by 6 points.
tolwfr2 1.0d-20
kptbounds2  0.5  0.25  0.75 # W point
            0.5  0.0  0.0 # L point
            0.0  0.0  0.0 # Gamma point
            0.5  0.0  0.5 # X point
            0.5  0.25 0.75 # W point
            0.0  0.0  0.0 # Gamma point in another cell.
# tolwfr2  1.0d-12
enunit2  0             # Will output the eigenenergies in Har
prteig2  1             #


#Definition of the unit cell
acell 3*9.48196         # This is equivalent to   9.48196 9.48196 9.48196
rprim  0.0  0.5  0.5   # FCC primitive vectors (to be scaled by acell)
       0.5  0.0  0.5   
       0.5  0.5  0.0

#Definition of the atom types
ntypat 2          # There is two types of atom
znucl 8 40        # The keyword "znucl" refers to the atomic number of the
                  # possible type(s) of atom. The pseudopotential(s)
                  # mentioned in the "files" file must correspond
                  # to the type(s) of atom. Here, the only type is Silicon.
                        

#Definition of the atoms
natom 3           # There are 3 atoms
typat 1 1 2       #
xred              # This keyword indicate that the location of the atoms
                  # will follow, one triplet of number for each atom
   0.0  0.0  0.0  # Triplet giving the REDUCED coordinate of atom 1.
   1/4  1/4  1/4  # Triplet giving the REDUCED coordinate of atom 2.
        3/4        3/4      3/4      #
#Definition of the planewave basis set
ecut 40.0         # Maximal kinetic energy cut-off, in Hartree

#Definition of the SCF procedure
nstep 60          # Maximal number of SCF cycles
diemac 3.0       # Although this is not mandatory, it is worth to
                  # precondition the SCF cycle. The model dielectric
                  # function used as the standard preconditioner
                  # is described in the "dielng" input variable section.
                  # Here, we follow the prescription for bulk silicon.
# add to conserve old < 6.7.2 behavior for calculating forces at each SCF step
optforces 1
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

你跨节点计算了没有?
嘟啦啦嘟啦啦嘟
2楼2012-03-14 16:54:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lei0819

金虫 (正式写手)

xcpu



我自己摸索的,都不知道怎么设置跨节点的?
向您请教该怎么设置?
3楼2012-03-14 16:57:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mailoliver

金虫 (小有名气)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
fzx2008: 金币+1, 谢谢指教 2012-03-30 20:35:40
lei0819: 金币+2, ★★★很有帮助, 谢谢提醒,这个参量可能设置有问题 2012-03-31 08:59:24
toldfe1  1.0d-20   
有必要吗?这会搞死自己!
一般1.0e-6 eV/atom就足够了。
4楼2012-03-30 19:57:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
4楼: Originally posted by mailoliver at 2012-03-30 19:57:51:
toldfe1  1.0d-20   
有必要吗?这会搞死自己!
一般1.0e-6 eV/atom就足够了。

5楼2012-03-30 20:36:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
WDD880227: 金币+1, 感谢交流提示 2012-03-31 08:50:34
lei0819: 金币+3, ★★★★★最佳答案, 谢谢您的回答,并行的另外参数设置能具体说说吗?我是生手,自己摸索太费时间了。麻烦您了。 2012-03-31 09:02:35
Abinit的收敛精度和通常理解的优点不一样。1.0-20是正常值。
此外,并行肯定是要快一些的。但是要注意的是,你这个并行只是对k点并行,且最好是cpu数目能被k点数目整除。
真正的并行比这复杂,还需要另外的参数。Abinit算起来真的很慢。
6楼2012-03-31 07:00:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)


fzx2008: 金币+1, 专家考核, 谢谢指教 2012-03-31 12:40:03
ABINIT并行是非常麻烦的,而且需要机器的支持,并不是一定每个机器都可以。
你看看para_kgp 这个参数的设置,就可以发现了,有许多相应的配置,看了说明就明白了
7楼2012-03-31 12:35:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lei0819 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 308求调剂 +3 是Lupa啊 2026-03-16 3/150 2026-03-16 10:07 by 求调剂zz
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +6 Ncdx123456 2026-03-13 7/350 2026-03-16 07:44 by L135790
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 梁成伟老师课题组欢迎你的加入 +6 一鸭鸭哟 2026-03-14 7/350 2026-03-15 22:12 by Winj1e
[考研] 材料工程专硕274一志愿211求调剂 +5 薛云鹏 2026-03-15 5/250 2026-03-15 20:38 by Logic2024
[考研] 0703化学调剂 ,六级已过,有科研经历 +4 曦熙兮 2026-03-15 4/200 2026-03-15 18:01 by JourneyLucky
[考研] 一志愿浙江大学0856材料与化工求调剂 +4 yansheng@211 2026-03-09 5/250 2026-03-14 02:10 by JourneyLucky
[考研] 求调剂 +6 yfihxh 2026-03-09 6/300 2026-03-14 01:18 by JourneyLucky
[考研] 考研材料与化工,求调剂 +8 戏精丹丹丹 2026-03-09 8/400 2026-03-14 01:14 by JourneyLucky
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 一志愿华中师范071000,325求调剂 +5 RuitingC 2026-03-12 5/250 2026-03-13 10:43 by hyswxzs
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
[考研] 一志愿:武汉理工,材料工程,英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
请填处理意见