| 1 | 1/1 | 返回列表 |
| 查看: 633 | 回復(fù): 0 | ||
moyi22334新蟲 (初入文壇)
|
[求助]
想固定x軸作優(yōu)化,重新編譯了VASP,為什么沒起作用?
|
|
按照網(wǎng)上的方法,修改了constr_cell_relax.F文件,并重新編譯了VASP,并沒有起作用,可能的原因是什么? 代碼如下: SUBROUTINE CONSTR_CELL_RELAX(FCELL) USE PREC REAL(Q) FCELL(3,3), SAVE(3) LOGICAL FILFLG INTEGER ICELL(3) INQUIRE(FILE='OPTCELL',EXIST=FILFLG) IF (FILFLG) THEN OPEN(67,FILE='OPTCELL',FORM='FORMATTED',STATUS='OLD') READ(67,"(3I1)" (ICELL(I),I=1,3)CLOSE(67) DO I=1,3 SAVE(I)=FCELL(I,I) ENDDO FCELL=0.0d0 DO I=1,3 IF (ICELL(I)==1) FCELL(I,I)=SAVE(I) ENDDO ENDIF ! just one simple example ! relaxation in x directions only ! SAVE=FCELL(1,1) ! FCELL=0 ! F90 style: set the whole array to zero ! FCELL(1,1)=SAVE RETURN END SUBROUTINE 方法來(lái)自這里:http://blog.sciencenet.cn/blog-671981-731203.html |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|