In [1]:

    

include("../src/msd.jl")


    

    
Out[1]:





run_lorentz (generic function with 2 methods)

In [34]:

    

@time stats = run_lorentz(1000000, 1000);


    

    




elapsed time: 3726.852819996 seconds (2004015797864 bytes allocated, 43.49% gc time)

In [35]:

    

stats2 = [filter(x->x[6]>10000, stats_object) for stats_object in stats];


    

In [8]:

    

using PyPlot


    

    




INFO: Loading help data...

In [40]:

    

for i in 1:6 #length(stats2)
    
    stats_object = stats[i]

    times = [data[1] for data in stats_object];
    means =[data[2] for data in stats_object]; 
    msds = [data[3] for data in stats_object];
    
    semilogx(times, msds./times, "o-", label="$i")
end

legend()

xlim(10, 300)
#ylim(0, 2)
ylim(0.5, 1.5)


    

    













    
Out[40]:





(0.5,1.5)

In [41]:

    

stats2[end]


    

    
Out[41]:





521-element Array{(Float64,Float64,Float64,Float64,Float64,Int64),1}:
 (5.0,0.03576175212958004,3.5757623419350946,0.058704697428415754,1.3474713246704768,371175)    
 (10.0,0.06598799840753888,10.569006391541562,0.11762380647866656,2.9187923444888,369197)       
 (15.0,0.09813510444972975,19.917042319872966,0.17170620348090992,4.151538911539102,367220)     
 (20.0,0.13204793641592535,31.011381382625054,0.2195295378702337,5.242865018183112,365279)      
 (25.0,0.1668331715559977,43.7657475307284,0.25947687110498663,6.240681420246132,363377)        
 (30.0,0.20344961054917635,58.1017867188667,0.29660709319764333,7.173356259087203,361409)       
 (35.0,0.23380147301785387,73.96405903240198,0.32878609689915383,8.05335463331433,359431)       
 (40.0,0.2664415381681061,91.27939795846133,0.3625571900484847,8.87231588487423,357454)         
 (45.0,0.3016445347494469,110.05592796203737,0.39514677079643323,9.64117301154943,355513)       
 (50.0,0.33778860587088033,130.38419122995037,0.42404598484141504,10.344841219143973,353611)    
 (55.0,0.37525303189498593,152.20079557144862,0.45460224430690416,11.00228315904787,351643)     
 (60.0,0.4078933453158353,175.5293156213118,0.48051384005101305,11.625928520320162,349665)      
 (65.0,0.4442006846768728,200.2742556207671,0.5064192607767224,12.20619804523279,347688)        
 ⋮                                                                                              
 (2550.0,447.3836793363734,2.2955276127109043e6,-0.0885038319240573,-0.45914917752656725,10153) 
 (2555.0,448.4334366821088,2.304724064135806e6,-0.08842322708263049,-0.45934739205554376,10139) 
 (2560.0,449.5403029877312,2.3147675271009663e6,-0.08842507608656278,-0.4602437836059763,10121) 
 (2565.0,451.0394973007621,2.3231661384599726e6,-0.08889129998218315,-0.4599844116231533,10114) 
 (2570.0,452.44847968723343,2.3312974286410473e6,-0.08903086029350495,-0.4591418247081158,10105)
 (2575.0,453.5543425196591,2.34066323937383e6,-0.08897655623986722,-0.4594293015856956,10091)   
 (2580.0,454.6132945451405,2.3499577390404725e6,-0.0888959434156215,-0.4596303997382307,10077)  
 (2585.0,455.72782322606304,2.3601179319136073e6,-0.08889687147509831,-0.4605302749371365,10059)
 (2590.0,457.2438862227771,2.368598618689869e6,-0.08936431124459782,-0.460280912563261,10052)   
 (2595.0,458.61736129245486,2.376587361161933e6,-0.08941009161769403,-0.4592506835804602,10044) 
 (2600.0,459.7332632410062,2.3860483917170907e6,-0.08935423765730106,-0.45953927735759503,10030)
 (2605.0,460.8024118490757,2.3954389404572453e6,-0.08927247062221405,-0.459743973365875,10016)  

In [42]:

    

stats2[1]


    

    
Out[42]:





221-element Array{(Float64,Float64,Float64,Float64,Float64,Int64),1}:
 (5.0,0.0005893829595243692,2.6514053105886295,0.0002461241361517818,1.089017667844641,26557753)    
 (10.0,0.0006609540779410161,6.5367101782713215,0.00030356972595024513,2.1043525159115815,26403832) 
 (15.0,0.0010697864813240697,10.813668897557063,-0.00032254706675602065,2.7819675226503353,26249459)
 (20.0,0.0006174726268212749,15.304498195782775,-0.0007502388896040342,3.3550990643792806,26095475) 
 (25.0,0.0004920213147661538,19.951633934314653,-8.430918498104511e-6,3.890296062610102,25942388)   
 (30.0,0.0010160016080403113,24.722952887589994,0.00030270381695561343,4.418318928785642,25788788)  
 (35.0,0.0011859775396365635,29.597152552709847,0.0007581086311486826,4.932619256524049,25634867)   
 (40.0,0.0015663976351611292,34.560255746787256,0.00017275112886179768,5.435335164366968,25480494)  
 (45.0,0.0010451374017777283,39.59648848317168,-0.0003718666404765007,5.928010231152214,25326510)   
 (50.0,0.0009266747793667578,44.69611480154076,-0.00013920313328020045,6.41441869978115,25173423)   
 (55.0,0.001353545787740194,49.869207198782895,0.00016206141948725924,6.896535346868266,25019823)   
 (60.0,0.001394659448708486,55.10212665616488,0.0010799079147023416,7.368677131219281,24865902)     
 (65.0,0.0018053879153457815,60.37975807686846,0.0010185580629246546,7.832119919904855,24711529)    
 ⋮                                                                                                  
 (1050.0,7.101385351964907,189431.27109755218,0.0007031018638675724,0.7517775299034688,13525)       
 (1055.0,7.611258828370341,194924.99463328026,-0.0009180243188271674,0.6949637648313418,13118)      
 (1060.0,7.850341535120583,200344.3094331561,-0.003281302712669728,0.6391584187734209,12755)        
 (1065.0,8.112127383276954,205591.2803649387,-0.00548444850405795,0.5948502508643054,12408)         
 (1070.0,8.403502415261482,210862.74333861264,-0.007594428897099004,0.5549186797687495,12054)       
 (1075.0,8.256747965464694,215978.11864575942,-0.006118399904004729,0.5126881688405329,11748)       
 (1080.0,8.681256835863032,221316.6028911052,-0.007147092426479946,0.47182473296002714,11445)       
 (1085.0,9.01153485914229,226552.27509768418,-0.009902151375600535,0.43102790944606983,11177)       
 (1090.0,9.198883320408346,232170.04207693166,-0.011508620755651303,0.3921542526016224,10890)       
 (1095.0,9.58200172293014,237405.4340185016,-0.013890539163192848,0.36264632557492904,10615)        
 (1100.0,9.619415279112829,242429.5041443338,-0.013410926112337422,0.33082552188979664,10379)       
 (1105.0,10.038514223662787,247669.28160559092,-0.014304986178829882,0.30010712661189354,10145)     

In [ ]: