本帖產(chǎn)生 1 個(gè) 模擬EPI ，點(diǎn)擊這里進(jìn)行查看

flowerstudy

木蟲 (正式寫手)

應(yīng)助: 4 (幼兒園)
金幣: 4398.6
帖子: 854
在線: 137小時(shí)
蟲號(hào): 162807

[交流] 【求助】求高人給我發(fā)一份GEMC的模擬程序

小弟，最近在學(xué)習(xí)MC，覺得GEMC對(duì)我比較有用，可是在網(wǎng)上，我怎么都找不著。
如果有人有的話，能否給我發(fā)一份呢？萬分感謝！
我的郵箱地址：zw198204@163.com

回復(fù)此樓

» 猜你喜歡

323求調(diào)劑已經(jīng)有6人回復(fù)
一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑已經(jīng)有4人回復(fù)
352求調(diào)劑已經(jīng)有3人回復(fù)
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑已經(jīng)有8人回復(fù)
277材料科學(xué)與工程080500求調(diào)劑已經(jīng)有7人回復(fù)
317求調(diào)劑已經(jīng)有18人回復(fù)
293求調(diào)劑已經(jīng)有5人回復(fù)
280分求調(diào)劑一志愿085802 已經(jīng)有7人回復(fù)
0854電子信息求調(diào)劑已經(jīng)有3人回復(fù)
263求調(diào)劑已經(jīng)有4人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

求助sybyl關(guān)于添加形式電荷的問題已經(jīng)有8人回復(fù)
正在做大學(xué)生創(chuàng)新實(shí)踐項(xiàng)目，車身方面，研究方向一直反復(fù)不定，跪求高人指點(diǎn)啊已經(jīng)有4人回復(fù)
苯環(huán)上的烷基取代。求高人指點(diǎn)。謝謝！已經(jīng)有23人回復(fù)
除掉反應(yīng)生成的乙醇，拉動(dòng)反應(yīng)平衡。求高人指點(diǎn) 已經(jīng)有15人回復(fù)
求各位高人給小弟支個(gè)招：我是否需要跨專業(yè)考研？已經(jīng)有24人回復(fù)
【求助】求將c語言改寫為C++，求教高人啊！已經(jīng)有6人回復(fù)
【求助】求高人給分析一下三方相的xrd，幫我算算晶胞參數(shù)啊,自己不在行這個(gè)啊，贊助50 已經(jīng)有8人回復(fù)
【求助】100金幣求助高人一個(gè)程序解方程已經(jīng)有15人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2010-11-28 18:52:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zyj8119

木蟲 (著名寫手)

模擬EPI: 10
應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
帖子: 2936
在線: 1329.4小時(shí)
蟲號(hào): 664177

★ ★
ghcacj(金幣+2):謝謝 2010-11-29 18:13:56
flowerstudy(金幣+1): 2010-11-29 20:38:15

引用回帖:

Originally posted by flowerstudy at 2010-11-28 18:52:27:
小弟，最近在學(xué)習(xí)MC，覺得GEMC對(duì)我比較有用，可是在網(wǎng)上，我怎么都找不著。
如果有人有的話，能否給我發(fā)一份呢？萬分感謝！
我的郵箱地址：zw198204@163.com

這個(gè)吉布斯系綜蒙特卡羅模擬的程序，比較難找，很多組都是組內(nèi)編寫，不外傳。

贊一下(1人)

回復(fù)此樓

2樓2010-11-29 00:01:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

qphll

金蟲 (正式寫手)

模擬EPI: 10
應(yīng)助: 18 (小學(xué)生)
金幣: 2708.7
帖子: 554
在線: 227.1小時(shí)
蟲號(hào): 89654

★ ★
ghcacj(金幣+2):謝謝 2010-11-29 18:14:05
flowerstudy(金幣+1): 2010-11-29 20:38:20

引用回帖:

沒有人會(huì)這樣給你的, 就算我有的話.

自己寫吧, 有啥困惑的, 拿出來討論.

研究這個(gè)路, 真沒有捷徑.

贊一下(1人)

回復(fù)此樓

3樓2010-11-29 09:50:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zyj8119

木蟲 (著名寫手)

模擬EPI: 10
應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
帖子: 2936
在線: 1329.4小時(shí)
蟲號(hào): 664177

引用回帖:

Originally posted by qphll at 2010-11-29 09:50:07:

沒有人會(huì)這樣給你的, 就算我有的話.

自己寫吧, 有啥困惑的, 拿出來討論.

研究這個(gè)路, 真沒有捷徑.

不過如果你知道GEMC的大概思路，可以幫人家寫個(gè)偽代碼，是吧？高人？

贊一下

回復(fù)此樓

4樓2010-12-02 10:47:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

qphll

金蟲 (正式寫手)

模擬EPI: 10
應(yīng)助: 18 (小學(xué)生)
金幣: 2708.7
帖子: 554
在線: 227.1小時(shí)
蟲號(hào): 89654

★ ★ ★
ghcacj(金幣+3):謝謝 2010-12-02 13:46:43

引用回帖:

Originally posted by zyj8119 at 2010-12-02 10:47:45:

不過如果你知道GEMC的大概思路，可以幫人家寫個(gè)偽代碼，是吧？高人？

我也就只有姚明的一半身高, 實(shí)在很低.

樓主你對(duì)于Gibbis MC的了解, 有哪些? 還有, 你打算用它來做什么? 目前, 你打算怎樣開始自己寫? 有哪些參考文獻(xiàn), 或者書?

我這個(gè)周末前看看能不能將偽碼弄出來...

如果食言, 莫怪我, 哈哈.

贊一下(1人)

回復(fù)此樓

5樓2010-12-02 12:43:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zyj8119

木蟲 (著名寫手)

模擬EPI: 10
應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
帖子: 2936
在線: 1329.4小時(shí)
蟲號(hào): 664177

引用回帖:

Originally posted by qphll at 2010-12-02 12:43:57:

我也就只有姚明的一半身高, 實(shí)在很低.

樓主你對(duì)于Gibbis MC的了解, 有哪些? 還有, 你打算用它來做什么? 目前, 你打算怎樣開始自己寫? 有哪些參考文獻(xiàn), 或者書?

我這個(gè)周末前看看能不能將偽碼弄出來... ...

偽代碼是可以翻譯成任何程序語言的。

贊一下

回復(fù)此樓

6樓2010-12-03 10:52:50

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

qphll

金蟲 (正式寫手)

模擬EPI: 10
應(yīng)助: 18 (小學(xué)生)
金幣: 2708.7
帖子: 554
在線: 227.1小時(shí)
蟲號(hào): 89654

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zh1987hs(金幣+10):謝謝 2010-12-04 23:13:56
ghcacj(金幣+3, 模擬EPI+1):精彩答疑 2010-12-06 12:16:56

答應(yīng)過的事情, 雖然覺得其實(shí)對(duì)樓主幫助也不大. 因?yàn)槲也豢赡芙o你寫一個(gè)完整的GibbsMC, 如果你自己不努力去理解, 那么什么都是白搭.

這是大體結(jié)構(gòu). 按照需要填充就是了. 應(yīng)該此貼不會(huì)再更新了.

文件名: GibbsMC.f90

文件內(nèi)容如下:

PROGRAM MC

  USE global
  IMPLICIT NONE

  ! Local
  INTEGER :: ibox, iblock, istep
  INTEGER :: ts,tf,rate
  DOUBLE PRECISION :: engtst, virtst, engerr

  ! Initialize the random number generator
  CALL seed_random(idum)

  ! Read the input file
  CALL read_input

  ! Get the start time
  CALL system_clock(count=ts, count_rate=rate)

  ! Perform minimization
  CALL relax

  ! Calculate energy and virial
  DO ibox = 1,n_box
CALL eng_total(ibox,energy(ibox),virial(ibox))
CALL write_config(ibox,0)
  ENDDO

  ! Loop over the blocks
  DO iblock = 1,n_blocks

! Reset basic statistics
CALL stats(1,0,0)

! Reset optional statistics
IF(opt .EQ. 1) CALL sample(1,0,0)

! Loop over the steps per block
DO istep = 1,n_steps

   ! Perform a trial move
   CALL trial_move(iblock,istep)

   ! Accumulate basic statistics
   CALL stats(2,iblock,istep)

   ! Accumulate optional statistics
   IF(opt .EQ. 1) call sample(2,iblock,istep)

! End loop over the steps per block
ENDDO

! Average basic statistics for iblock
CALL stats(3,iblock,0)

! Average optional statistics for iblock
IF(opt .EQ. 1) call sample(3,iblock,0)

! Update the max values (disp, volchng, etc)
CALL update_max

! Print status
IF(iblock .LE. n_equil) THEN
   IF(iblock .EQ. 1) WRITE(*,'(A)')'Equilibration stage... '
   IF(MOD(NINT(100.0d0*DBLE(iblock)),10*n_equil) .EQ. 0)&
      &WRITE(*,'(I0,A3)')NINT(100.0d0*DBLE(iblock)/DBLE(n_equil)),'%  '
   IF(iblock .EQ. n_equil) WRITE(*,*)''
ELSE
   IF(iblock .EQ. n_equil + 1)  WRITE(*,'(A)')'Production stage... '
   IF(MOD(NINT(100.0d0*DBLE(iblock-n_equil)),10*(n_blocks-n_equil)).EQ. 0)&
      &WRITE(*,'(I0,A3)')NINT(100.0d0*DBLE(iblock-n_equil)/DBLE(n_blocks-n_equil)),'%  '
   IF(iblock .EQ. n_blocks) WRITE(*,*)''
ENDIF

! Check the total and updated energy consistency for each box (every 10 blocks)
IF(MOD(iblock,1) .EQ. 0) THEN
   DO ibox = 1,n_box
      CALL eng_total(ibox,engtst,virtst)
      engerr = dABS((engtst - energy(ibox)))
      IF (engerr .GE. 1.0d-5) THEN
      WRITE(*,*) 'WARNING: LARGE ERROR BETWEEN THE TOTAL AND UPDATED ENERGY'
      WRITE(*,'(A7,I0,A6,I0,A10,F20.8,A10,F20.8)')&
      &'BLOCK: ',iblock,' BOX: ',ibox,' UPD ENG: ',energy(ibox),' TOT ENG: ',engtst
      ENDIF
   ENDDO
ENDIF

  ! End loop over the blocks
  ENDDO

  ! Get the finish time
  CALL system_clock(count=tf, count_rate=rate)

  ! Write the elapsed time to screen
  CALL time_elapsed(tf,ts,rate)

  ! Write basic statistics
  CALL stats(4,0,0)

  ! Write optional statistics
  IF(opt .EQ. 1) CALL sample(4,n_blocks+1,0)

  ! Write the final configurations
  DO ibox = 1,n_box
CALL write_config(ibox,n_blocks+1)
  ENDDO

END PROGRAM MC

贊一下(2人)

回復(fù)此樓

7樓2010-12-04 23:08:18

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zyj8119

木蟲 (著名寫手)

模擬EPI: 10
應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
帖子: 2936
在線: 1329.4小時(shí)
蟲號(hào): 664177

引用回帖:

Originally posted by qphll at 2010-12-04 23:08:18:
答應(yīng)過的事情, 雖然覺得其實(shí)對(duì)樓主幫助也不大. 因?yàn)槲也豢赡芙o你寫一個(gè)完整的GibbsMC, 如果你自己不努力去理解, 那么什么都是白搭.

這是大體結(jié)構(gòu). 按照需要填充就是了. 應(yīng)該此貼不會(huì)再更新了.

文件名: Gibb ...

高手啊，

回復(fù)此樓

8樓2010-12-05 08:58:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 293求調(diào)劑 +3	濤濤Wjt 2026-03-22	5/250	2026-03-22 22:21 by jiangpengfei
[考研] 招08考數(shù)學(xué) +4	laoshidan 2026-03-20	8/400	2026-03-22 19:56 by 小皮蛋醬
[考研] 一志愿武理材料工程348求調(diào)劑 +5	￣^￣゜汗 2026-03-19	7/350	2026-03-22 19:44 by 公瑾逍遙
[考研] 284求調(diào)劑 +5	Zhao anqi 2026-03-22	5/250	2026-03-22 17:38 by barlinike
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike
[考研] 0703化學(xué)調(diào)劑 +4	妮妮ninicgb 2026-03-21	4/200	2026-03-21 18:39 by 學(xué)員8dgXkO
[考研] 311求調(diào)劑 +3	勇敢的小吳 2026-03-20	3/150	2026-03-21 17:40 by ColorlessPI
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6	77chaselx 2026-03-20	6/300	2026-03-21 07:24 by JourneyLucky
[考研] 265求調(diào)劑 +3	Jack?k?y 2026-03-17	3/150	2026-03-21 03:17 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 材料專碩英一數(shù)二306 +7	z1z2z3879 2026-03-18	7/350	2026-03-20 23:48 by JourneyLucky
[考研] 329求調(diào)劑 +9	想上學(xué)吖吖 2026-03-19	9/450	2026-03-20 22:01 by luoyongfeng
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計(jì)算機(jī)專業(yè)調(diào)劑） +3	123阿標(biāo) 2026-03-17	3/150	2026-03-20 15:58 by 飛行琦
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	12/600	2026-03-20 12:25 by lidianxing
[考研] 求調(diào)劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學(xué)校友
[考研] 0703化學(xué)調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子