24小时热门版块排行榜

返回列表

【奖励】本帖被评价5次，作者y1ding增加金币 3.8 个

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

y1ding

铁杆木虫 (著名写手)

[资源] VASP5.3.2大内存崩溃的解决方案

当原子数很多时，同样体系，vasp5.2.2或者vasp4.6.38计算不出问题，vasp5.3.2计算会遇到大内存导致节点down或段错误。
原因是vasp5.3.2在第一个离子步之后，当ALGO/=ALL的情况下，会自作主张的重新生成一下WAVE，而在vasp5.2.2里是没有的。
解决办法：
（仅针对vasp5.3.2，其他版本没有测试）
main.F 2636行开始把以下内容中涉及 IF (.NOT.INFO%LONESW) THEN **** ENDIF 的注释掉!
   IF (NSTEP==1) THEN
!       IF (.NOT.INFO%LONESW) THEN
!          CALL ALLOCW(WDES,W_F,WTMP,WTMP)
!          CALL ALLOCW(WDES,W_G,WTMP,WTMP)
!          DEALLOCATE(CHAM, CHF)
!          ALLOCATE(CHAM(WDES%NB_TOT,WDES%NB_TOT,WDES%NKPTS,WDES%ISPIN), &
!             CHF (WDES%NB_TOT,WDES%NB_TOT,WDES%NKPTS,WDES%ISPIN))
!!          ! setup for scf subspace rotation
!!          INFO%LONESW=.TRUE.
!!          CALL SETUP_SUBROT_SCF(INFO,WDES,LATT_CUR,GRID,GRIDC,GRID_SOFT,SOFT_TO_C,IO%IU0,IO%IU5,IO%IU6)
!!          INFO%LONESW=.FALSE.
!       ENDIF
      CALL PEAD_ELMIN( &
      HAMILTONIAN,KINEDEN, &
      P,WDES,NONLR_S,NONL_S,W,W_F,W_G,LATT_CUR,LATT_INI, &
      T_INFO,DYN,INFO,IO,MIX,KPOINTS,SYMM,GRID,GRID_SOFT, &
      GRIDC,GRIDB,GRIDUS,C_TO_US,B_TO_C,SOFT_TO_C,E, &
      CHTOT,CHTOTL,DENCOR,CVTOT,CSTRF, &
      CDIJ,CQIJ,CRHODE,N_MIX_PAW,RHOLM,RHOLM_LAST, &
      CHDEN,SV,DOS,DOSI,CHF,CHAM,ECONV, &
      NSTEP,LMDIM,IRDMAX,NEDOS, &
      TOTEN,EFERMI,LDIMP,LMDIMP)
!       IF (.NOT.INFO%LONESW) THEN
!          CALL DEALLOCW(W_F)
!          CALL DEALLOCW(W_G)
!          DEALLOCATE(CHAM, CHF)
!          ALLOCATE(CHAM(1,1,1,1),CHF (1,1,1,1))
!       ENDIF
   ENDIF

PEAD_ELMIN 说明：
! This subroutine is called from main.F and takes care of one or all
! of the following:
!
!  ) calculate the macroscopic dipole moment in the cell
!  ) calculate the SCF and/or non-SCF response to an electric
! field specified through EFIELD_PEAD in the INCAR file
!  ) calculate the dielectric tensor excluding and/or
! including local field effects

以上注释后，仅会影响如上量的计算，其它量不影响。

补充说明：
vasp5.3.2这个问题，用以往网上说的
在FFLAGS 选项的赋值中增加  -heap-arrays  64
运行的时候加 ulimit -s unlimited
换用pgi+acml
都是无法解决的。

[ Last edited by y1ding on 2012-11-30 at 21:17 ]

回复此楼