24小时热门版块排行榜    

查看: 1792  |  回复: 5
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

imation

铁杆木虫 (正式写手)

[交流] castep多核运算的问题

大家好,我有一台服务器,cpu为至强clovertown双路四核,每个核心频率是2ghz,内存667,4条,单条1g,操作系统装过2003中文版,sp1和sp2,也装过英文版sp1,2008英文版sp1
我尝试着运行tutorial里面的co分子adsorption在pb表面,一个核心运行很正常,1个小时就有结果,可是8个核心就总是失败,出现的问题也是千奇百怪:
1)进程castep,共八个,可是运行一会儿以后,自动减少到5个,有一次甚至减少到1个,内存无释放,有时任务失败,有时一直运算,但结果遥遥无期
2)cpu利用率在开始计算时为满负荷,经过1到5分钟,降到1~3%,内存没有释放,过10分钟左右,任务失败
3)尝试改变核心数目,发现1到3个都可以正常运行得到结果,超过3个就失败
4)还有n多错误,其中出现最多一个(每次都是在第一个点时出现),代码如下:

*Warning* max. SCF cycles performed but system has not reached the groundstate.
  
Current total energy, E     =  -5951.893478460     eV
Current free energy (E-TS)  =  -5951.990514506     eV
(energies not corrected for finite basis set)

NB est. 0K energy (E-0.5TS)      =  -5951.941996483     eV

  ****************************************************************************
  Warning: electronic minimisation did not converge when finding ground state.
  ****************************************************************************

Writing model to 1.check
Error in geom_get_forces - electronic_minimisation of current_cell failed
Error in geom_get_forces - electronic_minimisation of current_cell failed
Error in geom_get_forces - electronic_minimisation of current_cell failed
Error in geom_get_forces - electronic_minimisation of current_cell failed
[1] MPI Abort by user Aborting program !
[1] Aborting program!
[2] MPI Abort by user Aborting program !
[2] Aborting program!
[0] MPI Abort by user Aborting program !
[0] Aborting program!
forrtl: severe (47): write to READONLY file, unit 60, file D:\PROGRA~1\Accelrys\MATERI~1.1\Gateway\root_default\dsd\jobs\4GLA1\killfile
Image              PC        Routine            Line        Source            
castepexe_mpi.exe  00AB6CC2  Unknown               Unknown  Unknown
castepexe_mpi.exe  00AB3F50  Unknown               Unknown  Unknown
castepexe_mpi.exe  00A42B9E  Unknown               Unknown  Unknown
castepexe_mpi.exe  00A427BB  Unknown               Unknown  Unknown
castepexe_mpi.exe  00A246D1  Unknown               Unknown  Unknown
castepexe_mpi.exe  009F8E7B  Unknown               Unknown  Unknown
castepexe_mpi.exe  009F8EAF  Unknown               Unknown  Unknown
castepexe_mpi.exe  004FB352  Unknown               Unknown  Unknown
castepexe_mpi.exe  004E5281  Unknown               Unknown  Unknown
castepexe_mpi.exe  004DC8D0  Unknown               Unknown  Unknown
castepexe_mpi.exe  00402353  Unknown               Unknown  Unknown
castepexe_mpi.exe  00ABE578  Unknown               Unknown  Unknown
castepexe_mpi.exe  00A904BB  Unknown               Unknown  Unknown
kernel32.dll       7C82F23B  Unknown               Unknown  Unknown

我尝试增加maximum iterations到1000,Max SCF cycles到5000,八核并行运算,算了2天半,正常,我stop运算,在输出中发现还是有些错误:
WARNING - user ionic constraints and symmetry specified
        - symmetry has precedence over constraints
        - may lead to a conflict?

HINT    - if convergence fails try switching symmetry OFF






BFGS: Warning - trial step suggests complex energy landscape
                 in which simple line minimization will fail.
  
               - This is usually an indication that the forces/streses
                 are not accurate enough. Consider increasing the cutoff
                 energy and/or the electronic convergence tolerance.
  
               - Proceeding with a bisection search to find root instead.




Warning: There are no empty bands for at least one kpoint and spin; this may
          slow the convergence and/or lead to an inaccurate groundstate.
          If this warning persists, you should consider increasing nextra_bands
          and/or reducing smearing_width in the param file.
          Recommend using nextra_bands of 14 to 29.  




BFGS: Warning - Repeated consecutive reset of inverse Hessian
BFGS:           without satisfying convergence criteria which
BFGS:           looks like BFGS has run out of search directions.
BFGS: Warning - Lets try allowing some uphill steps and see if
BFGS:           we can get around this barrier.
BFGS: Warning - It is possible that the system may now converge to
BFGS:           a stationary point OTHER than the desired minimum.
BFGS: Hint    - this may be an indication that either:
BFGS:           a) you are using a poor guess at geom_frequency_est
BFGS:           and/or geom_modulus_est, or
BFGS:           b) you are using unrealistic convergence criteria.
BFGS:           Suggest therefore that you consider changing them!


等等,就不一一列出了
我的问题是,为何单核1个小时就能算出来的project,多核要么不能运算,要么运算起来比单核还要慢得多?是不是多核并行处理的时候数据交换有问题啊?

谢谢,不知道描叙清楚了没有
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

imation

铁杆木虫 (正式写手)

呵呵,这儿高手挺多的,而且很热心

谢谢carlon,我对量子力学不熟悉,我自己主要用discover模块,这个castep程序是帮师兄做的,因为是新机器,一直没发现这个问题,所以出现问题以后,让服务器提供商来搞了几次都没成功

以后多多向大家请教
3楼2007-12-07 22:26:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

carlon

金虫 (小有名气)

★ ★
csfn(金币+2,VIP+0):感谢积极的交流 :-) 欢迎常来
多个cpu运行的话,容易出错。只要其中一个进程出错,计算就不会有结果
此时其他进程仍在运行,这可能就是为什么1个cpu1个小时算完,多cpu运行算很久没有结果的原因。
另外,在多个cpu运行过程中,上面还有很多warning信息
1。HINT    - if convergence fails try switching symmetry OFF
这里可能是限制对称性,导致不好收敛,或者给原子设置了constrain,有时候收敛困难,可以在计算的时候去掉symmetry,试试看。
2。BFGS: Warning - trial step suggests complex energy landscape
                 in which simple line minimization will fail.
平面波cutoff设置可能过小。
3。empty band空带数不足,如果体系是金属性的话,不容易收敛。
2楼2007-12-07 21:25:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cometring

木虫 (著名写手)

爱你就等于爱自己

这里确实高手很多
与超版斗,其乐无穷;与斑竹斗,其乐无穷;与虫子斗,其乐无穷。
4楼2007-12-08 14:48:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kuner

铁虫 (正式写手)

我也有这样的问题,单核比多核还快
6楼2008-10-19 00:23:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 YXM1FDCSE8Bn 2026-10-31 6/300 2026-11-02 12:51 by n0P4fnsXKbBK
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 YXM1FDCSE8Bn 2026-10-31 7/350 2026-11-02 12:46 by n0P4fnsXKbBK
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 YXM1FDCSE8Bn 2026-10-31 6/300 2026-11-02 12:44 by n0P4fnsXKbBK
[考博] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 12:25 by n0P4fnsXKbBK
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 YXM1FDCSE8Bn 2026-11-01 6/300 2026-11-02 09:57 by dcc2yOIRGGTT
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-02 09:29 by dcc2yOIRGGTT
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-02 06:07 by xzcv0Tehd8lQ
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-02 06:00 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 xvJS1BYDxvBS 2026-10-31 11/550 2026-11-02 06:00 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-02 05:34 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-02 03:39 by xzcv0Tehd8lQ
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:00 by xzcv0Tehd8lQ
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 02:47 by xzcv0Tehd8lQ
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:27 by xzcv0Tehd8lQ
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 16:52 by BCdvBAEpt3zm
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +9 zuesli 2026-10-27 15/750 2026-11-01 12:59 by Equinoxhua
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 12:27 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
信息提示
请填处理意见