应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » relax计算停住了。大家看看是怎么回事
251/3返回列表123下一页
查看: 1605  |  回复: 24
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

pariah

铜虫 (小有名气)

[求助] relax计算停住了。大家看看是怎么回事

输入文件如下，如果把mp网格取在3 3 3 1 1 1 那么程序能给出计算结果，如果是6 6 6 1 1 1 就会一直停在一个地方。
CODE:
&control
    calculation='relax'
    restart_mode='from_scratch',
    pseudo_dir = '/home/Xiaole/Desktop/Pwscf/espresso-4.3.2/pseudo/',
    outdir='/home/Xiaole/Desktop/Pwscf/espresso-4.3.2/work/work_bamgf4/tmp/',
    prefix='bamgf4'
    tprnfor = .true.
    tstress = .true.
/
&system
    ibrav=  0, celldm(1) =1.889725989, nat=  12, ntyp= 3, ecutwfc =40.0,
    occupations='smearing',smearing='gaussian',degauss=0.001
/
&electrons
    conv_thr    = 1.D-6
    mixing_beta = 0.5D0
/
&IONS
    ion_dynamics='bfgs'
/
ATOMIC_SPECIES
Ba  137.33 Ba.pbe-nsp-van.UPF
Mg   24.32 Mg.pbe-nsp-bpaw.UPF
  F   19.00  F.pbe-n-van.UPF
CELL_PARAMETERS {cubic}
    2.0594950000    7.2313000000    0.0000000000
   -2.0594950000    7.2313000000    0.0000000000
    0.0000000000    0.0000000000    5.8115100000
ATOMIC_POSITIONS {angstrom}
Ba      0.0000000000    -2.1481299780    -2.6947971870
Mg     -2.0594950000    -1.2510149000     0.0000000000
  F     -2.0594950000    -2.7999593600    -1.1448674700
  F     -2.0594950000    -2.3950065600     1.5144795060
  F     -2.0594950000    -0.3948289800    -1.7911073820
  F      0.0000000000    -1.1295290600    -0.0586962510
Ba      0.0000000000     2.1481299780     0.2109578130
Mg      0.0000000000    -5.9802851000    -2.9057550000
  F      0.0000000000    -4.4313406400     1.7608875300
  F      0.0000000000    -4.8362934400    -1.3912754940
  F      0.0000000000    -6.8364710200     1.1146476180
  F      0.0000000000     1.1295290600     2.8470587490
K_POINTS {automatic}
   6 6 6  0 0 0

输出的最后几行如下
        k(   61) = (  -0.1618520  -0.0230480  -0.0286787), wk =   0.0370370
        k(   62) = (  -0.1618520  -0.0230480  -0.0573574), wk =   0.0370370
        k(   63) = (  -0.1213890  -0.0115240  -0.0286787), wk =   0.0370370
        k(   64) = (  -0.1213890  -0.0115240  -0.0573574), wk =   0.0370370
        k(   65) = (   0.0000000   0.0460959  -0.0286787), wk =   0.0185185
        k(   66) = (   0.0000000   0.0460959  -0.0573574), wk =   0.0185185
        k(   67) = (  -0.2023150  -0.0115240  -0.0286787), wk =   0.0370370
        k(   68) = (  -0.2023150  -0.0115240  -0.0573574), wk =   0.0370370

     Dense  grid:    39989 G-vectors     FFT dimensions: (  60,  60,  45)
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-04-06 09:40:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanyanli0475

铜虫 (小有名气)

【答案】应助回帖

★
感谢参与，应助指数 +1
fzx2008: 金币+1, 谢谢回帖交流！ 2012-04-06 10:03:15
这是属于正常的，6 6 6 大一些，而且还是12个原子，计算会慢很多的。等待吧！如果报错，再想办法解决。
赞一下 回复此楼
什么时候能算完！！！！
2楼2012-04-06 09:59:29
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pariah

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by hanyanli0475 at 2012-04-06 09:59:29:
这是属于正常的，6 6 6 大一些，而且还是12个原子，计算会慢很多的。等待吧！如果报错，再想办法解决。

可是他一步自洽也没开始，这是在干什么呢。等了两三天还是停留在这里。
赞一下 回复此楼
3楼2012-04-06 10:56:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanyanli0475

铜虫 (小有名气)
引用回帖:
3楼: Originally posted by pariah at 2012-04-06 10:56:01:
可是他一步自洽也没开始，这是在干什么呢。等了两三天还是停留在这里。

你用几个核跑的？
赞一下 回复此楼
什么时候能算完！！！！
4楼2012-04-06 12:32:25
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pariah

铜虫 (小有名气)
引用回帖:
4楼: Originally posted by hanyanli0475 at 2012-04-06 12:32:25:
你用几个核跑的？

先用了6个独立核心。关键我觉得它并没有开始scf，而是卡住了。所以感觉等待也不会有什么结果。
赞一下 回复此楼
5楼2012-04-06 13:37:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
LZ用什么编译的PW？
赞一下 回复此楼
思想重于技巧，内涵重于表象
6楼2012-04-06 15:02:02
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pariah

铜虫 (小有名气)
引用回帖:
6楼: Originally posted by souledge at 2012-04-06 15:02:02:
LZ用什么编译的PW？

gfortran+mpich2 并行运行少k点的时候没有事情啊
赞一下 回复此楼
7楼2012-04-06 16:41:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
fzx2008: 金币+2, 专家考核, 谢谢指教 2012-04-06 16:49:14
引用回帖:
7楼: Originally posted by pariah at 2012-04-06 16:41:51:
gfortran+mpich2 并行运行少k点的时候没有事情啊

最新的gfortran似乎对QE 4.3.2有些问题～并且gfortran的速度也比较堪忧～如果是intel的平台，试试ifort吧～昨天测试了ifort和PGI在i5 750的效率，输入文件完全相同的前提下，计算结果也完全一致，但是使用ifort可以在1分钟内跑10个左右SCF迭代，而PGI编译的只能跑8个～
使用的破解版PGI 9，试用版的最新PGI 2012并行编译后的PW有问题，一运行就直接段错误……
赞一下 回复此楼
思想重于技巧，内涵重于表象
8楼2012-04-06 16:46:22
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pariah

铜虫 (小有名气)
引用回帖:
8楼: Originally posted by souledge at 2012-04-06 16:46:22:
最新的gfortran似乎对QE 4.3.2有些问题～并且gfortran的速度也比较堪忧～如果是intel的平台，试试ifort吧～昨天测试了ifort和PGI在i5 750的效率，输入文件完全相同的前提下，计算结果也完全一致，但是使用ifort可 ...

我记得当时用ifort编译过，但是make的时候好像mkl的库之类的出问题了。就用了gortran。。。请问如果我用ifort的话 是不是mpich2在安装的时候也要用ifort来编译呢？
赞一下 回复此楼
9楼2012-04-06 16:49:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
liliangfang: 金币+2, 专家考核, 谢谢指教 2012-04-06 18:07:12
pariah: 金币+5 2012-04-07 09:06:39
引用回帖:
9楼: Originally posted by pariah at 2012-04-06 16:49:13:
我记得当时用ifort编译过，但是make的时候好像mkl的库之类的出问题了。就用了gortran。。。请问如果我用ifort的话 是不是mpich2在安装的时候也要用ifort来编译呢？

可以在编译时使用内部库：
--with-internal-blas
--with-internal-lapack
PWscf自带了这些线性数学库，并且这两个库的质量还很不错。经常用这两个库编译静态的VASP，测试与使用动态MKL连接的VASP计算速度没有什么差别～
赞一下 回复此楼
思想重于技巧，内涵重于表象
10楼2012-04-06 16:54:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 pariah 的主题更新
251/3返回列表123下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭