mc3 3.2.1__tar.gz → 3.3.0__tar.gz

mc3-3.3.0/.github/workflows/python-package.yml

@@ -0,0 +1,47 @@
+# This workflow will install Python dependencies, run tests and lint with a variety of Python versions
+# For more information see: https://help.github.com/actions/language-and-framework-guides/using-python-with-github-actions
+
+name: Tests
+
+on:
+  push:
+    branches: [ "*" ]
+  pull_request:
+    branches: [ "main" ]
+  schedule:
+    - cron: '0 12 1 * *'
+
+jobs:
+  build:
+
+    runs-on: ${{ matrix.os }}
+    strategy:
+      fail-fast: false
+      matrix:
+        python-version: ["3.10", "3.11", "3.12", "3.13", "3.14"]
+        os: [ubuntu-latest]
+
+    steps:
+    - uses: actions/checkout@v6
+    - name: Set up Python ${{ matrix.python-version }}
+      uses: actions/setup-python@v6
+      with:
+        python-version: ${{ matrix.python-version }}
+    - name: Install dependencies
+      run: |
+        python -m pip install --upgrade pip
+        python -m pip install flake8 pytest
+        if [ -f requirements.txt ]; then pip install -r requirements.txt; fi
+
+    - name: Install package
+      run: |
+        python -m pip install -e .
+    - name: Lint with flake8
+      run: |
+        # stop the build if there are Python syntax errors or undefined names
+        flake8 . --count --select=E9,F63,F7,F82 --show-source --statistics
+        # exit-zero treats all errors as warnings. The GitHub editor is 127 chars wide
+        flake8 . --count --exit-zero --max-complexity=10 --max-line-length=127 --statistics
+    - name: Test with pytest
+      run: |
+        pytest tests -v

mc3-3.3.0/MANIFEST.in

@@ -0,0 +1,12 @@
+include MANIFEST.in
+include CONTRIBUTING.md
+include README.md
+include LICENSE
+include pyproject.toml
+include requirements.txt
+include pytest.ini
+include .readthedocs.yml
+include .github/workflows/python-package.yml
+include tests/*
+include docs/*
+include examples/*.py

{mc3-3.2.1 → mc3-3.3.0}/PKG-INFO

@@ -1,37 +1,19 @@
-Metadata-Version: 2.2
+Metadata-Version: 2.4
 Name: mc3
-Version: 3.2.1
+Version: 3.3.0
 Summary: Multi-core Markov-chain Monte Carlo package
 Author-email: Patricio Cubillos <patricio.cubillos@oeaw.ac.at>
-License: Copyright (c) 2015-2022 Patricio Cubillos and contributors.
-        
-        Permission is hereby granted, free of charge, to any person obtaining a copy
-        of this software and associated documentation files (the "Software"), to deal
-        in the Software without restriction, including without limitation the rights
-        to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
-        copies of the Software, and to permit persons to whom the Software is
-        furnished to do so, subject to the following conditions:
-        
-        The above copyright notice and this permission notice shall be included in
-        all copies or substantial portions of the Software.
-        
-        THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
-        IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
-        FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.  IN NO EVENT SHALL THE
-        AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
-        LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
-        OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
-        THE SOFTWARE.
-        
+License-Expression: MIT
 Project-URL: Homepage, https://github.com/pcubillos/mc3
-Requires-Python: >=3.9
+Requires-Python: >=3.10
 Description-Content-Type: text/markdown
 License-File: LICENSE
-Requires-Dist: numpy>=1.23.5
+Requires-Dist: numpy>=2.0.0
 Requires-Dist: scipy>=1.5.4
 Requires-Dist: matplotlib>=3.3.4
 Provides-Extra: test
 Requires-Dist: pytest>=6.0; extra == "test"
+Dynamic: license-file
 
 # mc3: Multi-core Markov-chain Monte Carlo
 > A Python implementation of the Markov-chain Monte Carlo algorithm.

mc3-3.3.0/docs/docs_requirements.txt

@@ -0,0 +1,6 @@
+sphinx==8.1.3
+sphinx-rtd-theme>=2.0.0
+nbsphinx
+sphinx_copybutton
+ipython>=8.17
+

{mc3-3.2.1 → mc3-3.3.0}/docs/get_started.rst

@@ -6,7 +6,7 @@ Getting Started
 System Requirements
 -------------------
 
-``mc3`` is compatible with Python 3.9+, and has been `tested
+``mc3`` is compatible with Python 3.10+, and is `tested
 <https://github.com/pcubillos/mc3/actions/workflows/python-package.yml>`_ to work on Unix/Linux and OS X machines.
 
 Install

{mc3-3.2.1 → mc3-3.3.0}/docs/index.rst

@@ -40,7 +40,7 @@ mc3: Multi-Core Markov-Chain Monte Carlo
 
 .. note::
 
-    Got  Python3.9+? Simply install as: ``pip install mc3``
+    Compatible with Python3.10.  Install as: ``pip install mc3``
 
 
 Features

{mc3-3.2.1 → mc3-3.3.0}/mc3/__init__.py

@@ -1,4 +1,4 @@
-# Copyright (c) 2015-2025 Patricio Cubillos and contributors.
+# Copyright (c) 2015-2026 Patricio Cubillos and contributors.
 # mc3 is open-source software under the MIT license (see LICENSE).
 
 from .sampler_driver import *

{mc3-3.2.1 → mc3-3.3.0}/mc3/__main__.py

@@ -1,4 +1,4 @@
-# Copyright (c) 2015-2025 Patricio Cubillos and contributors.
+# Copyright (c) 2015-2026 Patricio Cubillos and contributors.
 # mc3 is open-source software under the MIT license (see LICENSE).
 
 import sys

{mc3-3.2.1 → mc3-3.3.0}/mc3/chain.py

@@ -1,4 +1,4 @@
-# Copyright (c) 2015-2025 Patricio Cubillos and contributors.
+# Copyright (c) 2015-2026 Patricio Cubillos and contributors.
 # mc3 is open-source software under the MIT license (see LICENSE).
 
 import sys

{mc3-3.2.1 → mc3-3.3.0}/mc3/fit_driver.py

@@ -1,4 +1,4 @@
-# Copyright (c) 2015-2025 Patricio Cubillos and contributors.
+# Copyright (c) 2015-2026 Patricio Cubillos and contributors.
 # mc3 is open-source software under the MIT license (see LICENSE).
 
 __all__ = [

{mc3-3.2.1 → mc3-3.3.0}/mc3/mcmc_driver.py

@@ -1,4 +1,4 @@
-# Copyright (c) 2015-2025 Patricio Cubillos and contributors.
+# Copyright (c) 2015-2026 Patricio Cubillos and contributors.
 # mc3 is open-source software under the MIT license (see LICENSE).
 
 __all__ = [

{mc3-3.2.1 → mc3-3.3.0}/mc3/plots/__init__.py

@@ -1,4 +1,4 @@
-# Copyright (c) 2015-2025 Patricio Cubillos and contributors.
+# Copyright (c) 2015-2026 Patricio Cubillos and contributors.
 # mc3 is open-source software under the MIT license (see LICENSE).
 
 from .plot_functions import *

{mc3-3.2.1 → mc3-3.3.0}/mc3/plots/colors.py

@@ -1,4 +1,4 @@
-# Copyright (c) 2015-2025 Patricio Cubillos and contributors.
+# Copyright (c) 2015-2026 Patricio Cubillos and contributors.
 # mc3 is open-source software under the MIT license (see LICENSE).
 
 __all__ = [

{mc3-3.2.1 → mc3-3.3.0}/mc3/plots/plot_functions.py

@@ -1,4 +1,4 @@
-# Copyright (c) 2015-2025 Patricio Cubillos and contributors.
+# Copyright (c) 2015-2026 Patricio Cubillos and contributors.
 # mc3 is open-source software under the MIT license (see LICENSE).
 
 __all__ = [

{mc3-3.2.1 → mc3-3.3.0}/mc3/plots/posterior.py

@@ -1,4 +1,4 @@
-# Copyright (c) 2015-2025 Patricio Cubillos and contributors.
+# Copyright (c) 2015-2026 Patricio Cubillos and contributors.
 # mc3 is open-source software under the MIT license (see LICENSE).
 
 __all__ = [
@@ -140,6 +140,8 @@ def _histogram(
         hist_kw['facecolor'] = 'none'
 
     for i in range(npars):
+        if np.any(np.isnan(posterior[:,i])):
+            continue
         ax = axes[i]
         if clear:
             ax.clear()
@@ -533,6 +535,7 @@ class Marginal():
             bins=25, fontsize=11, linewidth=1.5,
             axes=None,
             show_texts=True, show_estimates=True,
+            alpha=0.6,
         ):
         self.source = source
         self.fig = None
@@ -544,6 +547,7 @@ class Marginal():
         self.bestp = bestp
         self.ranges = ranges
         self.theme = theme
+        self.alpha = alpha
 
         if nx is None or ny is None:
             # Default layout:
@@ -630,6 +634,7 @@ class Marginal():
             hpd_min, self.source.low_bounds, self.source.high_bounds,
             self.linewidth, self.theme,
             self.orientation,
+            alpha=self.alpha,
             top_pad=1.25,
         )
         _plot_marginal(self)
@@ -637,7 +642,7 @@ class Marginal():
         if savefile is not None:
             self.fig.savefig(savefile, dpi=300)
 
-    def overplot(self, post, labels=None, nlevels=4, alpha=0.4):
+    def overplot(self, posts, labels=None, nlevels=4, alpha=0.4, axes_mask=None):
         """
         Overplot additional posteriors in the same figure.
 
@@ -659,25 +664,39 @@ class Marginal():
             length of labels has to be one more than posts, because
             it also contains the label for self.
         """
-        if '_like' in post.statistics:
-            hpd_min = post.hpd_min
-        else:
+        all_posts = [self.source] + list(posts)
+        for i,post in enumerate(posts):
+            estimates = post.estimates
+            if not self.show_estimates:
+                estimates = [None for _ in estimates]
             hpd_min = None
+            if '_like' in self.statistics:
+                hpd_min = self.source.hpd_min
 
-        estimates = post.estimates
-        if not self.show_estimates:
-            estimates = [None for _ in estimates]
+            _histogram(
+                post.posterior, estimates, self.ranges,
+                self.hist_axes, self.bins,
+                post.pdf, post.xpdf,
+                hpd_min, post.low_bounds, post.high_bounds,
+                self.linewidth, post.theme,
+                post.orientation, alpha=alpha,
+                clear=False,
+            )
+        if not self.show_texts:
+            return
+        text_cols = [post.theme.color for post in all_posts]
+        for text in self.stats_texts:
+            text.set_visible(False)
+        self.stats_texts = []
 
-        _histogram(
-            post.posterior, estimates, post.ranges,
-            self.hist_axes, self.bins,
-            post.pdf, post.xpdf,
-            hpd_min, post.low_bounds, post.high_bounds,
-            self.linewidth, post.theme,
-            post.orientation,
-            top_pad=1.25,
-            clear=False,
-        )
+        for j in range(self.npars):
+            stats_texts = [
+                '' if post.mask[j] else rf'{post.tex_estimates[j]}'
+                for post in all_posts
+            ]
+            self.stats_texts += colors.rainbow_text(
+                self.hist_axes[j], stats_texts, self.fontsize, text_cols,
+            )
 
 
 class Figure(Marginal):
@@ -694,10 +713,11 @@ class Figure(Marginal):
             bins=25, nlevels=6, fontsize=None, linewidth=None,
             show_texts=True, show_estimates=True,
             show_colorbar=True,
+            figure=None,
             fignum=None,
         ):
         self.source = source
-        self.fig = None
+        self.fig = figure
         self.pair_axes = None
         self.hist_axes = None
         self.posterior = posterior
@@ -934,7 +954,7 @@ class Figure(Marginal):
                 post.posterior, estimates, self.ranges, self.hist_axes,
                 self.bins, post.pdf, post.xpdf,
                 hpd_min, post.low_bounds, post.high_bounds,
-                self.linewidth, post.theme, post.orientation, alpha=0.5,
+                self.linewidth, post.theme, post.orientation, alpha=alpha,
                 clear=False,
             )
         if not self.show_texts:
@@ -1044,7 +1064,7 @@ class StatisticsUpdate(ShareUpdate):
             )
 
 
-class Posterior(object):
+class Posterior():
     """
     Classification of posterior plotting tools.
 
@@ -1090,6 +1110,7 @@ class Posterior(object):
             show_texts=True, show_estimates=True,
             show_colorbar=True,
             seed=314159,
+            mask=None,
         ):
         self.figures = []
         nsamples, self.npars = np.shape(posterior)
@@ -1131,7 +1152,9 @@ class Posterior(object):
             self.bestp = bestp
         self.statistics = statistics
         self.quantile = quantile
-
+        if mask is None:
+            mask = np.ones(self.npars, bool)
+        self.mask = mask
 
     def plot(
             self, plot_marginal=True,
@@ -1143,6 +1166,7 @@ class Posterior(object):
             show_texts=None, show_estimates=None,
             show_colorbar=None,
             savefile=None,
+            nlevels=None,
         ):
         """
         Plot marginal histograms and pairwise posteriors.
@@ -1156,6 +1180,8 @@ class Posterior(object):
             show_texts = self.show_texts
         if show_colorbar is None:
             show_colorbar = self.show_colorbar
+        if nlevels is None:
+            nlevels = 6
 
         fig = Figure(
             self,
@@ -1174,7 +1200,8 @@ class Posterior(object):
             show_texts=show_texts,
             show_estimates=show_estimates,
             show_colorbar=show_colorbar,
-            # bins=25, nlevels=20,
+            # bins=25,
+            nlevels=nlevels,
         )
         self.figures.append(fig)
         fig.plot(figure=figure, savefile=savefile)
@@ -1185,7 +1212,9 @@ class Posterior(object):
             self, fignum=None, axes=None, quantile=None,
             nx=None, ny=None,
             savefile=None,
+            fontsize=11, linewidth=1.5,
             show_texts=None, show_estimates=None,
+            alpha=0.6,
         ):
         """
         Plot histogram of marginal posteriors.
@@ -1194,6 +1223,8 @@ class Posterior(object):
             show_estimates = self.show_estimates
         if show_texts is None:
             show_texts = self.show_texts
+        if quantile is None:
+            quantile = self.quantile
 
         fig = Marginal(
             self,
@@ -1201,8 +1232,12 @@ class Posterior(object):
             self.ranges, self.theme,
             statistics=self.statistics,
             nx=nx, ny=ny,
+            quantile=quantile,
+            fontsize=fontsize,
+            linewidth=linewidth,
             show_texts=show_texts,
             show_estimates=show_estimates,
+            alpha=alpha,
         )
         self.figures.append(fig)
         fig.plot(savefile=savefile, fignum=fignum, axes=axes)

{mc3-3.2.1 → mc3-3.3.0}/mc3/sampler_driver.py

@@ -1,4 +1,4 @@
-# Copyright (c) 2015-2025 Patricio Cubillos and contributors.
+# Copyright (c) 2015-2026 Patricio Cubillos and contributors.
 # mc3 is open-source software under the MIT license (see LICENSE).
 
 __all__ = [

{mc3-3.2.1 → mc3-3.3.0}/mc3/stats/__init__.py

@@ -1,15 +1,17 @@
-# Copyright (c) 2015-2025 Patricio Cubillos and contributors.
+# Copyright (c) 2015-2026 Patricio Cubillos and contributors.
 # mc3 is open-source software under the MIT license (see LICENSE).
 
 from .gelman import *
 from .stats import *
 from .time_averaging import *
+from .dwt import *
 from .prayer import *
 
 __all__ = (
     gelman.__all__
   + stats.__all__
   + time_averaging.__all__
+  + dwt.__all__
   + prayer.__all__
     )
 

mc3-3.3.0/mc3/stats/dwt.py

@@ -0,0 +1,229 @@
+# Copyright (c) 2015-2026 Patricio Cubillos and contributors.
+# mc3 is open-source software under the MIT license (see LICENSE).
+
+__all__ = [
+    'daubechies4',
+    'dwt_daub4',
+    'dwt_chisq',
+]
+
+import numpy as np
+
+
+def daubechies4(array, n, inverse):
+    """
+    Applies the Daubechies-4 discrete wavelet transform (DWT)
+    or its inverse (transpose filter) depending on isign.
+
+    Parameters
+    ----------
+    a : np.ndarray
+        Input signal (1D array of length n).
+        Modified in place.
+    n: int
+        Hierarchy level of the transform.
+    inverse: bool
+        Calculate forward (False) or inverse (True) DWT
+    """
+    C0 =  0.48296291314453410
+    C1 =  0.83651630373780772
+    C2 =  0.22414386804201339
+    C3 = -0.12940952255126034
+
+    if len(array) < 4:
+        return
+
+    dwt = np.zeros(n)
+    k = n >> 1  # n/2 using bit shift
+
+    if inverse:
+        dwt[0] = C2*array[k-1] + C1*array[n-1] + C0*array[0] + C3*array[k]
+        dwt[1] = C3*array[k-1] - C0*array[n-1] + C1*array[0] - C2*array[k]
+        j = 2
+        for i in range(0, k-1):
+            dwt[j]   = C2*array[i] + C1*array[i+k] + C0*array[i+1] + C3*array[i+k+1]
+            dwt[j+1] = C3*array[i] - C0*array[i+k] + C1*array[i+1] - C2*array[i+k+1]
+            j += 2
+    else:  # Forward transform
+        i = 0
+        for j in range(0, n-3, 2):
+            dwt[i]   = C0*array[j] + C1*array[j+1] + C2*array[j+2] + C3*array[j+3]
+            dwt[i+k] = C3*array[j] - C2*array[j+1] + C1*array[j+2] - C0*array[j+3]
+            i += 1
+        # Wrap-around (circular convolution)
+        dwt[i]   = C0*array[n-2] + C1*array[n-1] + C2*array[0] + C3*array[1]
+        dwt[i+k] = C3*array[n-2] - C2*array[n-1] + C1*array[0] - C0*array[1]
+
+    # Copy back into input array
+    array[0:n] = dwt
+
+
+def dwt_daub4(array, inverse=False):
+    """
+    1D discrete wavelet transform using the Daubechies 4-parameter wavelet
+    Implements the pyramid algorithm replacing array by its DWT or inverse DWT
+
+    Parameters
+    ----------
+    array: 1D float array
+        Input data array to which to apply the DWT (modified in place)
+    inverse: bool
+        Calculate forward (False) or inverse (True) DWT
+
+    Notes
+    -----
+    The input vector must have length 2**M with M an integer, otherwise
+    the output will zero-padded to the next size of the form 2**M.
+
+    Examples
+    --------
+    >>> import numpy as np
+    >>> import matplotlib.pyplot as plt
+    >>> import mc3.stats as ms
+
+    >>> # Calculate the inverse DWT for a unit vector:
+    >>> nx = 1000
+    >>> array = np.zeros(nx)
+    >>> array[10:20] = 1.0
+    >>> inverse = True
+    >>> idwt = ms.dwt_daub4(array, inverse=inverse)
+    >>> dwt_array = dawt(array, inverse=inverse)
+
+    >>> # Plot the inverse DWT
+    >>> fig = plt.figure(10)
+    >>> fig.set_size_inches(7,4)
+    >>> plt.clf()
+    >>> plt.subplots_adjust(0.1, 0.1, 0.99, 0.96)
+    >>> plt.plot(1-dwt_array/idwt, label='py DWT', color='k')
+    >>> plt.plot(0.1*array, label='signal', color='tomato')
+    >>> plt.plot(idwt, label='DWT', color='xkcd:blue')
+    >>> plt.plot(dwt_array, label='py DWT', color='limegreen', dashes=(8,3))
+    """
+    size = len(array)
+    M = np.ceil(np.log2(size))
+    dwt_size = 2**int(M)
+
+    dwt_array = np.zeros(dwt_size)
+    dwt_array[0:size] = array
+
+    if size < 4:
+        return dwt_array
+
+    if inverse:
+        # Inverse transform: start from smallest hierarchy to largest
+        n = 4
+        while n <= dwt_size:
+            daubechies4(dwt_array, n, inverse)
+            n <<= 1
+    else:
+        # Forward transform: start from largest hierarchy to smallest
+        n = dwt_size
+        while n >= 4:
+            daubechies4(dwt_array, n, inverse)
+            n >>= 1
+
+    return dwt_array
+
+
+def dwt_chisq(model, data, params, priors=None, priorlow=None, priorup=None):
+    """
+    Calculate -2*ln(likelihood) in a wavelet-base (a pseudo chi-squared)
+    based on Carter & Winn (2009), ApJ 704, 51.
+
+    Parameters
+    ----------
+    model: 1D ndarray
+        Model fit of data.
+    data: 1D ndarray
+        Data set array fitted by model.
+    params: 1D float ndarray
+        Model parameters (including the tree noise parameters: gamma,
+        sigma_r, sigma_w; which must be the last three elements in params).
+    priors: 1D ndarray
+        Parameter prior values.
+    priorlow: 1D ndarray
+        Left-sided prior standard deviation (param < prior).
+        A priorlow value of zero denotes a uniform prior.
+    priorup: 1D ndarray
+        Right-sided prior standard deviation (prior < param).
+        A priorup value of zero denotes a uniform prior.
+
+    Returns
+    -------
+    chisq: Float
+        Wavelet-based (pseudo) chi-squared.
+
+    Notes
+    -----
+    - If the residuals array size is not of the form 2**N, the routine
+    zero-pads the array until this condition is satisfied.
+
+    Examples
+    --------
+    >>> import mc3.stats as ms
+    >>> import numpy as np
+    >>> # Compute chi-squared for a given model fitting a data set:
+    >>> data = np.array([2.0, 0.0, 3.0, -2.0, -1.0, 2.0, 2.0, 0.0])
+    >>> model = np.ones(len(data))
+    >>> params = np.array([1.0, 0.1, 0.1])
+    >>> chisq = ms.dwt_chisq(model, data, params)
+    >>> print(chisq)
+    1693.22308882
+
+    >>> # Now, say this is a three-parameter model, with a Gaussian prior
+    >>> # on the last parameter:
+    >>> priors = np.array([1.0, 0.2, 0.3])
+    >>> plow   = np.array([0.0, 0.0, 0.1])
+    >>> pup    = np.array([0.0, 0.0, 0.1])
+    >>> chisq = ms.dwt_chisq(model, data, params, priors, plow, pup)
+    >>> print(chisq)
+    1697.2230888243134
+    """
+    if len(params) < 3:
+        raise ValueError('Wavelet chisq should have at least three parameters')
+
+    # g-factor of the covariance of the wavelet coefficients:
+    # g(gamma=1) = 1.0 / (2*ln(2))
+    g = 0.72134752
+
+    size = len(data)
+    M = int(np.ceil(np.log2(size)))
+
+    # Expand (zero padding) to res to a size proportional to 2**M (zero padding)
+    residuals = np.zeros(2**M)
+    residuals[0:size] = data - model
+    dwt_residuals = dwt_daub4(residuals, inverse=False)
+
+    gamma, sigma_r, sigma_w = params[-3:]
+    # Equation (34) of CW2009, square of sigma_S
+    sigma_S = g * sigma_r**2.0 * 2.0**-gamma + sigma_w**2.0
+    # Second term in right-hand side of equation 32
+    chisq = (
+        (dwt_residuals[0]**2.0 + dwt_residuals[1]**2.0) / sigma_S +
+        2.0*np.log(2*np.pi*sigma_S)
+    )
+
+    # Equation (33) of CW2009, sigma_W squared
+    for m in range(1, M): # Number of scales
+        sigma_W = sigma_r**2.0 * 2.0**(-gamma*m) + sigma_w**2.0
+        n = 2**m
+        res_sq = np.sum(dwt_residuals[n:2*n]**2.0)
+        chisq += res_sq/sigma_W + n*np.log(2*np.pi*sigma_W)
+
+    # Add priors contribution
+    if priors is None or priorlow is None or priorup is None:
+        return chisq
+
+    delta_priors = params - priors
+    for d_prior, e_low, e_up in zip(delta_priors, priorlow, priorup):
+        # No prior
+        if e_low<=0 or e_up<=0:
+            continue
+        # Gaussian prior
+        if d_prior > 0:
+            chisq += (d_prior/e_up)**2.0
+        else:
+            chisq += (d_prior/e_low)**2.0
+
+    return chisq
+

{mc3-3.2.1 → mc3-3.3.0}/mc3/stats/gelman.py

@@ -1,4 +1,4 @@
-# Copyright (c) 2015-2025 Patricio Cubillos and contributors.
+# Copyright (c) 2015-2026 Patricio Cubillos and contributors.
 # mc3 is open-source software under the MIT license (see LICENSE).
 
 __all__ = [

{mc3-3.2.1 → mc3-3.3.0}/mc3/stats/prayer.py

@@ -1,4 +1,4 @@
-# Copyright (c) 2015-2025 Patricio Cubillos and contributors.
+# Copyright (c) 2015-2026 Patricio Cubillos and contributors.
 # mc3 is open-source software under the MIT license (see LICENSE).
 
 __all__ = [