 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 775 | 回復: 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 ] |
回復此樓» 猜你喜歡
0703化學 305求調劑
已經有6人回復
-
一志愿211，0703化學310分求調劑
已經有3人回復
-
廣西大學材料導師推薦
已經有5人回復
-
化學調劑
已經有5人回復
-
資源與環(huán)境 調劑申請(333分)
已經有4人回復
-
311求調劑
已經有14人回復
-
269專碩求調劑
已經有4人回復
-
085600材料與化工306
已經有3人回復
-
一志愿西安交通大學材料工程專業(yè) 282分求調劑
已經有12人回復
-
材料求調劑
已經有5人回復
» 本主題相關商家推薦: (我也要在這里推廣)
» 本主題相關價值貼推薦，對您同樣有幫助:
- 分子動力學模擬長程庫侖勢的處理問題
已經有8人回復
- VASP分子動力學優(yōu)化的INCAR如下，怎么CONTCAR和POSCAR沒有什么變化��？
已經有13人回復
- 【討論】大家好，加水后跑完分子動力學，在log中發(fā)現(xiàn)了問題
已經有15人回復
- 【求助】分子動力學模擬尺度問題？
已經有11人回復
- 【求助】關于用分子動力學來做光譜的問題
已經有14人回復
- 【求助】分子動力學的溶劑化分析
已經有3人回復
yahoohoo
鐵桿木蟲 (著名寫手)
- 模擬EPI: 10
- 應助: 0 (幼兒園)
- 貴賓: 1.55
- 金幣: 7632.4
- 散金: 251
- 紅花: 4
- 帖子: 1176
- 在線: 167小時
- 蟲號: 74894
- 注冊: 2005-06-15
- 性別: GG
- 專業(yè): 理論和計算化學
2樓2011-08-23 13:43:20
wifloveice
金蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 1262.3
- 帖子: 12
- 在線: 13.6小時
- 蟲號: 1066005
- 注冊: 2010-07-29
- 專業(yè): 金屬材料的微觀結構
3樓2011-08-23 17:19:26
yahoohoo
鐵桿木蟲 (著名寫手)
- 模擬EPI: 10
- 應助: 0 (幼兒園)
- 貴賓: 1.55
- 金幣: 7632.4
- 散金: 251
- 紅花: 4
- 帖子: 1176
- 在線: 167小時
- 蟲號: 74894
- 注冊: 2005-06-15
- 性別: GG
- 專業(yè): 理論和計算化學
4樓2011-08-23 17:54:37
wifloveice
金蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 1262.3
- 帖子: 12
- 在線: 13.6小時
- 蟲號: 1066005
- 注冊: 2010-07-29
- 專業(yè): 金屬材料的微觀結構
5樓2011-08-23 23:01:18
yahoohoo
鐵桿木蟲 (著名寫手)
- 模擬EPI: 10
- 應助: 0 (幼兒園)
- 貴賓: 1.55
- 金幣: 7632.4
- 散金: 251
- 紅花: 4
- 帖子: 1176
- 在線: 167小時
- 蟲號: 74894
- 注冊: 2005-06-15
- 性別: GG
- 專業(yè): 理論和計算化學
6樓2011-08-26 18:30:48
相關版塊跳轉
我要訂閱樓主 wifloveice 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調劑 +7 | 26研0 2026-03-15 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 0703化學297求調劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
|
[考研] 280求調劑 +10 | 咕嚕曉曉 2026-03-18 | 11/550 |
|
|
|
[考研] 307求調劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
|
[考研] 298求調劑 +4 | 上岸6666@ 2026-03-20 | 4/200 |
|
|
|
[考研] 299求調劑 +4 | 某某某某位 2026-03-21 | 4/200 |
|
|
|
[考研] 南昌大學材料專碩311分求調劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
|
[考研] 機械專碩299求調劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
|
[考研] 299求調劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
|
[考研] 化學求調劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
|
[考研] 324分 085600材料化工求調劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿華中科技大學，080502，354分求調劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
|
[考研]
 271材料工程求調劑
+8
|
.6lL 2026-03-18 | 8/400 |
|
|
|
[考研] 295求調劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
|
[考研] 353求調劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿985，本科211，0817化學工程與技術319求調劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
|
[考研] 收復試調劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 生物學071000 329分求調劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 085601求調劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|