24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2548  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 3-溴吡啶-4-甲醛合成 10+3 雪糕也能爆炸 2026-03-31 9/450 2026-04-06 00:24 by nBu锂
[考研] 材料专硕(0856) 339分求调剂 +9 哈哈哈鹅哈哈哈 2026-04-05 9/450 2026-04-05 22:24 by dongzh2009
[考研] 315求调剂 +5 &123456789 2026-04-05 5/250 2026-04-05 19:55 by nepu_uu
[考研] 308求调剂 +4 maverick^_^ 2026-04-03 4/200 2026-04-05 19:08 by 蓝云思雨
[考研] 求调剂 +4 wos666 2026-04-03 4/200 2026-04-05 11:48 by arrow8852
[考研] 285求调剂 +11 哦呦呼o 2026-04-04 11/550 2026-04-05 08:15 by 544594351
[考研] 一志愿沪9,求生物学调剂,326分 +6 刘墨墨 2026-04-04 6/300 2026-04-04 19:44 by 唐沐儿
[考研] 278求调剂 +6 Yy7400 2026-04-03 6/300 2026-04-04 09:53 by zhangdingwa
[考研] 282求调剂 +20 ycy1201 2026-04-01 22/1100 2026-04-04 00:42 by userper
[考研] 294求调剂 +6 Grey_Ey 2026-04-03 6/300 2026-04-03 20:46 by 欣喜777
[考研] 335求调剂 +7 沈清璃 2026-04-03 7/350 2026-04-03 18:55 by lijunpoly
[考研] 机械专硕297 +3 Afksy 2026-04-03 3/150 2026-04-03 14:24 by 1753564080
[考研] 一志愿深大085601材料工程专业(专硕)300分可以调剂去哪 +8 10160315 2026-04-02 8/400 2026-04-03 09:36 by hypershenger
[考研] 298求B区调剂 +4 zzz,,r 2026-04-02 5/250 2026-04-02 12:17 by 土木硕士招生
[考研] 085900土木水利336分求调剂 +4 Zhangjiangj 2026-03-31 6/300 2026-04-02 11:40 by 1753564080
[考研] 07生物学求调剂 一志愿同济大学359分 +3 LAMC. 2026-03-30 3/150 2026-04-02 10:26 by 18828373951
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 22:25 by zhouyuwinner
[考研] 292求调剂 +17 木虫er12138 2026-04-01 17/850 2026-04-01 21:37 by 七度不信任
[考研] 0703一志愿南师大334求调剂 +4 seven7yu 2026-03-30 4/200 2026-04-01 16:10 by oooqiao
[考研] 生物考研337分求调剂 +4 cgxin 2026-03-30 6/300 2026-03-31 14:18 by 记事本2026
信息提示
请填处理意见