B zd.@sXdZddlmZmZddlZddlmZmZmZddl m Z GdddeZ dd Z dS) z&Base class and functions for samplers.)ABCabstractmethodN) create_seed dump_state load_state)Normalc@seZdZdZdZdZdZdZdZdZ dZ e ddZ e j ddZ e ddZd>dd Ze d d Zej d d Ze d dZej ddZe ddZej ddZe ddZej ddZe ddZe ddZeddZed?ddZe ddZej d dZd!d"Zd#d$Ze d%d&Ze d'd(Ze d)d*Ze d+d,Ze d-d.Z e d/d0Z!e d1d2Z"e d3d4Z#e d5d6Z$d7d8Z%d@d:d;Z&dAdRsz-BaseSampler.set_proposals..)copyrrunionr appendr)r rZdefault_proposalZdefault_proposal_argsZ given_paramsZmissing_paramsr r r set_proposals4s zBaseSampler.set_proposalscCs|jS)zThe model used.)_model)r r r r model[szBaseSampler.modelcCs ||_dS)N)r)r rr r r r`scCs|jS)zList of the chains used.)_chains)r r r r chainsdszBaseSampler.chainscCs ||_dS)N)r)r rr r r riscCs|jS)z+The seed used for the random bit generator.)_seed)r r r r seedmszBaseSampler.seedcCs|dkrt|}||_dS)z:Sets the seed. If the given seed is None, will create one.N)rr )r r!r r r r!rscCs|jS)z(The pool being used for parallelization.)_pool)r r r r poolyszBaseSampler.poolcCs$||_|dkrt}n|j}||_dS)N)r"map_map)r r#Zmap_funcr r r r#~s cCs t|jS)z The number of chains being used.)lenr)r r r r nchainsszBaseSampler.nchainscCs |jdjS)z?The number of iterations the chains have been run for. r)r iteration)r r r r niterationsszBaseSampler.niterationscCsdS)aConcatenates 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. Nr )r attrr r r _concatenate_dictsszBaseSampler._concatenate_dictscCsdS)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 sould have ``nchains x niterations`` as the last two dimensions. Nr )r r*itemr r r _concatenate_arraysszBaseSampler._concatenate_arrayscCs |dS)a$The start position of the chains. Will raise a ``ValueError`` is ``set_start`` hasn't been run. Returns ------- dict : Dictionary mapping parameters to arrays with shape ``nchains`` giving the starting position of each chain. start_position)r+)r r r r r.s zBaseSampler.start_positioncs2x,t|jD]\}fddD|_q WdS)zSets the starting position of all of the chains. Parameters ---------- positions : dict Dictionary mapping parameter names to arrays of values. The arrays must have shape ``[ntemps x] nchains``. csi|]}|df|qS).r )rr)ii positionsr r sz.BaseSampler.start_position..N) enumeraterr.)r r0chainr )r/r0r r.s  cCs`x0|jD]&}|jt||jt|d7_qWt|gt|j|j}t|t||_dS)zEvolves all of the chains by niterations. Parameters ---------- niterations : int The number of iterations to evolve the chains for. rN)rZ scratchlenmaxr&ziplistr% _evolve_chain)r r)cargsr r r runs &zBaseSampler.runcCsx|jD] }|qWdS)zClears all of the chains.N)rclear)r r3r r r r;s zBaseSampler.clearcCs |jdjS)z.The iteration that the last clear occurred at.r)r lastclear)r r r r r<szBaseSampler.lastclearcCs |dS)z0The history of positions from all of the chains.r0)r-)r r r r r0szBaseSampler.positionscCs |dS)zThe current position of the chains. This will default to the start position if the chains haven't been run yet. current_position)r+)r r r r current_positionsszBaseSampler.current_positionscCs |dS)z,The history of stats from all of the chains.stats)r-)r r r r r?szBaseSampler.statscCs |dS)zThe current stats of the chains. This will default to the stats of the start positions if the chains haven't been run yet. current_stats)r+)r r r r r@ szBaseSampler.current_statscCs |jdjr|d}nd}|S)z7The history of all of the blobs from all of the chains.rblobsN)rhasblobsr-)r rAr r r rAs  zBaseSampler.blobscCs |jdjr|d}nd}|S)zThe current blob data. This will default to the blob data of the start positions if the chains haven't been run yet. rZ current_blobN)rrBr+)r rAr r r current_blobss  zBaseSampler.current_blobscCs |dS)z;The history of all acceptance stats from all of the chains. acceptance)r-)r r r r rD,szBaseSampler.acceptancecCsdd|jDS)ziThe state of all of the chains. Returns a dictionary mapping chain ids to their states. cSsi|]}|j|jqSr )stateZchain_id)rr3r r r r17sz%BaseSampler.state..)r)r r r r rE1szBaseSampler.statecCs*x$t|jD]\}}|||q WdS)zSets the state of the all of the chains. Parameters ---------- state : dict Dictionary mapping chain id to the state to set the chain to. See :py:func:`chain.Chain.set_state` for details. N)r2r set_state)r rEr/r3r r r rF9s zBaseSampler.set_state sampler_statecCst|j|||ddS)aSave the sampler's state to an HDF file. Parameters ---------- fp : :py:class:h5py.File Open file handler to an hdf5 file. The file handler must have write permission. path : str, optional What group to write the state to in the hdf file. Default is the top-level. dsetname : str, optional The name of dataset to store the state to. Default is "sampler_state". )pathdsetnameN)rrE)r fprHrIr r r checkpointEszBaseSampler.checkpointcCs|t||ddS)a%Loads a state from an HDF file. Parameters ---------- fp : :py:class:h5py.File Open file handler to an hdf5 file. path : str, optional What group to store the file to in the hdf file. Default is the top-level ('/'). )rHN)rFr)r rJrHr r r set_state_from_checkpointVs z%BaseSampler.set_state_from_checkpoint)NNN)N)NrG)N)(__name__ __module__ __qualname____doc__r rrrr r"r%propertyr setterrrrrr!r#r'r)rr+r-r.r:r;r<r0r>r?r@rArCrDrErFrKrLr r r r rsT   &               rcCs&|\}}xt|D] }|qW|S)aEvolves a chain for some number of iterations. This is used by ``Sampler.run`` to evolve a collection of chains in a parallel environment. This is not a staticmethod of ``Sampler`` because such functions need to be picklable when using python's multiprocessing pool. Parameters ---------- niterations_chain : tuple of int, chain The number of iterations to run on the given chain. )rangestep)Zniterations_chainr)r3_r r r r7ds  r7) rPabcrrrZepsierrrZepsie.proposalsrrr7r r r r s M