| 查看: 765 | 回复: 5 | ||
wifloveice金虫 (初入文坛)
|
[求助]
分子動力學以OPENMP平行化時編譯問題
|
|
С???????4??????X????IVF???g????????OPENMP??г?????K??Lennard-Jones ????a???M?Мy??? ??????????Y????OPENMP?????\?????????c?????????????У? ??????????? ???P????? CPU??????m?_??100%?? ?????????(???????p)???????Y??n(?c???_??????\??r?Y?????)?? ????? ??????gOPENMP?a???????N?^?? ?? ??????????`???x?x M9832043@mail.npust.edu.tw SUBROUTINE SSEVAL use omp_lib IMPLICIT NONE INTEGER(KIND=4) IN,LIST,NPOINT,I,J,JBEGIN,JEND,NATOM,NABORS,KSORT,NSTEP,NATOM1 REAL(KIND=8) RL,X,Y,Z,XI,YI,ZI,RIJ ,ENERGY,VIRIAL LOGICAL LG,LUPDAT COMMON /INTGRS/ IN(10) COMMON /LGCLS / LG(4) COMMON /NABLST/ LIST(45000) ,NPOINT(360) COMMON /POSIT / X(360),Y(360),Z(360) COMMON /REALS / RL(26) COMMON /ZLOOP / XI,YI,ZI,RIJ ,I,J,JBEGIN,JEND EQUIVALENCE (IN(2),KSORT) ,(IN(3),NABORS) ,(IN(4),NATOM) EQUIVALENCE (IN(5),NATOM1),(IN(6),NSTEP) EQUIVALENCE (RL(6),ENERGY),(RL(24),VIRIAL),(LG(4),LUPDAT) ENERGY=0.D0 VIRIAL=0.D0 NABORS=0 LUPDAT=.FALSE. IF (MOD(NSTEP,KSORT) .EQ. 0) LUPDAT=.TRUE. !$OMP PARALLEL PRIVATE(I) !$OMP DO DO 300 I=1,NATOM1 JBEGIN=NPOINT(I) JEND =NPOINT(I+1)-1 IF (LUPDAT) THEN NPOINT(I)=NABORS+1 JBEGIN =I+1 JEND =NATOM ENDIF XI=X(I) YI=Y(I) ZI=Z(I) CALL SSINLP 300 CONTINUE !$OMP END DO !$OMP END PARALLEL IF (LUPDAT) NPOINT(NATOM)=NABORS+1 RETURN END *************************************************************** SUBROUTINE SSINLP IMPLICIT NONE INTEGER(KIND=4) IN,LIST,NPOINT,I,J,JBEGIN,JEND,JX REAL(KIND=8) FX,FY,FZ,Z1,X,Y,Z,RL,XI,YI,ZI,RIJ REAL(KIND=8) CUBE,CUBEH,ENERGY,ESHFT,FSHFT,RCUT,VIRIAL,FOR,ENR,R6INV,RSQINV,RSQ,ZIJ,YIJ,XIJ LOGICAL LG,LSHIFT,LUPDAT COMMON /FORCE / FX(360) ,FY(360) ,FZ(360) COMMON /LGCLS / LG(4) COMMON /NABLST/LIST(45000),NPOINT(360) COMMON /POSIT /X(360),Y(360),Z(360) COMMON /REALS /RL(26) COMMON /ZLOOP /XI,YI,ZI,RIJ,I,J,JBEGIN,JEND EQUIVALENCE (RL(2),CUBE), (RL(3),CUBEH), (RL(6),ENERGY) EQUIVALENCE (RL(7),ESHFT), (RL(10),FSHFT), (RL(11),RCUT) EQUIVALENCE (RL(24),VIRIAL), (LG(3),LSHIFT), (LG(4),LUPDAT) DO 400 JX=JBEGIN,JEND J=JX IF (.NOT.LUPDAT) J=LIST(JX) XIJ=XI-X(J) YIJ=YI-Y(J) ZIJ=ZI-Z(J) RSQ=XIJ*XIJ+YIJ*YIJ+ZIJ*ZIJ RIJ=DSQRT(RSQ) IF (LUPDAT) CALL SSUPDT IF(RIJ.LE.RCUT) THEN RSQINV=1.D0/RSQ R6INV =RSQINV*RSQINV*RSQINV ENR = 4.D0*R6INV*(R6INV-1.D0) FOR = RSQINV*48.D0*R6INV*(R6INV-0.5D0) IF (LSHIFT) ENR=ENR+ESHFT+RIJ*FSHFT IF (LSHIFT) FOR=FOR-RSQINV*RIJ*FSHFT FX(I)=FX(I)+FOR*XIJ FX(J)=FX(J)-FOR*XIJ FY(I)=FY(I)+FOR*YIJ FY(J)=FY(J)-FOR*YIJ FZ(I)=FZ(I)+FOR*ZIJ FZ(J)=FZ(J)-FOR*ZIJ ENERGY=ENERGY+ENR VIRIAL=VIRIAL-FOR*RSQ ENDIF 400 CONTINUE RETURN END [ Last edited by ghcacj on 2011-8-27 at 08:44 ] |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 分子动力学模拟长程库仑势的处理问题
已经有8人回复
- VASP分子动力学优化的INCAR如下，怎么CONTCAR和POSCAR没有什么变化啊？
已经有13人回复
- 【讨论】大家好，加水后跑完分子动力学，在log中发现了问题
已经有15人回复
- 【求助】分子动力学模拟尺度问题？
已经有11人回复
- 【求助】关于用分子动力学来做光谱的问题
已经有14人回复
- 【求助】分子动力学的溶剂化分析
已经有3人回复
yahoohoo
铁杆木虫 (著名写手)
- 模拟EPI: 10
- 应助: 0 (幼儿园)
- 贵宾: 1.55
- 金币: 7632.4
- 散金: 251
- 红花: 4
- 帖子: 1176
- 在线: 167小时
- 虫号: 74894
- 注册: 2005-06-15
- 性别: GG
- 专业: 理论和计算化学
2楼2011-08-23 13:43:20
wifloveice
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1262.3
- 帖子: 12
- 在线: 13.6小时
- 虫号: 1066005
- 注册: 2010-07-29
- 专业: 金属材料的微观结构
3楼2011-08-23 17:19:26
yahoohoo
铁杆木虫 (著名写手)
- 模拟EPI: 10
- 应助: 0 (幼儿园)
- 贵宾: 1.55
- 金币: 7632.4
- 散金: 251
- 红花: 4
- 帖子: 1176
- 在线: 167小时
- 虫号: 74894
- 注册: 2005-06-15
- 性别: GG
- 专业: 理论和计算化学
4楼2011-08-23 17:54:37
wifloveice
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1262.3
- 帖子: 12
- 在线: 13.6小时
- 虫号: 1066005
- 注册: 2010-07-29
- 专业: 金属材料的微观结构
5楼2011-08-23 23:01:18
yahoohoo
铁杆木虫 (著名写手)
- 模拟EPI: 10
- 应助: 0 (幼儿园)
- 贵宾: 1.55
- 金币: 7632.4
- 散金: 251
- 红花: 4
- 帖子: 1176
- 在线: 167小时
- 虫号: 74894
- 注册: 2005-06-15
- 性别: GG
- 专业: 理论和计算化学
6楼2011-08-26 18:30:48
