Release 3.0

.github/workflows/pypi-test.yml

@@ -69,10 +69,6 @@ jobs:
         working-directory: easyscience
         run: pixi run python -m pytest ../tests/unit/ --color=yes -v
 
-      - name: Run functional tests to verify the installation
-        working-directory: easyscience
-        run: pixi run python -m pytest ../tests/functional/ --color=yes -v
-
       - name: Run integration tests to verify the installation
         working-directory: easyscience
         run: pixi run python -m pytest ../tests/integration/ --color=yes -n auto

.github/workflows/test.yml

@@ -102,21 +102,6 @@ jobs:
             pixi run --environment $env unit-tests
           done
 
-      - name: Run functional tests
-        shell: bash
-        run: |
-          set -euo pipefail
-
-          for py_ver in $PY_VERSIONS; do
-            echo
-            echo "🔹🔸🔹🔸🔹 Python: $py_ver 🔹🔸🔹🔸🔹"
-
-            env="py-$(echo $py_ver | tr -d .)-env" # Converts 3.XX -> py-3XX-env
-
-            echo "Running tests in environment: $env"
-            pixi run --environment $env functional-tests
-          done
-
       - name: Run integration tests ${{ needs.env-prepare.outputs.pytest-marks }}
         shell: bash
         run: |
@@ -258,25 +243,6 @@ jobs:
             cd ..
           done
 
-      - name: Run functional tests
-        shell: bash
-        run: |
-          set -euo pipefail
-
-          for py_ver in $PY_VERSIONS; do
-            echo
-            echo "🔹🔸🔹🔸🔹 Python: $py_ver 🔹🔸🔹🔸🔹"
-
-            echo "Entering pixi project directory easyscience_py$py_ver"
-            cd easyscience_py$py_ver
-
-            echo "Running tests"
-            pixi run python -m pytest ../tests/functional/ --color=yes -v
-
-            echo "Exiting pixi project directory"
-            cd ..
-          done
-
       - name: Run integration tests ${{ needs.env-prepare.outputs.pytest-marks }}
         shell: bash
         run: |

docs/docs/api-reference/index.md

@@ -8,19 +8,16 @@ This section contains the reference detailing the functions and modules
 available in EasyScience.
 
 - [base_classes](base_classes.md) – Core abstract and helper base
-  classes used to build EasyScience objects (e.g. `ObjBase`,
-  `ModelBase`).
-- [fitting](fitting.md) – Fitting utilities and interfaces, including
-  `Fitter` and available minimizers.
+  classes used to build EasyScience objects (e.g. `EasyBase`,
+  `ModelBase`, `EasyList`, `SerializerBase`).
+- [fitting](fitting.md) – Fitting and sampling, including `Fitter`,
+  `Sampler` and the available minimizers.
 - [global_object](global_object.md) – Global singleton providing shared
-  services (logger, map, undo/redo stack, script manager).
-- [io](io.md) – Serialization and deserialization framework for
-  persisting EasyScience objects (serializers and components).
+  services (logger and object registry).
 - [job](job.md) – Job and experiment abstractions for running and
   organizing analyses.
-- [models](models.md) – Predefined model implementations (e.g.
-  polynomial models) used in fitting workflows.
 - [utils](utils.md) – Miscellaneous utility functions and helpers (class
   tools, decorators, type helpers).
 - [variable](variable.md) – Descriptor types and variable abstractions
-  (e.g. `DescriptorNumber`, `Parameter`, `DescriptorArray`).
+  (e.g. `DescriptorNumber`, `Parameter`, `DescriptorArray`). All of them
+  are `EasyBase` objects, serialized with `to_dict`/`from_dict`.

docs/docs/tutorials/fitting-bayesian.ipynb

@@ -28,7 +28,7 @@
     "\n",
     "where $\\theta$ are the model parameters, $d$ is the observed data, $p(d \\mid \\theta)$ is the likelihood, and $p(\\theta)$ is the prior. In `easyscience`, the `min`/`max` bounds of a `Parameter` are interpreted as a **uniform prior**, and a Gaussian likelihood is constructed from the data and supplied weights.\n",
     "\n",
-    "`easyscience` exposes a Bayesian Markov-chain Monte Carlo (MCMC) sampler through the `Sampler` class.  Under the hood this uses BUMPS' DREAM sampler, so the underlying minimizer must be switched to BUMPS.\n",
+    "`easyscience` exposes a Bayesian Markov-chain Monte Carlo (MCMC) sampler through the `Sampler` class.\n",
     "\n",
     "```{note}\n",
     "This tutorial focuses on Bayesian analysis with a simple QENS model for illustration. For dedicated QENS fitting with more sophisticated models, consider using [`EasyDynamics`](https://github.com/easyscience/easydynamics).\n",
@@ -138,12 +138,12 @@
   },
   {
    "cell_type": "markdown",
-   "id": "9",
+   "id": "7",
    "metadata": {},
    "source": [
     "## Defining parameters with priors\n",
     "\n",
-    "Create four `Parameter` objects, for the area $A$, $\\gamma$, $\\omega_0$ and $\\sigma$. The `min` and `max` arguments define a **uniform prior** on each parameter — the sampler will only consider values inside this range and will treat every value inside the range as equally plausible *a priori*.\n",
+    "Create four `Parameter` objects, for the area $A$, $\\gamma$, $\\omega_0$ and $\\sigma$. The `min` and `max` arguments define a **uniform prior** on each parameter — the sampler will only consider values inside this range and will treat every value inside the range as equally plausible *a priori*. The parameters are then collected in an `EasyList` container: this is the model object that both `Fitter` and `Sampler` take, so it is defined here, before the optional fit.\n",
     "\n",
     "| Parameter | Initial Value | Min | Max |\n",
     "| --- | --- | --- | --- |\n",
@@ -156,21 +156,24 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "10",
+   "id": "8",
    "metadata": {},
    "outputs": [],
    "source": [
+    "from easyscience import EasyList\n",
     "from easyscience import Parameter\n",
     "\n",
-    "area = Parameter(name='area', value=10, fixed=False, min=1, max=100)\n",
-    "gamma = Parameter(name='gamma', value=8e-3, fixed=False, min=1e-4, max=1e-2)\n",
-    "omega_0 = Parameter(name='omega_0', value=1e-3, fixed=False, min=0, max=2e-3)\n",
-    "sigma = Parameter(name='sigma', value=1e-3, fixed=False, min=1e-5, max=1e-1)"
+    "area = Parameter(display_name='area', value=10, fixed=False, min=1, max=100)\n",
+    "gamma = Parameter(display_name='gamma', value=8e-3, fixed=False, min=1e-4, max=1e-2)\n",
+    "omega_0 = Parameter(display_name='omega_0', value=1e-3, fixed=False, min=0, max=2e-3)\n",
+    "sigma = Parameter(display_name='sigma', value=1e-3, fixed=False, min=1e-5, max=1e-1)\n",
+    "\n",
+    "parameter_container = EasyList(area, gamma, omega_0, sigma)"
    ]
   },
   {
    "cell_type": "markdown",
-   "id": "77ce3f63",
+   "id": "9",
    "metadata": {},
    "source": [
     "## The model\n",
@@ -193,7 +196,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "36a0c4f4",
+   "id": "10",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -217,9 +220,7 @@
    "source": [
     "## Maximum-likelihood fit (optional, but recommended)\n",
     "\n",
-    "Perform a quick maximum-likelihood fit. This is **not** a prerequisite for sampling: `Sampler`\n",
-    "needs a *configured* `Fitter`, not a *fitted* one, and you can sample straight from the initial\n",
-    "parameter values.\n",
+    "Perform a quick maximum-likelihood fit.\n",
     "\n",
     "It is worth doing anyway, for two reasons:\n",
     "\n",
@@ -239,9 +240,6 @@
    "outputs": [],
    "source": [
     "from easyscience import Fitter\n",
-    "from easyscience import ObjBase\n",
-    "\n",
-    "parameter_container = ObjBase(name='params', A=area, gamma=gamma, omega_0=omega_0, sigma=sigma)\n",
     "\n",
     "mle_fitter = Fitter(parameter_container, intensity_model)\n",
     "mle_result = mle_fitter.fit(x=omega, y=intensity_obs, weights=1 / intensity_error)\n",
@@ -261,7 +259,7 @@
     "\n",
     "We now draw samples from the posterior distribution $p(\\theta \\mid d)$ using the BUMPS DREAM (DiffeRential Evolution Adaptive Metropolis) algorithm. DREAM is an ensemble MCMC method that runs multiple chains in parallel and automatically tunes the proposal distribution.\n",
     "\n",
-    "DREAM only works with the BUMPS minimizer. We reuse the ``mle_fitter`` created above — any configured `Fitter` would do, and it does not have to have been fitted — switch it to BUMPS, and create a `Sampler` instance bound to the fitter and data. Calling `sampler.sample()` returns a `SamplingResults` object with the following attributes:\n",
+    "Create a `Sampler` from the same `parameter_container` and `intensity_model` we gave the `Fitter`, bound to the data. Calling `sampler.sample()` returns a `SamplingResults` object with the following attributes:\n",
     "\n",
     "- `draws`: a `(n_samples, n_parameters)` array of posterior samples: each **row** is one complete draw from the joint posterior (one value for every parameter simultaneously), and each **column** holds all sampled values for a single parameter. Note this is a *trimmed* view of the chain rather than the raw buffer, so `n_samples` is smaller than `samples / thin` — see the note under [Extend the chain](#extend-the-chain-and-check-convergence);\n",
     "- `param_names`: the unique names of the parameters, in the same column order as `draws`;\n",
@@ -274,7 +272,9 @@
     "- `burn` (500): the number of initial *burn-in* generations to discard — the sampler needs time to find the typical set of the posterior, and early samples are not representative. Note this counts generations, not raw samples, so `burn=500` discards `500 × n_chains` raw samples;\n",
     "- `thin` (2): the *thinning* interval — only every second generation is kept, which reduces autocorrelation between consecutive draws;\n",
     "\n",
-    "First, we switch to the BUMPS minimizer:"
+    "```{note}\n",
+    "If you already have a `Fitter`, `Sampler.from_fitter(mle_fitter, omega, intensity_obs, weights=1 / intensity_error)` builds the same sampler from the model bound to it. The fitter is neither needed afterwards nor modified.\n",
+    "```"
    ]
   },
   {
@@ -283,22 +283,12 @@
    "id": "14",
    "metadata": {},
    "outputs": [],
-   "source": [
-    "from easyscience import AvailableMinimizers\n",
-    "\n",
-    "mle_fitter.switch_minimizer(AvailableMinimizers.Bumps)"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "5a219fcd",
-   "metadata": {},
-   "outputs": [],
    "source": [
     "from easyscience.fitting import Sampler\n",
     "\n",
-    "sampler = Sampler(mle_fitter, omega, intensity_obs, weights=1 / intensity_error)\n",
+    "sampler = Sampler(\n",
+    "    parameter_container, intensity_model, omega, intensity_obs, weights=1 / intensity_error\n",
+    ")\n",
     "results = sampler.sample(samples=10000, burn=500, thin=2)\n",
     "\n",
     "print(f'Drew {results.draws.shape[0]} samples for {results.draws.shape[1]} parameters.')\n",
@@ -307,7 +297,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "8766b170",
+   "id": "15",
    "metadata": {},
    "source": [
     "## Convergence diagnostics\n",
@@ -323,7 +313,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "3c49ab6f",
+   "id": "16",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -359,7 +349,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "15",
+   "id": "17",
    "metadata": {},
    "source": [
     "## Posterior summaries\n",
@@ -377,15 +367,15 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "ce3e38a8",
+   "id": "18",
    "metadata": {},
    "outputs": [],
    "source": []
   },
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "16",
+   "id": "19",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -407,7 +397,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "17",
+   "id": "20",
    "metadata": {},
    "source": [
     "## Visualise the joint posterior\n",
@@ -418,7 +408,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "18",
+   "id": "21",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -452,7 +442,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "19",
+   "id": "22",
    "metadata": {},
    "source": [
     "## Posterior-predictive band\n",
@@ -463,7 +453,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "20",
+   "id": "23",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -499,7 +489,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "03339658",
+   "id": "24",
    "metadata": {},
    "source": [
     "## Extend the chain and check convergence\n",
@@ -543,7 +533,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "293b140b",
+   "id": "25",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -573,7 +563,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "9ec4302c",
+   "id": "26",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -614,7 +604,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "b0f30be6",
+   "id": "27",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -641,7 +631,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "50b7213a",
+   "id": "28",
    "metadata": {},
    "source": [
     "### What is Gelman-Rubin R-hat?\n",
@@ -675,7 +665,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "3449e0a7",
+   "id": "29",
    "metadata": {},
    "outputs": [],
    "source": [

docs/docs/tutorials/fitting-qens.ipynb

@@ -246,10 +246,10 @@
    "source": [
     "from easyscience import Parameter\n",
     "\n",
-    "A = Parameter(name='A', value=10, fixed=False, min=0.01, max=100)\n",
-    "gamma = Parameter(name='gamma', value=8e-3, fixed=False, min=1e-4, max=10e-3)\n",
-    "omega_0 = Parameter(name='omega_0', value=0.001, fixed=False, min=0, max=0.002)\n",
-    "sigma = Parameter(name='sigma', value=0.001, fixed=False, min=1e-5, max=0.1)"
+    "A = Parameter(display_name='A', value=10, fixed=False, min=0.01, max=100)\n",
+    "gamma = Parameter(display_name='gamma', value=8e-3, fixed=False, min=1e-4, max=10e-3)\n",
+    "omega_0 = Parameter(display_name='omega_0', value=0.001, fixed=False, min=0, max=0.002)\n",
+    "sigma = Parameter(display_name='sigma', value=0.001, fixed=False, min=1e-5, max=0.1)"
    ]
   },
   {
@@ -293,10 +293,10 @@
    "metadata": {},
    "outputs": [],
    "source": [
+    "from easyscience import EasyList\n",
     "from easyscience import Fitter\n",
-    "from easyscience import ObjBase\n",
     "\n",
-    "params = ObjBase(name='params', A=A, gamma=gamma, omega_0=omega_0, sigma=sigma)\n",
+    "params = EasyList(A, gamma, omega_0, sigma)\n",
     "f = Fitter(params, intensity)\n",
     "\n",
     "res = f.fit(x=q, y=i, weights=1 / di)"