CUQIpy 1.4.0.post0.dev61__py3-none-any.whl → 1.4.0.post0.dev92__py3-none-any.whl

cuqi/__init__.py

@@ -1,6 +1,7 @@
 from . import data
 from . import density
 from . import diagnostics
+from . import algebra
 from . import distribution
 from . import experimental
 from . import geometry

cuqi/_version.py

@@ -8,11 +8,11 @@ import json
 
 version_json = '''
 {
- "date": "2025-11-06T09:26:50+0300",
+ "date": "2025-11-24T16:04:21+0100",
  "dirty": false,
  "error": null,
- "full-revisionid": "e8fd90f6c716b96b04543b4fe62dd637033ccd32",
- "version": "1.4.0.post0.dev61"
+ "full-revisionid": "d6622bf24d4986794de75feebd950e695d1da212",
+ "version": "1.4.0.post0.dev92"
 }
 '''  # END VERSION_JSON
 

cuqi/algebra/__init__.py

@@ -0,0 +1,2 @@
+from ._abstract_syntax_tree import VariableNode, Node
+from ._random_variable import RandomVariable

cuqi/{experimental/algebra/_randomvariable.py → algebra/_random_variable.py}

@@ -1,7 +1,7 @@
 from __future__ import annotations
 from typing import List, Any, Union
-from ._ast import VariableNode, Node
-from ._orderedset import _OrderedSet
+from ._abstract_syntax_tree import VariableNode, Node
+from ._ordered_set import _OrderedSet
 import operator
 import cuqi
 from cuqi.distribution import Distribution
@@ -58,7 +58,7 @@ class RandomVariable:
 
         from cuqi.testproblem import Deconvolution1D
         from cuqi.distribution import Gaussian, Gamma, GMRF
-        from cuqi.experimental.algebra import RandomVariable
+        from cuqi.algebra import RandomVariable
         from cuqi.problem import BayesianProblem
 
         import numpy as np
@@ -79,7 +79,7 @@ class RandomVariable:
     .. code-block:: python
 
         from cuqi.distribution import Gaussian, Gamma
-        from cuqi.experimental.algebra import RandomVariable, VariableNode
+        from cuqi.algebra import RandomVariable, VariableNode
 
         # Define the variables
         x = VariableNode('x')

cuqi/distribution/_distribution.py

@@ -419,5 +419,5 @@ class Distribution(Density, ABC):
     @property
     def rv(self):
         """ Return a random variable object representing the distribution. """
-        from cuqi.experimental.algebra import RandomVariable
+        from cuqi.algebra import RandomVariable
         return RandomVariable(self)

cuqi/distribution/_joint_distribution.py

@@ -62,16 +62,16 @@ class JointDistribution:
         posterior = joint(y=y_obs)
         
     """
-    def __init__(self, *densities: [Density, cuqi.experimental.algebra.RandomVariable]):
+    def __init__(self, *densities: [Density, cuqi.algebra.RandomVariable]):
         """ Create a joint distribution from the given densities. """
 
         # Check if all RandomVariables are simple (not-transformed)
         for density in densities:
-            if isinstance(density, cuqi.experimental.algebra.RandomVariable) and density.is_transformed:
+            if isinstance(density, cuqi.algebra.RandomVariable) and density.is_transformed:
                 raise ValueError(f"To be used in {self.__class__.__name__}, all RandomVariables must be untransformed.")
 
         # Convert potential random variables to their underlying distribution
-        densities = [density.distribution if isinstance(density, cuqi.experimental.algebra.RandomVariable) else density for density in densities]
+        densities = [density.distribution if isinstance(density, cuqi.algebra.RandomVariable) else density for density in densities]
 
         # Ensure all densities have unique names
         names = [density.name for density in densities]

cuqi/experimental/__init__.py

@@ -1,4 +1,2 @@
 """ Experimental module for testing new features and ideas.  """
-from . import algebra
-from . import geometry
 from ._recommender import SamplerRecommender

cuqi/geometry/__init__.py

@@ -16,6 +16,8 @@ from ._geometry import (
     StepExpansion
 )
 
+from ._product_geometry import _ProductGeometry
+
 
 # TODO: We will remove the use of identity geometries in the future
 _identity_geometries = [_DefaultGeometry1D, _DefaultGeometry2D, Continuous1D, Continuous2D, Discrete, Image2D]

cuqi/{experimental/geometry/_productgeometry.py → geometry/_product_geometry.py}

@@ -17,7 +17,7 @@ class _ProductGeometry(Geometry):
     .. code-block:: python
         import numpy as np
         from cuqi.geometry import Continuous1D, Discrete
-        from cuqi.experimental.geometry import _ProductGeometry
+        from cuqi.geometry import _ProductGeometry
         geometry1 = Continuous1D(np.linspace(0, 1, 100))
         geometry2 = Discrete(["sound_speed"])
         product_geometry = _ProductGeometry(geometry1, geometry2)

cuqi/model/_model.py

@@ -229,7 +229,7 @@ class Model(object):
             self._domain_geometry = self._create_default_geometry(value)
         elif isinstance(value, tuple) and self.number_of_inputs > 1:
             geometries = [item if isinstance(item, Geometry) else self._create_default_geometry(item) for item in value]
-            self._domain_geometry = cuqi.experimental.geometry._ProductGeometry(*geometries)
+            self._domain_geometry = cuqi.geometry._ProductGeometry(*geometries)
         elif value is None:
             raise AttributeError(
                 "The parameter 'domain_geometry' is not specified by the user and it cannot be inferred from the attribute 'forward'."
@@ -273,7 +273,7 @@ class Model(object):
         consistent with the forward operator."""
         if (
             not isinstance(
-                self.domain_geometry, cuqi.experimental.geometry._ProductGeometry
+                self.domain_geometry, cuqi.geometry._ProductGeometry
             )
             and self.number_of_inputs > 1
         ):
@@ -448,13 +448,13 @@ class Model(object):
         # If len of kwargs is larger than 1, the geometry needs to be of type
         # _ProductGeometry
         if (
-            not isinstance(geometry, cuqi.experimental.geometry._ProductGeometry)
+            not isinstance(geometry, cuqi.geometry._ProductGeometry)
             and len(kwargs) > 1
         ):
             raise ValueError(
                 "The input is specified by more than one argument. This is only "
                 + "supported for domain geometry of type "
-                + f"{cuqi.experimental.geometry._ProductGeometry.__name__}."
+                + f"{cuqi.geometry._ProductGeometry.__name__}."
             )
 
         # If is_par is bool, make it a tuple of bools of the same length as
@@ -464,7 +464,7 @@ class Model(object):
         # Set up geometries list
         geometries = (
             geometry.geometries
-            if isinstance(geometry, cuqi.experimental.geometry._ProductGeometry)
+            if isinstance(geometry, cuqi.geometry._ProductGeometry)
             else [geometry]
         )
 
@@ -691,7 +691,7 @@ class Model(object):
            is_par and
            len(args) == 1 and
            args[0].shape == (self.domain_dim,) and
-           isinstance(self.domain_geometry, cuqi.experimental.geometry._ProductGeometry)):
+           isinstance(self.domain_geometry, cuqi.geometry._ProductGeometry)):
             # Split the stacked input
             split_args = np.split(args[0], self.domain_geometry.stacked_par_split_indices)
             # Convert split args to CUQIarray if input is CUQIarray
@@ -757,7 +757,7 @@ class Model(object):
             return self._handle_case_when_model_input_is_distributions(kwargs)
 
         # If input is a random variable, we handle it separately
-        elif all(isinstance(x, cuqi.experimental.algebra.RandomVariable)
+        elif all(isinstance(x, cuqi.algebra.RandomVariable)
                for x in kwargs.values()):
             if partial_arguments:
                 raise ValueError(
@@ -769,7 +769,7 @@ class Model(object):
         # We use NotImplemented to indicate that the operation is not supported from the Model class
         # in case of operations such as "@" that can be interpreted as both __matmul__ and __rmatmul__
         # the operation may be delegated to the Node class.
-        elif any(isinstance(args_i, cuqi.experimental.algebra.Node) for args_i in args):
+        elif any(isinstance(args_i, cuqi.algebra.Node) for args_i in args):
             return NotImplemented
 
         # if input is partial, we create a new model with the partial input
@@ -803,7 +803,7 @@ class Model(object):
         if len(distributions) == 1:
             return list(distributions)[0].dim == self.domain_dim
         elif len(distributions) > 1 and isinstance(
-            self.domain_geometry, cuqi.experimental.geometry._ProductGeometry
+            self.domain_geometry, cuqi.geometry._ProductGeometry
         ):
             return all(
                 d.dim == self.domain_geometry.par_dim_list[i]
@@ -851,7 +851,7 @@ class Model(object):
         # Create a partial domain geometry with the geometries corresponding
         # to the non-default arguments that are not in kwargs (remaining
         # unspecified inputs)
-        partial_domain_geometry = cuqi.experimental.geometry._ProductGeometry(
+        partial_domain_geometry = cuqi.geometry._ProductGeometry(
             *[
                 self.domain_geometry.geometries[i]
                 for i in range(self.number_of_inputs)
@@ -864,7 +864,7 @@ class Model(object):
 
         # Create a domain geometry with the geometries corresponding to the
         # non-default arguments that are specified
-        substituted_domain_geometry = cuqi.experimental.geometry._ProductGeometry(
+        substituted_domain_geometry = cuqi.geometry._ProductGeometry(
             *[
                 self.domain_geometry.geometries[i]
                 for i in range(self.number_of_inputs)
@@ -967,7 +967,7 @@ class Model(object):
     def _handle_case_when_model_input_is_random_variables(self, kwargs):
         """ Private function that handles the case of the input being a random variable. """
         # If random variable is not a leaf-type node (e.g. internal node) we return NotImplemented
-        if any(not isinstance(x.tree, cuqi.experimental.algebra.VariableNode) for x in kwargs.values()):
+        if any(not isinstance(x.tree, cuqi.algebra.VariableNode) for x in kwargs.values()):
             return NotImplemented        
 
         # Extract the random variable distributions and check dimensions consistency with domain geometry
@@ -1157,7 +1157,7 @@ class Model(object):
         domain_geometries = (
             self.domain_geometry.geometries
             if isinstance(
-                self.domain_geometry, cuqi.experimental.geometry._ProductGeometry
+                self.domain_geometry, cuqi.geometry._ProductGeometry
             )
             else [self.domain_geometry]
         )
@@ -1185,7 +1185,7 @@ class Model(object):
         # Create list of domain geometries
         geometries = (
             self.domain_geometry.geometries
-            if isinstance(self.domain_geometry, cuqi.experimental.geometry._ProductGeometry)
+            if isinstance(self.domain_geometry, cuqi.geometry._ProductGeometry)
             else [self.domain_geometry]
         )
 
@@ -1197,7 +1197,7 @@ class Model(object):
         # stacked, split it
         if (
             isinstance(
-                self.domain_geometry, cuqi.experimental.geometry._ProductGeometry
+                self.domain_geometry, cuqi.geometry._ProductGeometry
             )
             and not isinstance(grad, (list, tuple))
             and isinstance(grad, np.ndarray)
@@ -1206,7 +1206,7 @@ class Model(object):
 
         # If the domain geometry is not a _ProductGeometry, turn grad into a
         # list of length 1, so that we can iterate over it
-        if not isinstance(self.domain_geometry, cuqi.experimental.geometry._ProductGeometry):
+        if not isinstance(self.domain_geometry, cuqi.geometry._ProductGeometry):
             grad = [grad]
 
         # apply the gradient of each geometry component

cuqi/sampler/_gibbs.py

@@ -180,33 +180,39 @@ class HybridGibbs:
         for sampler in self.samplers.values():
             sampler.validate_target()
 
-    def sample(self, Ns) -> 'HybridGibbs':
+    def sample(self, Ns, Nt=1) -> 'HybridGibbs':
         """ Sample from the joint distribution using Gibbs sampling
 
         Parameters
         ----------
         Ns : int
             The number of samples to draw.
-
+        Nt : int, optional, default=1
+            The thinning interval. If Nt >= 1, every Nt'th sample is stored. The larger Nt, the fewer samples are stored.
+            
         """
         for idx in tqdm(range(Ns), "Sample: "):
 
             self.step()
 
-            self._store_samples()
+            if (Nt > 0) and (idx % Nt == 0):
+                self._store_samples()
 
             # Call callback function if specified
             self._call_callback(idx, Ns)
 
         return self
 
-    def warmup(self, Nb, tune_freq=0.1) -> 'HybridGibbs':
+    def warmup(self, Nb, Nt=1, tune_freq=0.1) -> 'HybridGibbs':
         """ Warmup (tune) the samplers in the Gibbs sampling scheme
 
         Parameters
         ----------
         Nb : int
             The number of samples to draw during warmup.
+        
+        Nt : int, optional, default=1
+            The thinning interval. If Nt >= 1, every Nt'th sample is stored. The larger Nt, the fewer samples are stored.
 
         tune_freq : float, optional
             Frequency of tuning the samplers. Tuning is performed every tune_freq*Nb steps.
@@ -221,9 +227,10 @@ class HybridGibbs:
 
             # Tune the sampler at tuning intervals (matching behavior of Sampler class)
             if (idx + 1) % tune_interval == 0:
-                self.tune(tune_interval, idx // tune_interval) 
-                
-            self._store_samples()
+                self.tune(tune_interval, idx // tune_interval)
+
+            if (Nt > 0) and (idx % Nt == 0):
+                self._store_samples()
 
             # Call callback function if specified
             self._call_callback(idx, Nb)

cuqi/sampler/_pcn.py

@@ -24,7 +24,14 @@ class PCN(Sampler):  # Refactor to Proposal-based sampler?
     def validate_target(self):
         if not isinstance(self.target, cuqi.distribution.Posterior):
             raise ValueError(f"To initialize an object of type {self.__class__}, 'target' need to be of type 'cuqi.distribution.Posterior'.")
-        if not isinstance(self.prior, (cuqi.distribution.Gaussian, cuqi.distribution.Normal)):
+        if not isinstance(
+            self.prior,
+            (
+                cuqi.distribution.Gaussian,
+                cuqi.distribution.Normal,
+                cuqi.distribution.GMRF,
+            ),
+        ):
             raise ValueError("The prior distribution of the target need to be Gaussian")
 
     def step(self):
@@ -46,7 +53,7 @@ class PCN(Sampler):  # Refactor to Proposal-based sampler?
             self.current_point = x_star
             self.current_likelihood_logd = loglike_eval_star
             acc = 1
-        
+
         return acc
 
     @property
@@ -56,7 +63,7 @@ class PCN(Sampler):  # Refactor to Proposal-based sampler?
     @property
     def likelihood(self):
         return self.target.likelihood
-        
+
     def _loglikelihood(self, x):
         return self.likelihood.logd(x)
 

cuqi/sampler/_sampler.py

@@ -236,7 +236,7 @@ class Sampler(ABC):
 
             # Store samples
             self._acc.append(acc)
-            if (Nt > 0) and ((idx + 1) % Nt == 0):
+            if (Nt > 0) and (idx % Nt == 0):
                 self._samples.append(self.current_point)
 
             # display acc rate at progress bar
@@ -285,7 +285,7 @@ class Sampler(ABC):
 
             # Store samples
             self._acc.append(acc)
-            if (Nt > 0) and ((idx + 1) % Nt == 0):
+            if (Nt > 0) and (idx % Nt == 0):
                 self._samples.append(self.current_point)
 
             # display acc rate at progress bar

{cuqipy-1.4.0.post0.dev61.dist-info → cuqipy-1.4.0.post0.dev92.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: CUQIpy
-Version: 1.4.0.post0.dev61
+Version: 1.4.0.post0.dev92
 Summary: Computational Uncertainty Quantification for Inverse problems in Python
 Maintainer-email: "Nicolai A. B. Riis" <nabr@dtu.dk>, "Jakob S. Jørgensen" <jakj@dtu.dk>, "Amal M. Alghamdi" <amaal@dtu.dk>, Chao Zhang <chaz@dtu.dk>
 License:                                  Apache License

{cuqipy-1.4.0.post0.dev61.dist-info → cuqipy-1.4.0.post0.dev92.dist-info}/RECORD

@@ -1,8 +1,12 @@
-cuqi/__init__.py,sha256=TIdTvq0th4ywHVNGFovmEzneXxkpGWC-LlUdrxbPPLM,509
+cuqi/__init__.py,sha256=gs3F4QTPPvZe1sGbPcHn9U4eIyYJHXCJuwK4-eLe0qs,531
 cuqi/_messages.py,sha256=fzEBrZT2kbmfecBBPm7spVu7yHdxGARQB4QzXhJbCJ0,415
-cuqi/_version.py,sha256=S2ddEKfXvlEJwt4HQ_HYcSdVnIse2ZS69CQGg95CvD8,509
+cuqi/_version.py,sha256=THAxuqXPxO5VvVxshGbGP7u9M3U8CCa9aQarrDWeJDw,509
 cuqi/config.py,sha256=wcYvz19wkeKW2EKCGIKJiTpWt5kdaxyt4imyRkvtTRA,526
 cuqi/diagnostics.py,sha256=5OrbJeqpynqRXOe5MtOKKhe7EAVdOEpHIqHnlMW9G_c,3029
+cuqi/algebra/__init__.py,sha256=2uWq2_H1XpXg1TLheTzLL3ySG5tFiNQkZCd9TcZeaEI,99
+cuqi/algebra/_abstract_syntax_tree.py,sha256=PdPz19cJMjvnMx4KEzhn4gvxIZX_UViE33Mbttj_5Xw,9873
+cuqi/algebra/_ordered_set.py,sha256=fKysh4pmI4xF7Y5Z6O86ABzg20o4uBs-v8jmLBMrdpo,2849
+cuqi/algebra/_random_variable.py,sha256=owEkaiF1lD0UUGLS9BOTsLKxL98NRTtmNnyVd9J0Y94,19717
 cuqi/array/__init__.py,sha256=-EeiaiWGNsE3twRS4dD814BIlfxEsNkTCZUc5gjOXb0,30
 cuqi/array/_array.py,sha256=UJoWyQjDpl8g2q4wuuz_ufpr2JoYL0BbCxrxZpZR99I,4303
 cuqi/data/__init__.py,sha256=1aGgPmtG_Kqbb880vLnPksGvyYQB_6o2mz_q-4KGYaU,173
@@ -19,12 +23,12 @@ cuqi/distribution/_beta.py,sha256=QlibnuHNcvWjl-du5aRc9QuzS3n4PsyD_8Nc47w-E0Q,29
 cuqi/distribution/_cauchy.py,sha256=Qwi21WkwUBnBkLbhR-yCGO0tQ_U_3mmvR0pDMPPPB5c,3296
 cuqi/distribution/_cmrf.py,sha256=tCbEulM_O7FB3C_W-3IqZp9zGHkTofCdFF0ybHc9UZI,3745
 cuqi/distribution/_custom.py,sha256=11lnAG7CS15bpV6JOOP--2YyKnAduvXY1IjAOJ1tqO0,10504
-cuqi/distribution/_distribution.py,sha256=FWwJy6kM-MYEfN0Dv7DbLL7JHl_VDNAvNe-GqiLiRBc,18290
+cuqi/distribution/_distribution.py,sha256=ZT8k0J-DNoEptdVBbsTYHBBd17V-6DfCBwo1BatC_0Y,18277
 cuqi/distribution/_gamma.py,sha256=VcvBJS51N-MxuX42r9L2j2QYRlzhdgAtQ6Wa5IFO_YE,3536
 cuqi/distribution/_gaussian.py,sha256=3L1L_3W6i6YuPQ8vnFmju5QsvkLlg4VsgCnj11lYBUE,32977
 cuqi/distribution/_gmrf.py,sha256=OwId8qQWEtmC2fxVhL4iBHZnc8ZCrZzfV6yGXDE3k30,9522
 cuqi/distribution/_inverse_gamma.py,sha256=oPJuiYp3O1m547pmmIz9OWesky9YpwLTHT7-9MmcYss,3159
-cuqi/distribution/_joint_distribution.py,sha256=ALOnQsIrzE8Rx_FYOs4f276u4QZQeN_e0CLC7CJpb-E,20396
+cuqi/distribution/_joint_distribution.py,sha256=WoiaaZ_4nLFJaKUy-Ks5U0ED4DGb3BCKChcuCqJPCgU,20357
 cuqi/distribution/_laplace.py,sha256=5exLvlzJm2AgfvZ3KUSkjfwlGwwbsktBxP8z0iLMik8,1401
 cuqi/distribution/_lmrf.py,sha256=rdGoQ-fPe1oW6Z29P-l3woq0NX3_RxUQ2rzm1VzemNM,3290
 cuqi/distribution/_lognormal.py,sha256=8_hOFQ3iu88ujX8vxmfVEZ0fdmlhTY98PlG5PasPjEg,2612
@@ -34,16 +38,11 @@ cuqi/distribution/_posterior.py,sha256=6LxXAIBzFxyEqx5bUfgqTXVKgXxqNhRItec6FcRVL
 cuqi/distribution/_smoothed_laplace.py,sha256=p-1Y23mYA9omwiHGkEuv3T2mwcPAAoNlCr7T8osNkjE,2925
 cuqi/distribution/_truncated_normal.py,sha256=_ez3MmO6qpBeP6BKCUlW3IgxuF7k--A7jPGPUhtYK0g,4240
 cuqi/distribution/_uniform.py,sha256=fVgj_4SBav8JMc1pNAO1l_CZ9ZwdoMIpN9iQ3i9_Z0Q,3255
-cuqi/experimental/__init__.py,sha256=HCI1aG9xn4XPnBfGbCRrFaq3JfmVyBH09IuvdDancsY,154
+cuqi/experimental/__init__.py,sha256=OZ4WHvYNiQWGdRikzSDhDZcajr9ib35NblQpIOHnQCg,109
 cuqi/experimental/_recommender.py,sha256=ulW69nDGwkB8X93OMgMVLcVz2L0T3SZqHKG7jZ9wNm8,7907
-cuqi/experimental/algebra/__init__.py,sha256=btRAWG58ZfdtK0afXKOg60AX7d76KMBjlZa4AWBCCgU,81
-cuqi/experimental/algebra/_ast.py,sha256=PdPz19cJMjvnMx4KEzhn4gvxIZX_UViE33Mbttj_5Xw,9873
-cuqi/experimental/algebra/_orderedset.py,sha256=fKysh4pmI4xF7Y5Z6O86ABzg20o4uBs-v8jmLBMrdpo,2849
-cuqi/experimental/algebra/_randomvariable.py,sha256=isbFtIWsWXF-yF5Vb56nLy4MCkQM6akjd-dQau4wfbE,19725
-cuqi/experimental/geometry/__init__.py,sha256=kgoKegfz3Jhr7fpORB_l55z9zLZRtloTLyXFDh1oF2o,47
-cuqi/experimental/geometry/_productgeometry.py,sha256=IlBmmKsWE-aRZHp6no9gUXGRfkHlgM0CdPBx1hax9HI,7199
-cuqi/geometry/__init__.py,sha256=Tz1WGzZBY-QGH3c0GiyKm9XHN8MGGcnU6TUHLZkzB3o,842
+cuqi/geometry/__init__.py,sha256=1BcuUHP5tNzfsTLz1qB0WACRIVJU3AtdgfEqgI89oo8,891
 cuqi/geometry/_geometry.py,sha256=W-oQTZPelVS7fN9qZj6bNBuh-yY0eqOHJ39UwB-WmQY,47562
+cuqi/geometry/_product_geometry.py,sha256=tdFUrQzoKPtqbd5aLxXP3DMOkBNrM9B-DTUQPNBxgFE,7186
 cuqi/implicitprior/__init__.py,sha256=6Fl4Lmld8ikg9sW9tReKRGTCJC6_WCTExHaYuIv34nM,323
 cuqi/implicitprior/_regularizedGMRF.py,sha256=BUeT4rwJzary9K56fkxCNGCeKZd-2VSgOT8XNHxFPRE,6345
 cuqi/implicitprior/_regularizedGaussian.py,sha256=9BSKHGEW0OT9OIt_42strDzxBM8mB6A-blcf0kEguHw,21836
@@ -65,7 +64,7 @@ cuqi/legacy/sampler/_sampler.py,sha256=hw5BNhnFkERGtdvZ4Gd_QQZYshNiU_Jx-Yl1--x5H
 cuqi/likelihood/__init__.py,sha256=QXif382iwZ5bT3ZUqmMs_n70JVbbjxbqMrlQYbMn4Zo,1776
 cuqi/likelihood/_likelihood.py,sha256=I12qQF3h_Z8jj7zb_AYD-SEUn34VCU7VxcTcH25Axao,7074
 cuqi/model/__init__.py,sha256=jgY2-jyxEMC79vkyH9BpfowW7_DbMRjqedOtO5fykXQ,62
-cuqi/model/_model.py,sha256=VwNTBlqQ94k_SQB_Y4ADqMOsrrNCcLMii3n-oTAiMv0,75120
+cuqi/model/_model.py,sha256=SWGCzZ2-rmosngulNOofgyMC1K3-WmDH3_Ws3qs8Pmg,74912
 cuqi/operator/__init__.py,sha256=0pc9p-KPyl7KtPV0noB0ddI0CP2iYEHw5rbw49D8Njk,136
 cuqi/operator/_operator.py,sha256=yNwPTh7jR07AiKMbMQQ5_54EgirlKFsbq9JN1EODaQI,8856
 cuqi/pde/__init__.py,sha256=pExaUnEvdNtgL9sQachcKkbDilE87rRiUTKk4hd2zSM,158
@@ -78,14 +77,14 @@ cuqi/sampler/_conjugate.py,sha256=QeQI-ZH4OK_wIIygWLM7dNkH0TSVixdb-VUBFc2prw4,22
 cuqi/sampler/_conjugate_approx.py,sha256=LXPykdG-2KKCCIsxeu4G4Iu6BcAtujtAVOJTGV_aRbc,3525
 cuqi/sampler/_cwmh.py,sha256=0CVKuFuD5-99ipXF-EQ5T5RpYORWriugzttNu6xM7Ik,7320
 cuqi/sampler/_direct.py,sha256=tyVoC9lOGrdR4T0bQUlRewC8uDaenIV61Jk3NEttVVw,767
-cuqi/sampler/_gibbs.py,sha256=H8ZYoNGyQun2u9dU86HDfr5Z1zZO9AnwkeATgS3Jer8,13971
+cuqi/sampler/_gibbs.py,sha256=8V5yBEmBsjiUTVKPTd8nT94AJ4hMraI6LitDDqv7Nrg,14405
 cuqi/sampler/_hmc.py,sha256=yQ0mq51sRMLICRXIdDZLNzlJLt-03OecDOsMkADb-sI,19462
 cuqi/sampler/_langevin_algorithm.py,sha256=2FubHq6p5MruHx8hwLRQXHoHqf5POO4JR6UAkBWV6r4,15029
 cuqi/sampler/_laplace_approximation.py,sha256=tfsfbHhMA_q0pcNF5X8SYT6bWWzY5BSC-hKKRdSaGUI,5906
 cuqi/sampler/_mh.py,sha256=WsVmpvkhvYuFoUvXevOYibIwmhvaKY4YXrkEsF_5YvY,2931
-cuqi/sampler/_pcn.py,sha256=vilU7YjU2WrK5b4bh4VAv8a8TTgj1f9Osivcf3kUlFQ,3364
+cuqi/sampler/_pcn.py,sha256=C8JzklKRXrF-XjEGjyMJ12gZS1ufxQwkTOvhcXlHdZc,3470
 cuqi/sampler/_rto.py,sha256=AQFi8Kqa3uQZsZxdpPb2B9Is_lyfz1v8l9JAY3e4J2w,17186
-cuqi/sampler/_sampler.py,sha256=7_a9i6A7AX3NNz7qK1jTsEYt6bFCUR5WK464KfH_Kvc,21034
+cuqi/sampler/_sampler.py,sha256=rwlaidODyOcG3RQIOaFkPwk0aLZAXdw-3xA35izNgqQ,21022
 cuqi/samples/__init__.py,sha256=vCs6lVk-pi8RBqa6cIN5wyn6u-K9oEf1Na4k1ZMrYv8,44
 cuqi/samples/_samples.py,sha256=hUc8OnCF9CTCuDTrGHwwzv3wp8mG_6vsJAFvuQ-x0uA,35832
 cuqi/solver/__init__.py,sha256=KYgAi_8VoAwljTB3S2I87YnJkRtedskLee7hQp_-zp8,220
@@ -95,8 +94,8 @@ cuqi/testproblem/_testproblem.py,sha256=EJWG_zXUtmo6GlHBZFqHlRpDC_48tE0XZEu0_C66
 cuqi/utilities/__init__.py,sha256=d5QXRzmI6EchS9T4b7eTezSisPWuWklO8ey4YBx9kI0,569
 cuqi/utilities/_get_python_variable_name.py,sha256=wxpCaj9f3ZtBNqlGmmuGiITgBaTsY-r94lUIlK6UAU4,2043
 cuqi/utilities/_utilities.py,sha256=R7BdNysrE36a4D729DvfrTisWY4paP5nfqdkQxSX3Mg,18431
-cuqipy-1.4.0.post0.dev61.dist-info/licenses/LICENSE,sha256=kJWRPrtRoQoZGXyyvu50Uc91X6_0XRaVfT0YZssicys,10799
-cuqipy-1.4.0.post0.dev61.dist-info/METADATA,sha256=8RwtR0AsvJQF8RQuPlj-VpzmMlNRt0bws6uMjfUQf9c,18623
-cuqipy-1.4.0.post0.dev61.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-cuqipy-1.4.0.post0.dev61.dist-info/top_level.txt,sha256=AgmgMc6TKfPPqbjV0kvAoCBN334i_Lwwojc7HE3ZwD0,5
-cuqipy-1.4.0.post0.dev61.dist-info/RECORD,,
+cuqipy-1.4.0.post0.dev92.dist-info/licenses/LICENSE,sha256=kJWRPrtRoQoZGXyyvu50Uc91X6_0XRaVfT0YZssicys,10799
+cuqipy-1.4.0.post0.dev92.dist-info/METADATA,sha256=O6rws-Zn0aPz13jiQkcpRcXcA5HQV2fRG54KOpAwr4Y,18623
+cuqipy-1.4.0.post0.dev92.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+cuqipy-1.4.0.post0.dev92.dist-info/top_level.txt,sha256=AgmgMc6TKfPPqbjV0kvAoCBN334i_Lwwojc7HE3ZwD0,5
+cuqipy-1.4.0.post0.dev92.dist-info/RECORD,,

cuqi/experimental/algebra/__init__.py

@@ -1,2 +0,0 @@
-from ._ast import VariableNode, Node
-from ._randomvariable import RandomVariable

cuqi/experimental/geometry/__init__.py

@@ -1 +0,0 @@
-from ._productgeometry import _ProductGeometry