o h=(@sdZddlZddlmZmZgdZGdddZeZeZe ej e ej ddZ e ej d d Z e ej d d Ze eje eje eje ejd dZe eje eje eje ejddZe eje ejddZe eje eje eje ejddZe ejddZe ejddZe ejddZe eje ej ddZ!e ej"e ej"ddZ#e ej$dd Z%e ej$d!d"Z&e ej'd#d$Z(e ej'd%d&Z)dS)'aF PyTorch provides two global :class:`ConstraintRegistry` objects that link :class:`~torch.distributions.constraints.Constraint` objects to :class:`~torch.distributions.transforms.Transform` objects. These objects both input constraints and return transforms, but they have different guarantees on bijectivity. 1. ``biject_to(constraint)`` looks up a bijective :class:`~torch.distributions.transforms.Transform` from ``constraints.real`` to the given ``constraint``. The returned transform is guaranteed to have ``.bijective = True`` and should implement ``.log_abs_det_jacobian()``. 2. ``transform_to(constraint)`` looks up a not-necessarily bijective :class:`~torch.distributions.transforms.Transform` from ``constraints.real`` to the given ``constraint``. The returned transform is not guaranteed to implement ``.log_abs_det_jacobian()``. The ``transform_to()`` registry is useful for performing unconstrained optimization on constrained parameters of probability distributions, which are indicated by each distribution's ``.arg_constraints`` dict. These transforms often overparameterize a space in order to avoid rotation; they are thus more suitable for coordinate-wise optimization algorithms like Adam:: loc = torch.zeros(100, requires_grad=True) unconstrained = torch.zeros(100, requires_grad=True) scale = transform_to(Normal.arg_constraints['scale'])(unconstrained) loss = -Normal(loc, scale).log_prob(data).sum() The ``biject_to()`` registry is useful for Hamiltonian Monte Carlo, where samples from a probability distribution with constrained ``.support`` are propagated in an unconstrained space, and algorithms are typically rotation invariant.:: dist = Exponential(rate) unconstrained = torch.zeros(100, requires_grad=True) sample = biject_to(dist.support)(unconstrained) potential_energy = -dist.log_prob(sample).sum() .. note:: An example where ``transform_to`` and ``biject_to`` differ is ``constraints.simplex``: ``transform_to(constraints.simplex)`` returns a :class:`~torch.distributions.transforms.SoftmaxTransform` that simply exponentiates and normalizes its inputs; this is a cheap and mostly coordinate-wise operation appropriate for algorithms like SVI. In contrast, ``biject_to(constraints.simplex)`` returns a :class:`~torch.distributions.transforms.StickBreakingTransform` that bijects its input down to a one-fewer-dimensional space; this a more expensive less numerically stable transform but is needed for algorithms like HMC. The ``biject_to`` and ``transform_to`` objects can be extended by user-defined constraints and transforms using their ``.register()`` method either as a function on singleton constraints:: transform_to.register(my_constraint, my_transform) or as a decorator on parameterized constraints:: @transform_to.register(MyConstraintClass) def my_factory(constraint): assert isinstance(constraint, MyConstraintClass) return MyTransform(constraint.param1, constraint.param2) You can create your own registry by creating a new :class:`ConstraintRegistry` object. N) constraints transforms)ConstraintRegistry biject_to transform_tocs2eZdZdZfddZd ddZddZZS) rz5 Registry to link constraints to transforms. csi|_tdSN) _registrysuper__init__)self __class__]/var/www/html/ai/venv/lib/python3.10/site-packages/torch/distributions/constraint_registry.pyr TszConstraintRegistry.__init__Ncs\|dur fddSttjrtttr ttjs'td|j<|S)a Registers a :class:`~torch.distributions.constraints.Constraint` subclass in this registry. Usage:: @my_registry.register(MyConstraintClass) def construct_transform(constraint): assert isinstance(constraint, MyConstraint) return MyTransform(constraint.arg_constraints) Args: constraint (subclass of :class:`~torch.distributions.constraints.Constraint`): A subclass of :class:`~torch.distributions.constraints.Constraint`, or a singleton object of the desired class. factory (Callable): A callable that inputs a constraint object and returns a :class:`~torch.distributions.transforms.Transform` object. Ncs |Sr)register)factory constraintr rrks z-ConstraintRegistry.register..zLExpected constraint to be either a Constraint subclass or instance, but got ) isinstancer Constrainttype issubclass TypeErrorrr rrrrrrXs   zConstraintRegistry.registercCs@z |jt|}W||Stytdt|jddw)aq Looks up a transform to constrained space, given a constraint object. Usage:: constraint = Normal.arg_constraints['scale'] scale = transform_to(constraint)(torch.zeros(1)) # constrained u = transform_to(constraint).inv(scale) # unconstrained Args: constraint (:class:`~torch.distributions.constraints.Constraint`): A constraint object. Returns: A :class:`~torch.distributions.transforms.Transform` object. Raises: `NotImplementedError` if no transform has been registered. zCannot transform z constraintsN)rrKeyErrorNotImplementedError__name__rrrr__call__{s zConstraintRegistry.__call__r)r __module__ __qualname____doc__r rr __classcell__rrr rrOs   #rcCstjSr)ridentity_transformrrrr_transform_to_realsr%cCt|j}t||jSr)rbase_constraintrIndependentTransformreinterpreted_batch_ndimsrbase_transformrrr_biject_to_independent r,cCr&r)rr'rr(r)r*rrr_transform_to_independentr-r.cCtSr)r ExpTransformr$rrr_transform_to_positivesr1cCttt|jdgS)N)rComposeTransformr0AffineTransform lower_boundr$rrr_transform_to_greater_thans  r7cCr2)N)rr4r0r5 upper_boundr$rrr_transform_to_less_thans  r:cCslt|jtjo |jdk}t|jtjo|jdk}|r |r tS|j}|j|j}ttt||gS)Nrr3) rr6numbersNumberr9rSigmoidTransformr4r5)r lower_is_0 upper_is_1locscalerrr_transform_to_intervals rBcCr/r)rStickBreakingTransformr$rrr_biject_to_simplexrDcCr/r)rSoftmaxTransformr$rrr_transform_to_simplexrErGcCr/r)rLowerCholeskyTransformr$rrr_transform_to_lower_choleskyrErIcCr/r)rPositiveDefiniteTransformr$rrr_transform_to_positive_definiterKcCr/r)rCorrCholeskyTransformr$rrr_transform_to_corr_choleskyrLrNcCtdd|jD|j|jS)NcSg|]}t|qSrr.0crrr z"_biject_to_cat..r CatTransformcseqdimlengthsr$rrr_biject_to_cat r\cCrO)NcSrPrrrRrrrrUrVz%_transform_to_cat..rWr$rrr_transform_to_catr]r_cCtdd|jD|jS)NcSrPrrQrRrrrrUrVz$_biject_to_stack..rStackTransformrYrZr$rrr_biject_to_stackrccCr`)NcSrPrr^rRrrrrU#rVz'_transform_to_stack..rar$rrr_transform_to_stack rdre)*r!r;torch.distributionsrr__all__rrrrrealr% independentr,r.positive nonnegativer1 greater_thangreater_than_eqr7 less_thanr:intervalhalf_open_intervalrBsimplexrDrGlower_choleskyrIpositive_definitepositive_semidefiniterK corr_choleskyrNcatr\r_stackrcrerrrrshCI