B zd!@sddlZddlZddlZddlmZddlmZddlm Z ddgZ ej ddd Gd dde Z ej dd Gd dde ZdS) N) instance_of) RandomState)utilEnsembleEnsembleConfigurationT)slotsfrozenc@s,eZdZeZeZeZeZdS)rN) __name__ __module__ __qualname__attribadaptation_lagadaptation_time scale_factor evaluatorrr]/work/yifan.wang/ringdown/master-ringdown-env/lib/python3.7/site-packages/ptemcee/ensemble.pyr s)rc@sneZdZdZejeeedZeje j e j dZ eje j e jdZeje j ddZeje j ddZejeeddZejeeeedZejed Zejedd d Zejedd Zejedd Zejedd Zeje j ddd Zeje j ddd Zeje j ddd Zeje j ddd Z ej!d dZ"e j!ddZ#ddZ$ddZ%ddZ&ddZ'dddZ(ddZ)ddZ*dS) rz This contains as little contextual information as it can. It represents an ensemble.py that performs steps in the parameter space. )type validator)r converterN)rdefaultF)rrr)rrfactory)rr)rinitr)rrcCst|std|jdS)Nz{} must be callable.)callable ValueErrorformatname)self attributevaluerrr _is_callable5szEnsemble._is_callablecCst|t|jkrtddS)Nz>Number of temperatures not consistent with starting positions.)lenxr)rr r!rrr_is_consistent:szEnsemble._is_consistentcCs|jj\|_|_|_t|j|jf|_t|jd|j|_ |j dksR|j dkrx| |j\}}| |||_ ||_ |j tj krtddS)Nrz;Attempting to start with samples outside posterior support.)r$shapentempsnwalkersndimnpZonesjumps_proposedfullswaps_proposedlogPlogl _evaluate_tempered_likelihoodinfanyr)rr/Zlogprrr__attrs_post_init__?szEnsemble.__attrs_post_init__cCs||j|j|j||j|j|j|_|j|j}|jr|jdkr| |j |j |}|j |7_ |j|j |j|d7_|j d7_ dS)Nr)betas) _stretchr$r.r/_temperature_swapsswaps_acceptedr-adaptiver'_get_ladder_adjustmenttimer5r1)rratiosZdbetasrrrstepNs z Ensemble.stepc Csdt|j|jf|_|jd}|j}|j}t|jj}x$dD]}|} |dd} |dd| ddddf} |dd| ddddf} t |j j | |||fd} t |||f}xt |D]|}|j jd||d}| ||ddf| |ddf|df| |ddddf| ||ddf||ddddf<qW||\}}|||}|t| ||dd| ddf}t|j j dd||fd}||k}|}|d|f|ddf| d|f|ddf<|d ||dd| ddfd |<|d ||dd| ddfd |<|||f|jdd| ddf<q@WdS) zI Perform the stretch-move proposal on each ensemble.py. )rrrN)lowhighsizer)r@rA)rB)r*Zzerosr'r(jumps_acceptedr)log_configrexp_randomuniformemptyrangerandintreshaper0r1flatten)rr$r.r/wdtZlogajZj_updateZj_sampleZx_updateZx_samplezykjsZy_loglZy_logpZy_logPZ logp_acceptZlogrZacceptsrrrr6]s6  B(0**zEnsemble._stretchcCsj|jdd}|d|jf}t|}tj||jj |}t |t d||d}t t |dS)zi Evaluate the log likelihood and log prior functions at the specified walker positions. NrBr>)r>)r&rLr)r# itertoolschain from_iterable_mapperrErr*ZfromiterfloattupleZrollaxis)rr$r&valueslengthresultsarrayrrrr0s zEnsemble._evaluatec CsP|dkr|j}tjdd||dddf}WdQRXtj |t|<|S)a Compute tempered log likelihood. This is usually a mundane multiplication, except for the special case where beta == 0 *and* we're outside the likelihood support. Here, we find a singularity that demands more careful attention; we allow the likelihood to dominate the temperature, since wandering outside the likelihood support causes a discontinuity. Nignore)invalid)r5r*Zerrstater2isnan)rr/r5ZloglTrrrr1s zEnsemble._tempered_likelihoodc Cs|}|jj||jj}||jj}||dd|dd}td|dd}|t|9}dt|d|d|dd<||S)z Execute temperature adjustment according to dynamics outlined in `arXiv:1501.05823 `_. NrBrr)copyrErrr*diffrFZcumsum) rr;Zbetas0r<r5ZdecaykappaZdSsZdeltaTsrrrr:s "zEnsemble._get_ladder_adjustmentcCs|j}t|j}t|d|_xt|dddD]}|j|}|j|d}||} |j|} |j|} t |jj |d} | ||| f||d| f} | | k}t ||j|d<t ||| |ddf}t ||| |f}t ||| |f}||d| |ddf||| |ddf<||d| |f||| |f<||d| |f| ||d| |f||| |f<|||d| |ddf<|||d| |f<|| |||d| |f<q4W|S)z| Perform parallel-tempering temperature swaps on the state in ``x`` with associated ``logP`` and ``logl``. rrrB)rAN) r(r#r5r*rIr8rJrGZ permutationrDrHsumrc)rr$r.r/r(r'iZbiZbi1ZdbetaZipermZi1permZracceptZpacceptselZx_tempZ logl_tempZ logP_temprrrr7s.     , 8 zEnsemble._temperature_swaps)N)+r r r __doc__r rrrrEr*ZndarrayrZ_ladderr5r_r$r.r/boolr9rrGmaprYintr;r(r'r)r+rCr-r8rr"r%r4r=r6r0r1r:r7rrrrrs4 * )r numpyr*rVZattr.validatorsrZnumpy.random.mtrandrr__all__sobjectrrrrrrs