24小时热门版块排行榜    

查看: 650  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 博三一直没文章怎么办 +16 133456 2024-05-17 26/1300 2024-05-19 01:20 by 020100099
[基金申请] 信息三处送审了吗? +3 957083516 2024-05-12 3/150 2024-05-18 23:29 by 临川树荫
[基金申请] 大佬们,哪位大神知道E09有没有送啊? +4 牛棕色左右 2024-05-15 7/350 2024-05-18 22:48 by 牛棕色左右
[教师之家] 来用亲身体会,一起说说年轻老师的辛苦 +18 zylfront 2024-05-17 23/1150 2024-05-18 18:25 by zylfront
[基金申请] 连续3年国社科青年本子都没中,今年最后一次青年了,写本子完全浪费时间 +11 ddkk3000 2024-05-15 16/800 2024-05-18 16:13 by kcmn1000
[基金申请] 申请基金代表性成果 +11 lancet0903 2024-05-17 11/550 2024-05-18 13:37 by pinyuan33
[论文投稿] 推荐转投( transfer pending)是否有用? 50+3 lily5289 2024-05-17 6/300 2024-05-18 13:06 by lily5289
[论文投稿] 中文投稿 +7 liujinyu80 2024-05-13 7/350 2024-05-18 11:09 by dgfj564
[教师之家] 博士去高校就是为了有寒暑假吗? +16 wenwen0825 2024-05-16 22/1100 2024-05-18 09:57 by 鱼翔浅底1
[硕博家园] 耐高温垫片求购 +6 Sexyflea 2024-05-16 9/450 2024-05-18 08:31 by nanbeiY
[教师之家] “直接受聘正高专业技术职务”怎么理解 +7 ZHONGWU_U 2024-05-17 8/400 2024-05-17 17:36 by jurkat.1640
[论文投稿] 投稿成功后又想撤回 +5 otani 2024-05-16 5/250 2024-05-17 16:02 by topedit
[基金申请] 怎么成为NSFC评审专家? +12 phamacy 2024-05-12 13/650 2024-05-17 15:40 by oooooo?o
[教师之家] 普通高校有编制讲师与211高校无编制专职科研岗,应该怎么选 +17 1137025804 2024-05-15 22/1100 2024-05-17 12:10 by cleverly
[基金申请] 12个本子最多让给1a2b +16 地球e村长 2024-05-13 19/950 2024-05-16 17:48 by sars518
[基金申请] 青年拔尖出炉 +5 babu2015 2024-05-14 6/300 2024-05-15 21:27 by msjy
[硕博家园] 导师这样说到底有没有机会啊……… +7 理科生yes 2024-05-14 14/700 2024-05-15 18:45 by 理科生yes
[基金申请] 2024工材口的国自然开始涵评了吗? +3 Thanks1989 2024-05-13 3/150 2024-05-14 22:27 by weiyongscut
[论文投稿] 第一次投SCI,一审给了Revise +5 慎独的小花卷 2024-05-13 12/600 2024-05-14 20:58 by 慎独的小花卷
[论文投稿] 七个月了,还在selecting for review +3 g9522 2024-05-14 6/300 2024-05-14 19:00 by lizhengke06
信息提示
请填处理意见