 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 785 | 回復: 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 ] |
回復此樓» 猜你喜歡
309求調(diào)劑
已經(jīng)有8人回復
-
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷
已經(jīng)有8人回復
-
求調(diào)劑
已經(jīng)有5人回復
-
求化學調(diào)劑
已經(jīng)有5人回復
-
085600，材料與化工321分求調(diào)劑
已經(jīng)有10人回復
-
0703 化學 求調(diào)劑，一志愿山東大學 342 分
已經(jīng)有5人回復
-
一志愿南開大學0710生物學359求調(diào)劑
已經(jīng)有3人回復
-
085600，專業(yè)課化工原理，320分求調(diào)劑
已經(jīng)有4人回復
-
材料與化工272求調(diào)劑
已經(jīng)有16人回復
-
290求調(diào)劑
已經(jīng)有3人回復
» 本主題相關商家推薦: (我也要在這里推廣)
» 本主題相關價值貼推薦，對您同樣有幫助:
- 分子動力學模擬長程庫侖勢的處理問題
已經(jīng)有8人回復
- VASP分子動力學優(yōu)化的INCAR如下，怎么CONTCAR和POSCAR沒有什么變化�。�
已經(jīng)有13人回復
- 【討論】大家好，加水后跑完分子動力學，在log中發(fā)現(xiàn)了問題
已經(jīng)有15人回復
- 【求助】分子動力學模擬尺度問題？
已經(jīng)有11人回復
- 【求助】關于用分子動力學來做光譜的問題
已經(jīng)有14人回復
- 【求助】分子動力學的溶劑化分析
已經(jīng)有3人回復
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
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
wifloveice
金蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 1262.3
- 帖子: 12
- 在線: 13.6小時
- 蟲號: 1066005
- 注冊: 2010-07-29
- 專業(yè): 金屬材料的微觀結構
5樓2011-08-23 23:01:18
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +5 | 家佳佳佳佳佳 2026-03-29 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 291求調(diào)劑 +7 | Y-cap 2026-03-29 | 8/400 |
|
|
|
[考研] 322求調(diào)劑：一志愿湖南大學 材料與化工（085600），已過六級。 +4 | XX小鄧 2026-03-29 | 4/200 |
|
|
|
[考研] 一志愿南昌大學324求調(diào)劑 +5 | hanamiko 2026-03-29 | 5/250 |
|
|
|
[考研] 297求調(diào)劑 +11 | 田洪有 2026-03-26 | 11/550 |
|
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
|
[考研] 343求調(diào)劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|
|
[考研] 283求調(diào)劑 +7 | A child 2026-03-28 | 7/350 |
|
|
|
[考研] 085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富 +4 | 來財aa 2026-03-25 | 4/200 |
|
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學-環(huán)境工程專碩-298 +4 | 手機用戶 2026-03-26 | 4/200 |
|
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +10 | A-哆啦Z夢 2026-03-23 | 16/800 |
|
|
|
[考研] 一志愿陜師大生物學071000，298分，求調(diào)劑 +5 | SYA！ 2026-03-23 | 5/250 |
|
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
|
[考研] 一志愿武理085500機械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
|
[考研] 300分，材料，求調(diào)劑，英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|