# RNG = rand48
# L = 200 d = 3 N = 8000000
# 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 = 825628.8000035 825628.8000035   0.0000000
# Thermalization complete:  sweep   1897 H = 4926273.6000124
# Accumulation complete:  acc   1000 H = 4974101.4000245
# rho_infty max dev =   0.9996162 at k = 2780
# k=  1 <counts> 5670642.1760000
# k=  2 <counts> 458075.4760000
# k=  3 <counts> 95858.1520000
# k=  4 <counts> 44288.9170000
# k=  5 <counts> 22636.3350000
# k=  6 <counts> 13686.0540000
# k=  7 <counts> 8763.9710000
# k=  8 <counts> 6063.2450000
# k=  9 <counts> 4383.2200000
# k= 10 <counts> 3294.2990000
# k*counts sum = 8000000.0000000
# rho_infty_T = 0.0003838
#
# <H>    = 4953890.7018188
# <D>    = 4953890.7018188
# <V>    =   0.0000000
# <h>    =   0.6192363
# <d>    =   0.6192363
# <v>    =   0.0000000
#
# <n2>   = 458075.4760000
# <lmax> = 2000.4590000
# <winding #> =   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
# <fS>   =   0.0000000
# bulkflow box length = 5
# <bulkflow> =   0.8088800
# <ellx>         =   3.2370000
# <elly>         =   7.2040000
# s_ellx         =  19.0149520
# s_elly         =  62.3592580
# cov_ellx_elly  = -12.0983480
# corr_ellx_elly =  -0.0102030
# <jumplen>   0.3194902
# s_jumplen   0.4969165
# Site (0,0,0) mean jump length   0.3194902 stddev   0.4969165
# Site (0,0,100) mean jump length   0.2882247 stddev   0.4908427
# Site (0,100,100) mean jump length   0.4418799 stddev   0.5724877
# Site (100,100,100) mean jump length   0.1964975 stddev   0.4223609
#
# 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 100,100,100
# <dot0>             =  -0.0320000
# <dot1>             =  -0.0210000
# <dot2>             =   0.0010000
# <dot3>             =   0.0010000
# s_dot0             =   0.2665465
# s_dot1             =   0.2421108
# s_dot2             =   0.2169010
# s_dot3             =   0.1703761
#
# <cos_angle0>       =  -0.2555150
# <cos_angle1>       =  -0.1501697
# <cos_angle2>       =   0.0058520
# <cos_angle3>       =   0.0272727
# s_cos_angle0       =   0.6228523
# s_cos_angle1       =   0.5586632
# s_cos_angle2       =   0.6547367
# s_cos_angle3       =   0.6695644
# cos_angle0 num defined = 120/1000 =   0.1200000
# cos_angle1 num defined = 110/1000 =   0.1100000
# cos_angle2 num defined = 84/1000 =   0.0840000
# cos_angle3 num defined = 55/1000 =   0.0550000
#
# <kinkdot>          =  -0.1990000
# s_kinkdot          =   0.4707660
# kink cos_angle num defined = 300/1000 =   0.3000000
#
# GRU quantities:
# alpha_GRU  = 1.7000
# <lmax>/|X| = 0.0002501
# phi(alpha) = 0.0003838
# GRU quot   = 0.6515485
#
# Metropolis selves:  0/8000000000 (  0.000%)
# Metropolis keeps:   7105709801/8000000000 ( 88.821%)
# Metropolis changes: 894290199/8000000000 ( 11.179%)
# Elapsed seconds: 10362.771791