love_yy 发表于 2024-11-13 23:48
看了论坛里面也还是没有修复bug.
https://theory.cm.utexas.edu/henkelman/forum/viewtopic.php?t=16161

亲测 5楼有效，http://bbs.keinsci.com/thread-46112-1-1.html这个贴也可以参考乐平老师的方法

看了论坛里面也还是没有修复bug.
https://theory.cm.utexas.edu/hen ... ewtopic.php?t=16161

由于VTST中chain.F移除了SUBROUTINE chain_stress这部分。解决方式：在VASP官方的chain.F中复制过来，加在 END SUBROUTINE chain_force后面。并在55行PUBLIC :: chain_force后面加入chain_stress。

这是官方源文件的子函数：

1.     SUBROUTINE chain_stress( stress )
   
2. 
   
3.        implicit none
   
4.        REAL(q),intent( inout ) :: stress(3,3)
   
5.        REAL(q)                 :: stress_all(3,3,2)
   
6.        INTEGER                 :: node
   
7. 
   
8.        IF ( images == 0 ) RETURN
   
9.        IF ( spring == -1000 ) RETURN
   
10. 
    
11. #if defined(MPI) || defined (MPI_CHAIN)
    
12.        node = comm_chain%node_me
    
13. 
    
14. !======================================================================
    
15. ! Parallel tempering return
    
16. !======================================================================
    
17.        IF ( LTEMPER ) THEN
    
18.           RETURN
    
19. !======================================================================
    
20. ! VCA average stresses over two images
    
21. !======================================================================
    
22.        ELSE IF ( LVCAIMAGES ) THEN
    
23.           stress_all(:,:,1:2)  = 0
    
24.           stress_all(:,:,node) = stress
    
25. #ifndef old_vca
    
26.           ! first core in image merges the data
    
27.           IF ( comm%node_me == 1 ) THEN
    
28.              !! 3x3 stress components and two molecular dynamics streams
    
29.              CALLMPI_C( M_sum_d( comm_chain, stress_all(1,1,1), 3*3*2))
    
30.              stress(:,:)=stress_all(:,:,1)*VCAIMAGES+stress_all(:,:,2)*(1-VCAIMAGES)
    
31.           ELSE
    
32.              stress=0
    
33.           ENDIF
    
34.           ! now sum over all cores in one image
    
35.           CALLMPI_C( M_sum_d( comm, stress( 1, 1 ), 3*3 ) )
    
36. #else
    
37.           !! 3x3 stress components and two molecular dynamics streams
    
38.           CALLMPI_C( M_sum_d( comm_chain, stress_all( 1, 1, 1 ), 3*3*2 ) )
    
39.           stress = 0
    
40.           stress( : , : ) = stress_all( :, :, 1 ) * VCAIMAGES + stress_all( :, :, 2 ) * ( 1 - VCAIMAGES )
    
41.           RETURN
    
42. #endif
    
43.        END IF
    
44. #endif
    
45.     END SUBROUTINE chain_stress
    

复制代码

老师好，我怀疑是系统问题，我再centos7.9中可以正常编译，在rocky8.9中会出现报错。

6.4.2我编译没问题，建议检查一下这里