24小时热门版块排行榜    

查看: 802  |  回复: 0

cavediger

金虫 (正式写手)

[交流] 关于vasp的bug

以下内容转自:完了,不小心把窗口关了

最近在vasp测试和实施过程中,发现了vasp的一些小bug,总结如下:

1. vasp.4.6.28使用fftw3时,会报段错误。这是因为传递给ftw3函数的一个整形参数长度不够造成的。
涉及到两个源文件:fftw3d.F(串行版),fftmpiw.F(MPI并行版),修正如下:

lib@node203:~/source/vasp> diff vasp.4.6.28/vasp.4.6/fftw3d.F vasp.4.6.28/vasp.4.6/fftw3d.F.orig
64,65c64
<        integer grid(3), isign
<        integer(8) plan
---
>        integer grid(3), plan, isign
91,92c90
<        integer grid(3), isign
<        integer(8) plan
---
>        integer grid(3), isign, plan
120,121c118
<        integer grid(3)
<        integer(8) plan
---
>        integer grid(3), plan


lib@node203:~/source/vasp> diff vasp.4.6.28/vasp.4.6/fftmpiw.F vasp.4.6.28/vasp.4.6/fftmpiw.F.orig
56c56
<       INTEGER(8) :: planx, plany, planz
---
>       INTEGER :: planx, plany, planz


2. vasp.5.2.2,如果用的是PGI/gfortran编译,在计算有Spin-Orbital Coupling (SOC)算例时,会报如下错误:
PGFIO-F-235/formatted write/unit=8/edit descriptor does not match item type.
File name = OUTCAR    formatted, sequential access   record = 1493
In source file relativistic.f, at line number 491
这是由于输出结果时,格式参数与实际数据类型不符合,Intel编译器可能容错比较好,自动做了类型转换,不会报错

修正如下:
lib@node203:~/source/vasp> diff vasp.5.2.2/vasp.5.2/relativistic.F vasp.5.2.2/vasp.5.2/relativistic.F.orig
568c568
<                    -SIN(BETA),0.0,COS(BETA), &
---
>                    -SIN(BETA),0,COS(BETA), &
570c570
<                    -SIN(ALPHA),COS(ALPHA),0.0, &
---
>                    -SIN(ALPHA),COS(ALPHA),0, &

修正后的源代码:
vasp.4.6.28:
vasp.5.2.2:  没权限,拿不到这个!
回复此楼
我们都老得太快,却聪明得太迟。当我一路跌跌撞撞登上山巅,蓦然发现,本来就有一条笔直的路,直达山顶....
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 cavediger 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 情人节自我反思:在爱情中有过遗憾吗? +5 瞬息宇宙 2026-02-15 6/300 2026-02-18 12:51 by 月下雪林
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 11:09 by lqtl9djx19
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:54 by lqtl9djx19
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:39 by lqtl9djx19
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:53 by lqtl9djx19
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:38 by lqtl9djx19
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 07:55 by lotyj5cz79
[基金申请] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:40 by lotyj5cz79
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:38 by lotyj5cz79
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:23 by lotyj5cz79
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:08 by lotyj5cz79
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 3/150 2026-02-18 06:53 by lotyj5cz79
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 00:40 by tk2gfblvuz
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 00:23 by tk2gfblvuz
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-17 23:40 by tk2gfblvuz
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +3 successhe 2026-02-16 4/200 2026-02-17 20:56 by successhe
[基金申请] 今年春晚有几个节目很不错,点赞! +5 瞬息宇宙 2026-02-16 6/300 2026-02-17 12:49 by jymy19840415
[微米和纳米] 球磨粉体时遇到了大的问题,请指教! 10+3 6sbiam 2026-02-12 15/750 2026-02-16 15:03 by tgzxzqj
[基金申请] 过年走亲戚时感受到了所开私家车的鄙视链 +3 瞬息宇宙 2026-02-15 5/250 2026-02-16 14:23 by aspect3000
[硕博家园] 江汉大学解明教授课题组招博士研究生/博士后 +3 cleverlyy 2026-02-12 3/150 2026-02-12 21:02 by qsdf1
信息提示
请填处理意见