以下内容转自:完了,不小心把窗口关了
最近在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: 没权限,拿不到这个! |