B zdn @s`dZddlZddlZddlmZmZddlmZddlm Z m Z ddl m Z Gdd d e Z dS) z,Classes for parallel tempered Markov chains.N)create_bit_generators array2dict)ParallelTemperedChain) ChainData detect_dtypes) BaseSamplerc@seZdZdZdddZdddZedd Zed d Zed d Z eddZ e j ddZ ddZ dddZ eddZeddZdS)ParallelTemperedSamplera\Evolves a collection of parallel tempered Markov chains. Parameters ---------- parameters : tuple or list Names of the parameters to sample. model : object Model object. nchains : int The number of chains to create per temperature. Must be greater than zero. betas : numpy.ndarray of floats The betas (= 1 / temperatures) to use. All betas must be between [0,1]. Must provide at least one. If one is not included in the betas, a warning will be printed. swap_interval : int, optional How often to calculate temperature swaps. Default is 1 (= swap on every iteration). proposals : list, optional List of proposals to use. Any parameters that do not have a proposal provided will use the ``default_propsal``. default_proposal : an epsie.Proposal class, optional The default proposal to use for parameters not in ``proposals``. Default is :py:class:`epsie.proposals.Normal`. default_proposal_args : dict, optional Dictionary of arguments to pass to the default proposal. seed : int, optional Seed for the random number generator. If None provided, will create one. pool : Pool object, optional Specify a process pool to use for parallelization. Default is to use a single core. rNc Csd||_||_|||| | |_| |_t|ttfr:|g}t|tj sPt |}| ||||dS)N) parametersmodelZ set_proposalsseedpool isinstancefloatintnumpyZndarrayarray create_chains) selfr r nchainsbetas swap_interval proposalsadaptive_annealerZdefault_proposalZdefault_proposal_argsr r re/work/yifan.wang/ringdown/master-ringdown-env/lib/python3.7/site-packages/epsie/samplers/ptsampler.py__init__@s  z ParallelTemperedSampler.__init__cs@|dkrtdt|jd}fddt|D_dS)aPCreates a list of :py:class:`chain.ParallelTemperedChain`. Parameters ---------- nchains : int The number of parallel tempered chains to create. betas : array Array of inverse temperatures to use for each parallel tempered chain. swap_interval : int, optional How often to calculate temperature swaps. Default is 1 (= swap on every iteration). adaptive_annealer : object, optional Adaptive anneler adjusting temperatures on the go. Default is `None`. rznchains must be >= 1)r c s8g|]0\}}tjjddjD||dqS)cSsg|]}t|qSr)copydeepcopy).0prrr jszDParallelTemperedSampler.create_chains...)rrrZ bit_generatorZchain_id)rr r r)rZcidbg)rrrrrrr!hsz9ParallelTemperedSampler.create_chains..N) ValueErrorrr enumerateZ_chains)rrrrrZbitgensr)rrrrrrRs z%ParallelTemperedSampler.create_chainscstfddtjDS)zThe betas used for each chain. Shape is (nchains, ntemps). If temperatures are being dynamically adjusted each chain will have a different set of temperatures. csg|]}j|jqSr)chainsr)ri)rrrr!vsz1ParallelTemperedSampler.betas..)rrranger)rr)rrrpszParallelTemperedSampler.betascCs |jjdS)z(Returns the number of temperatures used.r)rshape)rrrrntempsxszParallelTemperedSampler.ntempscCs d|jS)z1The temperatures used. Shape is (nchains, ntemps)g?)r)rrrr temperatures}sz$ParallelTemperedSampler.temperaturescCs |jdjS)z6Returns the swap interval used for parallel tempering.r)r%r)rrrrrsz%ParallelTemperedSampler.swap_intervalcCsx|jD] }||_qWdS)zSets the swap interval to use.N)r%r)rintervalcrrrrs cCsjt|jd|}tt|t||jd}||jx$t |jD]\}}t||||<qDWt |j j S)a:Concatenates dictionary attributes over all of the chains. This is a convenience function used by properties such as ``current_positions`` to gather all of the dictionary attributes from the chains. Parameters ---------- attr : str The name of the attribute to get from the chains. The attribute is assumed to return a dictionary. Returns ------- dict : Dictionary mapping parameters to arrays. The arrays have shape ``ntemps x nchains``. r)Zdtypesr)) getattrr%rlistkeysrr)extendrr$rdataT)rattrdoutiichainrrr_concatenate_dictss  z*ParallelTemperedSampler._concatenate_dictscsJdkr"ttfdd|j}nttfdd|j}tj|ddS)aConcatenates the given attribute over all of the chains. This is a convenience function used by properties such as ``positions`` to gather all of the array attributes from the chains. Parameters ---------- attr : str The name of the attribute to get from the chains. The attribute is assumed to return a (structred) array. item : str or array index, optional Get a particular item from the (structred) array from each chain before concatenating. Returns ------- array : The returned array has shape ``ntemps x nchains x niterations``. Ncs t|S)N)r-)x)r3rrz=ParallelTemperedSampler._concatenate_arrays..cst|S)N)r-)r9)r3itemrrr:r;r)Zaxis)r.mapr%rstack)rr3r<Zarrsr)r3r<r_concatenate_arrayssz+ParallelTemperedSampler._concatenate_arrayscCs|jdkrdS|dS)zThe history of the temperature acceptance from all of the chains. If ntemps is 1, just returns None. .. note:: This does not return a structured array, since there is only one field. rNtemperature_acceptance)r)r?)rrrrr@s z.ParallelTemperedSampler.temperature_acceptancecCs|jdkrdS|dS)zThe history of all the temperature swaps from all of the chains. If ntemps is 1, just returns None. .. note:: This does not return a structured array, since there is only one field. rNtemperature_swaps)r)r?)rrrrrAs z)ParallelTemperedSampler.temperature_swaps)rNNNNNN)rN)N)__name__ __module__ __qualname____doc__rrpropertyrr)r*rsetterr8r?r@rArrrrr s!        r )rErrZepsierrZ epsie.chainrZepsie.chain.chaindatarrbaserr rrrrs