24小时热门版块排行榜

返回列表

cavediger

金虫 (正式写手)

1ST强帖: 9
应助: 6 (幼儿园)
金币: 709
散金: 5
红花: 17
帖子: 545
在线: 230.2小时
虫号: 1255502
注册: 2011-04-05
专业: 当代宗教

[交流] 关于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:  没权限，拿不到这个！

回复此楼