Skip to content

Commit 13a45ea

Browse files
jinwoolim8180jinwoolim8180
committed
multiprocessing
1 parent 8696f56 commit 13a45ea

6 files changed

Lines changed: 259 additions & 5 deletions

‎montecarlo2d.py‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -3,7 +3,7 @@
33
import numpy as np
44

55

6-
@numba.jit
6+
@numba.jit(nopython=True)
77
def _metropolis(spin, L, beta):
88
for _ in range(L ** 2):
99
# randomly sample a spin

‎montecarlo3d.py‎

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -3,7 +3,7 @@
33
import numpy as np
44

55

6-
@numba.jit
6+
@numba.jit(nopython=True)
77
def _metropolis3d(spin, L, beta):
88
for _ in range(L ** 3):
99
# randomly sample a spin
@@ -66,7 +66,7 @@ def simulate(self, beta):
6666
M = self._calc_magnetization(spin)
6767
E1 += E
6868
M1 += M
69-
E2 += E ** 2
70-
M2 += M ** 2
69+
E2 += self.n1 * E ** 2
70+
M2 += self.n1 * M ** 2
7171

72-
return self.n1 * E1, self.n1 * M1, (self.n1 * E2 - self.n2 * E1 * E1) * beta ** 2, (self.n1 * M2 - self.n2 * M1 * M1) * beta
72+
return self.n1 * E1, self.n1 * M1, (E2 - self.n2 * E1 * E1) * beta ** 2, (M2 - self.n2 * M1 * M1) * beta

‎result/plot_L30_D2_EQ1000_MC1000.png‎

81.1 KB
Loading

‎result/plot_L30_D3_EQ1000_MC1000.png‎

76.1 KB
Loading

‎result/result_L30_D2_EQ1000_MC1000.csv‎

Lines changed: 127 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,127 @@
1+
T,E,M,C,X
2+
1.5,-0.9763666666666666,0.987051111111111,0.04476399999996122,0.02346983999996155
3+
1.54,-0.9720466666666666,0.9844088888888888,0.05133291917317823,0.0299617431456776
4+
1.58,-0.96766,0.9819111111111111,0.05280704837186684,0.03224303797472928
5+
1.62,-0.9603088888888889,0.9771533333333333,0.0666864451557188,0.04716608779147093
6+
1.66,-0.955771111111111,0.9744911111111111,0.08236575377814709,0.05901354484605516
7+
1.7,-0.9496266666666666,0.9705977777777777,0.07604117800845221,0.058712781699341576
8+
1.74,-0.9440288888888888,0.9665022222222222,0.0868161521410031,0.07453518263085651
9+
1.78,-0.9340555555555555,0.9601488888888888,0.08779106874834686,0.07887768039952615
10+
1.82,-0.9242666666666666,0.9529733333333332,0.12307370285393907,0.13665874481065965
11+
1.8599999999999999,-0.9165288888888888,0.9480422222222221,0.1187268354723301,0.12525648506568596
12+
1.9,-0.9034288888888888,0.9370777777777777,0.12971394398273153,0.18369070175433413
13+
1.94,-0.8881977777777778,0.9213555555555555,0.16435833185719353,0.42240985108816825
14+
1.98,-0.8824555555555555,0.9215377777777777,0.15797251414253555,0.23914508641974183
15+
2.02,-0.8637466666666667,0.9021444444444444,0.17017220969620123,0.420056710671059
16+
2.06,-0.848991111111111,0.8907377777777777,0.19762872613399649,0.43744935490828485
17+
2.1,-0.83632,0.8779333333333332,0.20344637339379926,0.5164721693121183
18+
2.14,-0.8235133333333333,0.8689955555555555,0.2074256257218145,0.5355050799583946
19+
2.18,-0.7890799999999999,0.8161777777777777,0.2801864826193185,1.5387392456676863
20+
2.2199999999999998,-0.7598199999999999,0.7789066666666666,0.27695482780077846,1.7843001721721543
21+
2.26,-0.7369355555555555,0.7360422222222222,0.3979127913958329,4.499274412979348
22+
2.3,-0.6902266666666667,0.5823177777777777,0.5745932333543429,24.31520222415459
23+
2.34,-0.6405733333333333,0.02345111111111111,0.3086041801608621,61.56387310731243
24+
2.38,-0.6339355555555555,0.4502422222222222,0.4327861018917395,17.749538971055088
25+
2.42,-0.5931488888888888,0.23945555555555553,0.28329585183161377,19.11895863177227
26+
2.46,-0.5787288888888888,0.09983777777777778,0.2980154493577503,43.454817833785
27+
2.5,-0.5645377777777777,0.11395333333333332,0.21933243448888792,22.69257094933333
28+
2.54,-0.5274777777777777,0.06235111111111111,0.20013486749195425,14.096543545056866
29+
2.58,-0.5160133333333333,0.006684444444444444,0.1642807016138177,12.74842366236003
30+
2.62,-0.5037088888888889,0.004686666666666666,0.16798518540100524,13.075954400339267
31+
2.66,-0.4865,0.021506666666666664,0.14670722859781304,7.906182649958228
32+
2.7,-0.47946666666666665,0.055431111111111105,0.15119304983996187,6.458259858436214
33+
2.74,-0.46950888888888886,0.01296,0.11782499452406561,6.962160220600161
34+
2.7800000000000002,-0.4576733333333333,0.040897777777777776,0.12580369718613216,4.4248837474020775
35+
2.8200000000000003,-0.4451733333333333,0.04827555555555555,0.10276725181496661,4.961623048069345
36+
2.8600000000000003,-0.44162888888888885,0.03795777777777778,0.1266129346178305,5.269175507381506
37+
2.9000000000000004,-0.43448666666666663,0.024611111111111108,0.12831266428854557,4.4691607279693475
38+
2.94,-0.42046666666666666,0.001071111111111111,0.10757549580678884,5.029988928193499
39+
2.98,-0.4113222222222222,0.010484444444444444,0.09995144112627327,4.388179680835198
40+
3.02,-0.40340666666666664,0.014224444444444444,0.09997656584262707,4.83326004267844
41+
3.06,-0.39650444444444444,0.01741111111111111,0.0946229944798072,3.5471863108206243
42+
3.1,-0.39226666666666665,0.0005333333333333333,0.08785367094462025,3.226584086021505
43+
3.14,-0.38335777777777774,0.0021777777777777776,0.09139242926240111,2.306510403397027
44+
3.1799999999999997,-0.3755711111111111,0.007146666666666667,0.0796280698284612,3.0116916058700207
45+
3.2199999999999998,-0.37030888888888885,0.002331111111111111,0.07096977525729793,2.9086799158040026
46+
3.26,-0.3635133333333333,0.005313333333333333,0.06932700975656694,2.109513454669393
47+
3.3,-0.3552755555555555,0.010195555555555555,0.07555152494643338,2.156504048484848
48+
3.34,-0.35169999999999996,0.017055555555555553,0.06769406297186092,1.7691754823685961
49+
3.38,-0.34721555555555556,0.019157777777777777,0.061496442701585025,1.9866842327416172
50+
3.42,-0.34043777777777773,0.0010088888888888888,0.06509132572605486,1.8036385756985052
51+
3.46,-0.33295555555555556,0.014093333333333333,0.06348027368475652,1.6725356917148362
52+
3.5,-0.32984444444444444,0.008546666666666666,0.060594321995463864,1.6506873447619046
53+
3.54,-0.3267822222222222,0.00278,0.0605237234936743,2.074832390458255
54+
3.58,-0.3221177777777778,0.02558,0.055540918406624146,1.86962616883923
55+
3.62,-0.31558888888888886,0.005286666666666666,0.05877204094299128,1.6400925217925109
56+
3.66,-0.3110711111111111,0.007462222222222222,0.053756847522072534,1.4507636138433515
57+
3.7,-0.3073244444444444,0.0009888888888888888,0.05477814592971389,1.2925633333333333
58+
3.74,-0.3023488888888889,0.00045111111111111106,0.04921366893216444,1.4285131206179438
59+
3.7800000000000002,-0.29788,0.00548,0.04887939680673387,1.126223214579659
60+
3.82,-0.29323777777777776,0.0013355555555555555,0.04651406089623512,1.2739334485165792
61+
3.86,-0.2916022222222222,0.006122222222222222,0.04714463893139557,1.2862532815198617
62+
3.9,-0.2864311111111111,0.015826666666666666,0.04601088874278608,1.1933346051282052
63+
3.94,-0.2821022222222222,0.00258,0.04467551455704674,1.1856255826283133
64+
3.98,-0.2786333333333333,0.002862222222222222,0.04408440162846524,1.1016382534896703
65+
4.02,-0.27635333333333334,0.005844444444444444,0.04142474993743253,0.9277491431730239
66+
4.0600000000000005,-0.27312,0.011813333333333332,0.03895340661829409,1.0139071089217295
67+
4.1,-0.27246444444444445,0.009739999999999999,0.04174185180778621,0.9124729658536586
68+
4.140000000000001,-0.2691311111111111,0.001631111111111111,0.03794594428496972,0.9575558174986578
69+
4.18,-0.26464,0.008195555555555555,0.04082279760587468,0.9461229728867624
70+
4.220000000000001,-0.2616222222222222,0.0022533333333333333,0.03872856504670657,0.8784410784623485
71+
4.26,-0.2584022222222222,0.000671111111111111,0.036938104824390444,0.8053738090766823
72+
4.300000000000001,-0.25456888888888884,0.007642222222222222,0.03414430959677895,0.7990659214470281
73+
4.34,-0.25210666666666665,0.001751111111111111,0.03635557632001772,0.8902847475678443
74+
4.38,-0.24946666666666664,1.3333333333333333e-05,0.033443005775526044,0.7317604911212583
75+
4.42,-0.24699333333333331,0.006944444444444444,0.03336618164611262,0.7692884615384615
76+
4.46,-0.24413555555555555,0.001902222222222222,0.032845953735915076,0.7368164065769804
77+
4.5,-0.24084888888888886,0.00046666666666666666,0.03447003039780521,0.7258902716049382
78+
4.54,-0.23981777777777777,0.0005977777777777777,0.029117859414655525,0.7743696906510034
79+
4.58,-0.23719555555555555,0.0023377777777777776,0.03180703732660402,0.6386718059194565
80+
4.62,-0.23378888888888888,0.0013155555555555555,0.028994745434139223,0.7037609773929773
81+
4.66,-0.23342888888888888,0.0013822222222222222,0.029824980096233848,0.6836634391988554
82+
4.7,-0.23121333333333333,0.0023866666666666667,0.028429934912730756,0.7146123158392436
83+
4.74,-0.2262511111111111,0.0035377777777777777,0.028326030174809707,0.5596788898265352
84+
4.78,-0.22573555555555555,0.0016866666666666666,0.02666803588483683,0.6563727745234774
85+
4.82,-0.22550888888888887,0.0008133333333333332,0.026530297534669,0.5900507551867218
86+
4.859999999999999,-0.22104666666666664,0.005168888888888889,0.02596547104946708,0.5970532455418381
87+
4.9,-0.21803555555555554,0.00452,0.02489639390994501,0.6031535999999998
88+
4.9399999999999995,-0.21654222222222222,0.004922222222222222,0.025155644431331548,0.5270344354475933
89+
4.98,-0.21344444444444444,0.00646,0.026600198204688265,0.5733764846050868
90+
5.02,-0.2116511111111111,0.006251111111111111,0.02338833263669519,0.6045737326250553
91+
5.0600000000000005,-0.21254666666666666,0.0027133333333333332,0.023307952180335916,0.58459741677646
92+
5.1,-0.20996222222222222,0.004073333333333333,0.023797293519586416,0.5724227546840959
93+
5.140000000000001,-0.2092133333333333,0.0017444444444444443,0.023419764787422095,0.518115672287073
94+
5.18,-0.20400888888888888,0.0003733333333333333,0.02435456703255931,0.5016179974259973
95+
5.220000000000001,-0.20504888888888886,0.0027266666666666663,0.023274274371257557,0.5412153861217538
96+
5.26,-0.20177555555555554,0.0033511111111111107,0.020320695848011654,0.5418554592310942
97+
5.300000000000001,-0.20101777777777777,1.9999999999999998e-05,0.020853009928404546,0.446742909014675
98+
5.34,-0.20196666666666666,0.0019888888888888886,0.022600615491559412,0.4220074698293799
99+
5.38,-0.19727777777777777,0.003324444444444444,0.020199397619036524,0.4647814820322181
100+
5.42,-0.19701777777777776,0.0036066666666666664,0.019676117336819256,0.43550050922509226
101+
5.46,-0.1942488888888889,0.006644444444444444,0.019968402456314566,0.488460643060643
102+
5.5,-0.1907911111111111,0.0014844444444444443,0.018118238971533456,0.43836103111111113
103+
5.54,-0.19022,0.0014488888888888889,0.019344590702342058,0.45723617809867634
104+
5.58,-0.19159777777777776,0.0012,0.01783365221698356,0.4524925527678216
105+
5.62,-0.1876911111111111,0.0002711111111111111,0.017898073310452952,0.4400752202451562
106+
5.66,-0.18654666666666667,0.0005155555555555556,0.01823874883636384,0.4468364161758932
107+
5.7,-0.1837822222222222,0.0001288888888888889,0.018757661639478767,0.39856891695906427
108+
5.74,-0.18344666666666665,0.0025977777777777774,0.017791941547588775,0.4199337506775067
109+
5.78,-0.18243777777777775,0.005502222222222222,0.017520755392190123,0.41288844598231444
110+
5.82,-0.18004666666666666,0.003428888888888889,0.015817343126950826,0.44714177243222597
111+
5.86,-0.17824222222222222,0.002942222222222222,0.017180206460697827,0.43000608570345084
112+
5.9,-0.17878666666666665,0.003013333333333333,0.016689186887548304,0.4127135698681732
113+
5.94,-0.17629333333333333,0.0017577777777777776,0.015750191275518577,0.38154382267115594
114+
5.98,-0.17620666666666665,0.00033999999999999997,0.016927122490551245,0.36834454923820137
115+
6.0200000000000005,-0.17282,0.0002666666666666666,0.016590648251368263,0.33995836101882604
116+
6.0600000000000005,-0.17143555555555554,0.0005622222222222222,0.01601981456671516,0.3962574191419142
117+
6.1000000000000005,-0.1707911111111111,0.0008422222222222222,0.013482787064409203,0.3829692961748633
118+
6.14,-0.16882444444444444,0.002182222222222222,0.01520018001959341,0.3295589402823018
119+
6.18,-0.17061555555555555,0.001791111111111111,0.015377960129356799,0.3631789744696153
120+
6.22,-0.16696444444444444,0.0006288888888888889,0.014463140143069046,0.3445808581636299
121+
6.26,-0.16476666666666664,0.0035199999999999997,0.014836209526596235,0.35235849769258076
122+
6.3,-0.16486666666666666,2.222222222222222e-05,0.01546607317824245,0.33635971781305113
123+
6.34,-0.1641022222222222,0.004415555555555555,0.01334632358876203,0.3410282726954083
124+
6.38,-0.16177999999999998,0.0013466666666666666,0.013297541297746763,0.33170377150818525
125+
6.42,-0.16236444444444442,0.001711111111111111,0.013965520089629865,0.3267671166493596
126+
6.46,-0.16114888888888887,0.000551111111111111,0.013205203517505013,0.33029136291709665
127+
6.5,-0.15991777777777777,0.0003422222222222222,0.011855319500328686,0.3120808779487179

0 commit comments

Comments
 (0)