| 查看: 954 | 回复: 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 ] |
回复此楼» 猜你喜欢
综述论文作为代表作会不会影响评审专家的印象分？
已经有12人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有4人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有29人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有9人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
帮忙看看fileCode
已经有13人回复
-
为什么网上很多人说本周 12号出结果
已经有7人回复
-
读博的好处
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 分子动力学模拟长程库仑势的处理问题
已经有8人回复
- VASP分子动力学优化的INCAR如下，怎么CONTCAR和POSCAR没有什么变化啊？
已经有13人回复
- 【讨论】大家好，加水后跑完分子动力学，在log中发现了问题
已经有15人回复
- 【求助】分子动力学模拟尺度问题？
已经有11人回复
- 【求助】关于用分子动力学来做光谱的问题
已经有14人回复
- 【求助】分子动力学的溶剂化分析
已经有3人回复
yahoohoo
铁杆木虫 (著名写手)
- 模拟EPI: 10
- 应助: 0 (幼儿园)
- 贵宾: 1.55
- 金币: 7636.4
- 散金: 251
- 红花: 4
- 帖子: 1179
- 在线: 167.5小时
- 虫号: 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
- 金币: 7636.4
- 散金: 251
- 红花: 4
- 帖子: 1179
- 在线: 167.5小时
- 虫号: 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
- 金币: 7636.4
- 散金: 251
- 红花: 4
- 帖子: 1179
- 在线: 167.5小时
- 虫号: 74894
- 注册: 2005-06-15
- 性别: GG
- 专业: 理论和计算化学
6楼2011-08-26 18:30:48
100