24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2394  |  回复: 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 的主题更新
信息提示
请填处理意见