³ò •AâJc@smddkTddkTddkiZdZded„Zd„Zd„Zd„Ze d jo d GHndS( iÿÿÿÿ(t*Ni'gð?c CsI||||ƒ} dG| GH| } |} | } | } dG|GdGHd}g}d}xç||joÙ|| |ƒ\}}||||ƒ} || | ttiƒƒjoJ|i| |fƒ|} | } |d7}| | jo| } |} qn|i| | fƒ|d7}|d|djo|G| G| GHqUqUW|||fS(Ns starting lnp=tdoingslinks of MCMC...ii(tlogtrandomtuniformtappend(tnlinkstdatatparamst prior_infot lnposteriort step_infotproposaltbetat output_periodtlnptoldlnpt oldparamst bestparamstbestlnptnaccepttchaintlinkt newparams((s2/home/dstn/astrometry/projects/holmes/hogg_mcmc.pytmetropolis_chain$s6   !   c Cs§t|ƒ}|\}}}}t|ƒ}x@t|ƒD]2}|tiƒ|||dd…|f7}q7W|d7}||||f} |iƒ||} | | fS(Ni(tlentzerostrangeRRtcopy( RR tnpartwtvtQtnproposetdparamstit newstep_infoR((s2/home/dstn/astrometry/projects/holmes/hogg_mcmc.pytgaussian_proposalHs   0 cCs|\}}t|djƒ}t|ƒ}|tid|ƒ}||cd7<||f}|iƒ}||c||tiƒ7<||fS(Nii(tfindRRtrandintRtnormal(RR R"R!tindxRR$R((s2/home/dstn/astrometry/projects/holmes/hogg_mcmc.pytaxis_aligned_proposalTs    cCs2t|ƒ}|tid|ƒ\}}||fS(Ni(RRR'(RRRtnlinktnewlnpR((s2/home/dstn/astrometry/projects/holmes/hogg_mcmc.pytequal_posterior_proposal`s t__main__s hello world( tnumpytpylabt numpy.randomRt OUTPUTPERIODRR%R*R-t__name__(((s2/home/dstn/astrometry/projects/holmes/hogg_mcmc.pyss   "