±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 2206  |  »Ø¸´: 3
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

18772958660

гæ (³õÈëÎÄ̳)

[ÇóÖú] ÔõÑù°ÑÓÃFortran±àдµÄ³ÌÐò¸ÄдΪ²¢ÐгÌÐò ÒÑÓÐ2È˲ÎÓë

ÓÐû¶®Fortran²¢Ðбà³ÌµÄ£¬Ä¿Ç°ÎÒÃæÁÙµÄÒ»¸öÎÊÌ⣺ÓÃFortran ±àдµÄ³ÌÐòÔÚ·þÎñÆ÷ÉÏÅÜÁ˼¸Ìì¾Í²»³ö½á¹û£¬Ïë°ÑËü¸Ä³É²¢ÐеÄÓöàºË¼ÆË㣬µ«Ã»Ñ§¹ý²¢Ðеģ¬ËùÒÔºÜÍ·ÌÛ£¬ÐèҪһλ´óÉñ°ï°ïÎÒ¡£»òÕßÊÇÓÅ»¯Ò»Ï³ÌÐòµÄËã·¨£¬ÏÂÃæÒ»¶Î³ÌÐò£¬Óõ¥ºË¼ÆËãÅܵÄÌ«¾Ã£¬ÔõÑù¸Ä½øÒ»ÏÂËã·¨£º
subroutine initial_condition()
  use global
  implicit none
    real :: r1, p1 = 0.5
        do i = 1,n
                call random_number(r1)
                if(r1 < p1)then
                        s(i) = 1.0
                else
                        s(i) = -1.0
                end if
                sum1 = sum1 + s(i)
        end do
        return
end

subroutine mcs()
  use global
  implicit none
    integer :: ii
    real :: r2, r3, sum2 = 0.0
    real :: ws(n)
        do ii = 1,n
                sum2 = 0.0
                call random_number(r2)
                i = n*r2 + 1
                sum1 = sum1 - s(i)
                do j = 1,n
                        if(ad(i,j) == 1)then         
                                sum2 = sum2 + s(j)
                        end if
                end do
                if(sum2<0)then
                        ws(i) = (1-(1-2*q)*s(i)*(-1.0))/2
                else if(sum2>0)then
                        ws(i) = (1-(1-2*q)*s(i)*(1.0))/2
                else
                        ws(i) = 0.5
                end if
                call random_number(r3)
                if(r3 <= ws(i))then
                        s(i) = -s(i)  
                else
                        s(i) = s(i)
                end if
                sum1 = sum1 + s(i)
        end do
        return
end

!**************************Ising dynamics*******************************
subroutine dynamics()
use global
implicit none
integer :: mm, t = 0
real :: sum3 = 0.0, sum4 = 0.0, sum5 = 0.0
real :: m(time)
do mm = 1,n1
q = q + 0.005
call initial_condition()
do while(t < time)
        call mcs()       
        t = t + 1
        m(t) = abs(sum1/n)
        if(t>200000)then
        sum3 = sum3 + m(t)**4
        sum4 = sum4 + m(t)**2
        sum5 = sum5 + m(t)
        end if
end do
MNq(mm) = sum5/(time-200000)
UNq(mm) = 1 - (sum3/(time-200000))/(3*((sum4/(time-200000))**2))
t = 0
sum1 = 0.0
sum3 = 0.0
sum4 = 0.0
sum5 = 0.0
end do
return
end
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ÌÇÒ¯Ñýºï

гæ (³õÈëÎÄ̳)

4Â¥2018-06-07 09:32:54
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 4 ¸ö»Ø´ð

yongwu22

гæ (³õÈëÎÄ̳)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï
jjdg: ½ð±Ò+1, ¸Ðл²ÎÓë 2017-04-25 18:25:18
Äã´òËãÔÚʲôÑùµÄÓ²¼þ»·¾³ÅܳÌÐò°¡£¬Êǹ²ÏíÄڴ滹ÊÇ·Ö²¼ÄÚ´æ°¡£¿
ÓÐÊÂûʳ£À´Íæ°¡~~
2Â¥2017-04-25 13:47:35
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

yongwu22

гæ (³õÈëÎÄ̳)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

Ö÷º¯ÊýÊÇʲôÑùµÄ°¡£¿Ñ­»·ÄÚÓк¯Êýµ÷Ó㬲¢ÐеÄÄѶȱȽϴ󰡣¬¶øÇÒÐèÒªÖªµÀһЩ±äÁ¿µÄÖµ£¬²ÅÄÜÈ·¶¨º¯ÊýµÄÈȵ㣬Ҳ¾ÍÊdzÌÐòÔËÐÐʱÄÄ×îºÄʱ¼ä
ÓÐÊÂûʳ£À´Íæ°¡~~
3Â¥2017-04-25 13:51:45
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 316Çóµ÷¼Á +10 ÖÛ×Ô¹£ 2026-04-01 11/550 2026-04-02 09:41 by sanrepian
[¿¼ÑÐ] 337Çóµ÷¼Á +10 ¡¶Ê÷¡· 2026-03-29 10/500 2026-04-02 07:51 by ÉÐË®¸óÖ÷
[¿¼ÑÐ] 279Çóµ÷¼Á +7 Īxiao 2026-04-01 7/350 2026-04-01 22:05 by ¿Í¶ûÃÀµÂ
[¿¼ÑÐ] »·¾³¹¤³Ì297·ÖÇóµ÷¼ÁÒ»Ö¾Ô¸º¼¸ßÔº +10 GENJIOW 2026-03-31 11/550 2026-04-01 21:38 by GENJIOW
[¿¼ÑÐ] ÉúÎïÓëÒ½Ò©¿¼Ñе÷¼Á +5 Ìúº©º©123425 2026-03-31 5/250 2026-04-01 18:01 by syh9288
[¿¼ÑÐ] 0703Ò»Ö¾Ô¸ÄÏʦ´ó334Çóµ÷¼Á +4 seven7yu 2026-03-30 4/200 2026-04-01 16:10 by oooqiao
[¿¼ÑÐ] Ó¢Ò»Êý¶þÉúÎïÐÅϢѧ287·Ö£¬±¾¿ÆÉúÎï¿ÆÑ§£¬Çóµ÷¼Á +3 ±ÌË®xyz 2026-03-29 4/200 2026-04-01 15:35 by plum
[¿¼ÑÐ] ²ÄÁÏ0856 Ó¢Ò»Êý¶þ 323 Çóµ÷¼Á +9 Ô¬sy 2026-04-01 9/450 2026-04-01 14:30 by wangjy2002
[¿¼ÑÐ] 349Çóµ÷¼Á +6 zwjjjjjj 2026-03-31 6/300 2026-04-01 09:16 by JourneyLucky
[¿¼ÑÐ] 322Çóµ÷¼Á +4 ìäÙÒXX 2026-03-31 4/200 2026-04-01 08:21 by JourneyLucky
[¿¼ÑÐ] ÕÅ·¼Ãú-Öйúũҵ´óѧ-»·¾³¹¤³Ìר˶-298 +9 ÊÖ»úÓû§ 2026-03-26 9/450 2026-03-31 18:09 by 544594351
[¿¼ÑÐ] 336²ÄÁÏÇóµ÷¼Á +10 ³ÂäÞÓ¨ 2026-03-26 12/600 2026-03-31 17:59 by jp9609
[¿¼ÑÐ] 318Çóµ÷¼Á +10 ³Â³¿79 2026-03-30 10/500 2026-03-31 17:37 by 544594351
[¿¼ÑÐ] 262Çóµ÷¼Á +7 ZZ..000 2026-03-30 8/400 2026-03-31 10:05 by cal0306
[¿¼ÑÐ] 304²ÄÁÏÇóµ÷¼Á +5 ÖÓllll 2026-03-26 5/250 2026-03-30 19:37 by Ô´_2020
[¿¼ÑÐ] 296Çóµ÷¼Á +10 ±Ë°¶t 2026-03-29 10/500 2026-03-30 10:50 by ̽123
[¿¼ÑÐ] 298Çóµ÷¼Á +3 ÖÖÊ¥´Í 2026-03-29 3/150 2026-03-29 12:06 by longlotian
[¿¼ÑÐ] 352·Ö »¯¹¤Óë²ÄÁÏ +5 º£Äɰٴ¨Ly 2026-03-27 5/250 2026-03-28 03:39 by fmesaito
[¿¼ÑÐ] 0856µ÷¼Á +5 ÇóÇóÈÃÎÒÓÐÊé¶Á° 2026-03-26 6/300 2026-03-27 15:12 by caszguilin
[¿¼ÑÐ] µ÷¼Á +4 èÖèÖyoyo 2026-03-26 4/200 2026-03-26 20:43 by fmesaito
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û