24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2547  |  回复: 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锂
[考研] 285求调剂 +7 AZMK 2026-04-04 9/450 2026-04-06 00:06 by 永字号
[考研] 找调剂 +10 楚乔乔 2026-04-01 10/500 2026-04-05 22:19 by syh9288
[考研] 工科08专硕机械275求调剂 +3 AaAa7420 2026-04-02 3/150 2026-04-05 13:26 by jp9609
[考研] 323分(计算机视觉和大模型项目)能直接上手 +3 chaoxiicy 2026-04-01 3/150 2026-04-05 00:50 by chongya
[考研] 311分 22408 求调剂 +3 bing_bot 2026-04-03 3/150 2026-04-05 00:43 by chongya
[考研] 材料调剂 +11 一样YWY 2026-04-02 13/650 2026-04-04 23:10 by 无际的草原
[考研] 求生物学学硕调剂——364分 +7 云朵遛弯指南 2026-04-04 7/350 2026-04-04 22:49 by zhyzzh
[考研] 求调剂 +3 ffyyu 2026-04-02 3/150 2026-04-04 19:03 by 蓝云思雨
[考研] 280求调剂 +21 咕噜晓晓 2026-04-02 22/1100 2026-04-04 11:12 by 猪会飞
[考研] 材料专业383求调剂 +8 郭阳阳阳成 2026-04-03 8/400 2026-04-04 10:29 by Rednal.
[考研] 材料调剂 +11 吴棂颖! 2026-04-03 11/550 2026-04-04 09:56 by 小小树2024
[考研] 学硕机械工程303求调剂 +6 无名所以叫吴明 2026-03-30 7/350 2026-04-03 16:48 by asdfzly
[基金申请] esi高被引论文是不是能对中标有所加分和帮助呢 +5 redcom 2026-04-01 6/300 2026-04-03 15:15 by Howard28
[考研] 专硕085601求调剂 +7 suyifei 2026-04-03 8/400 2026-04-03 14:00 by 欣喜777
[考研] 302求调剂 +9 zyx上岸! 2026-04-02 9/450 2026-04-02 23:07 by 马儿快快地跑
[考研] 一志愿北交大材料工程总分358 +8 cs0106 2026-04-01 9/450 2026-04-02 10:36 by 不吃魚的貓
[考研] 322求调剂 +8 三水sss 2026-04-01 8/400 2026-04-01 10:19 by 唐沐儿
[考研] 0710生物学求调剂! +6 叙述文 2026-03-31 6/300 2026-04-01 09:39 by JourneyLucky
[考研] 一志愿中海洋320化学工程与技术学硕求调剂 +8 披星河 2026-03-30 8/400 2026-03-31 08:53 by lbsjt
信息提示
请填处理意见