# RNG = rand48 # L = 50 d = 3 N = 125000 # Boundary conditions: periodic. # T = 6.8000000 beta = 0.1470588 alpha0 = 0.0000000 # GRU algorithm for Metropolis sweeps. # Interactions: none. # Initial permutation: identity. # Site selection for Metropolis sweeps: sequential. # Cycle reversals allowed. # Thermalization is detected by 40 turning points of system energy # smoothed over a 200-point window. # Terminate after 1000 Metropolis accumulations. # 1 Metropolis sweep per accumulation. # # Initial HDV = 0.0000000 0.0000000 0.0000000 # Sweep 1 HDV = 12648.0000000 12648.0000000 0.0000000 # Thermalization complete: sweep 880 H = 75092.4000000 # Accumulation complete: acc 1000 H = 77448.6000000 # rho_infty max dev = 0.9934160 at k = 476 # k= 1 88679.2880000 # k= 2 7164.0590000 # k= 3 1505.8680000 # k= 4 697.3290000 # k= 5 357.6330000 # k= 6 216.2560000 # k= 7 138.3830000 # k= 8 95.6870000 # k= 9 69.3760000 # k= 10 52.2870000 # k*counts sum = 125000.0000000 # rho_infty_T = 0.0065840 # # = 77220.0656000 # = 77220.0656000 # = 0.0000000 # = 0.6177605 # = 0.6177605 # = 0.0000000 # # = 7164.0590000 # = 504.5260000 # = 0.0000000 # P[Wx = 0 ] = 1000 / 1000 = 1.00000000 # P[Wx = 1 ] = 0 / 1000 = 0.00000000 # P[Wx = 2 ] = 0 / 1000 = 0.00000000 # P[Wx = 3 ] = 0 / 1000 = 0.00000000 # P[Wx = 4 ] = 0 / 1000 = 0.00000000 # P[Wx = 5 ] = 0 / 1000 = 0.00000000 # P[Wx = 6 ] = 0 / 1000 = 0.00000000 # P[Wx = 7 ] = 0 / 1000 = 0.00000000 # P[Wx = 8 ] = 0 / 1000 = 0.00000000 # P[Wx = 9 ] = 0 / 1000 = 0.00000000 # P[Wx = 10 ] = 0 / 1000 = 0.00000000 # = 0.0000000 # bulkflow box length = 5 # = 1.0152400 # = 10.9300000 # = 34.4850000 # s_ellx = 51.1802747 # s_elly = 157.3769783 # cov_ellx_elly = -319.9970500 # corr_ellx_elly = -0.0397285 # 0.2791324 # s_jumplen 0.4827889 # Site (0,0,0) mean jump length 0.2791324 stddev 0.4827889 # Site (0,0,25) mean jump length 0.2424510 stddev 0.4541118 # Site (0,25,25) mean jump length 0.1632843 stddev 0.3851470 # Site (25,25,25) mean jump length 0.4044722 stddev 0.5499111 # # cos_angle0 sites 0,0,0 and 0,0,1 # cos_angle1 sites 0,0,0 and 0,1,1 # cos_angle2 sites 0,0,0 and 1,1,1 # cos_angle3 sites 0,0,0 and 25,25,25 # = -0.0170000 # = 0.0010000 # = -0.0100000 # = -0.0100000 # s_dot0 = 0.1695280 # s_dot1 = 0.2430186 # s_dot2 = 0.2047973 # s_dot3 = 0.1730027 # # = -0.1961643 # = 0.0278105 # = -0.0832883 # = -0.1109689 # s_cos_angle0 = 0.5371512 # s_cos_angle1 = 0.6154028 # s_cos_angle2 = 0.6000366 # s_cos_angle3 = 0.5532922 # cos_angle0 num defined = 73/1000 = 0.0730000 # cos_angle1 num defined = 75/1000 = 0.0750000 # cos_angle2 num defined = 85/1000 = 0.0850000 # cos_angle3 num defined = 69/1000 = 0.0690000 # # = -0.1220000 # s_kinkdot = 0.4597035 # kink cos_angle num defined = 258/1000 = 0.2580000 # # GRU quantities: # alpha_GRU = 1.7000 # /|X| = 0.0040362 # phi(alpha) = 0.0065840 # GRU quot = 0.6130313 # # Metropolis selves: 0/125000000 ( 0.000%) # Metropolis keeps: 111052620/125000000 ( 88.842%) # Metropolis changes: 13947380/125000000 ( 11.158%) # Elapsed seconds: 121.700865