24小时热门版块排行榜    

查看: 653  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 选择 +8 459582015 2024-05-31 9/450 2024-06-02 12:47 by mch20120521
[基金申请] 三个评委收到的基金包里面的项目是相同的吗? +3 河马の史诗 2024-06-02 3/150 2024-06-02 12:17 by lyfbangong
[基金申请] 入职高校3年发表10+SCI,尽人事听天命 +31 kaoyan250 2024-05-27 42/2100 2024-06-02 12:09 by ECOblance
[考博] 导师不让硕转博,让我去国外读博,能理解吗? +12 萧山幽谷 2024-05-29 20/1000 2024-06-02 12:01 by yuan0806
[基金申请] 化学B02口青基 代表作都是什么水平的?向大佬求助 +10 arthas_007 2024-06-01 12/600 2024-06-02 11:24 by duxin_30
[教师之家] 职能部门工作人员态度不好是普遍的吗?怎么让他们态度好一些? +4 河西夜郎 2024-06-01 4/200 2024-06-02 09:14 by 摩天思瑞
[教师之家] 博士高校求职 安建大vs西科大 +4 chengmy19 2024-06-01 10/500 2024-06-02 06:54 by icm639
[硕博家园] 每到中夜,情难自抑 +34 sioc-sunj 2024-05-28 58/2900 2024-06-02 02:53 by csyky2007
[硕博家园] 实验室太吵闹,无法安静学习,怎么办? +7 utahh 2024-05-31 12/600 2024-06-02 01:23 by dking22
[教师之家] 在大地上我们只过一生---看完我的阿勒泰上头了好几天,完结那天晚上几乎失眠 +11 瞬息宇宙 2024-05-27 13/650 2024-06-01 22:15 by otani
[基金申请] 化学口B0109(高分子合成),拿青年基金一般需要怎样的文章水平? +13 salmon95 2024-05-30 22/1100 2024-06-01 21:22 by salmon95
[硕博家园] 每天学术时间不能保证,能保证的只有: +10 hahamyid 2024-05-27 10/500 2024-06-01 21:11 by 小小芝麻影
[教师之家] 研究生在毕业答辩时挂了,遗憾 +10 spacr2022 2024-05-27 13/650 2024-06-01 17:45 by keyaner23
[硕博家园] 论大家对6070后普通教授导师的看法 +6 SNaiL1995 2024-05-28 10/500 2024-06-01 16:55 by 篮球易飞
[基金申请] 九部门发文:不得将专利授权数量作为人才评价、项目评审、职称评定、高校评价等的条件 +16 sjtu2012 2024-05-28 19/950 2024-06-01 15:24 by 天外飞去来
[硕博家园] 哈工大硕博招生!博士每月入学! +4 nailooo 2024-05-30 5/250 2024-06-01 06:47 by anevay
[基金申请] 信息学部函评结束了吗? +6 ducan21 2024-05-28 7/350 2024-05-29 12:10 by WORLD0256
[论文投稿] 真急着毕业,CPB主编终审17天了,邮件催稿了两次,就是一点动静没有 5+3 kkkk夏 2024-05-28 6/300 2024-05-29 11:18 by hitsdu
[基金申请] E05青基有几个评审 +4 KYXY123 2024-05-28 4/200 2024-05-28 19:25 by popt2t
[硕博家园] 课题做到一半担心做不出来怎么办 +4 mkpooooo 2024-05-27 7/350 2024-05-27 11:56 by kanmiaolucky
信息提示
请填处理意见