PySDM 3.0.0rc41__tar.gz → 3.0.0rc42__tar.gz

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/.github/dependabot.yml

@@ -4,10 +4,21 @@ updates:
   directory: "/"
   schedule:
     interval: daily
+  ignore:
+  - dependency-name: "setuptools-scm"
+    update-types:
+      - "version-update:semver-patch" 
+  - dependency-name: "setuptools"
+    update-types:
+      - "version-update:semver-patch" 
 - package-ecosystem: pip
   directory: "/examples"
   schedule:
     interval: daily
+  ignore:
+  - dependency-name: "*"
+    update-types:
+      - "version-update:semver-patch" 
 - package-ecosystem: "github-actions"
   directory: "/"
   schedule:

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/.github/workflows/docs_code_snippets.yml

@@ -91,11 +91,11 @@ jobs:
         with:
           release: R2024b
           cache: true
-      - uses: matlab-actions/run-command@v3.2
+      - uses: matlab-actions/run-command@v3.3
         with:
           command: docs_code
       - run: sed -i -e 's/CPU/GPU/g' docs_code.m
-      - uses: matlab-actions/run-command@v3.2
+      - uses: matlab-actions/run-command@v3.3
         with:
           command: docs_code
 

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/.github/workflows/precommit.yml

@@ -24,7 +24,7 @@ jobs:
       - uses: actions/checkout@v4.1.6
       - uses: actions/setup-python@v7.0.0
         with:
-          python-version: "3.11"
+          python-version: "3.13"
       - run: |
           pip install pre-commit
           pre-commit run --all-files --show-diff-on-failure

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/.pre-commit-config.yaml

@@ -23,7 +23,7 @@ repos:
         files: '.json'
 
   - repo: https://github.com/lk16/detect-missing-init
-    rev: v0.1.6
+    rev: v0.1.9
     hooks:
     - id: detect-missing-init
       args: ['--create', "--track", "--python-folders", "tests,PySDM,examples/PySDM_examples"]

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pysdm
-Version: 3.0.0rc41
+Version: 3.0.0rc42
 Summary: Pythonic particle-based (super-droplet) warm-rain/aqueous-chemistry cloud microphysics package with box, parcel & 1D/2D prescribed-flow examples in Python, Julia and Matlab
 Author-email: "https://github.com/open-atmos/PySDM/graphs/contributors" <sylwester.arabas@agh.edu.pl>
 License-Expression: GPL-3.0-only
@@ -39,7 +39,7 @@ Requires-Dist: pytest; extra == "nonunit-tests"
 Requires-Dist: PySDM-examples; extra == "nonunit-tests"
 Requires-Dist: PyPartMC; extra == "nonunit-tests"
 Provides-Extra: ci-version-pins
-Requires-Dist: PyPartMC==2.1.0; extra == "ci-version-pins"
+Requires-Dist: PyPartMC==2.1.2; extra == "ci-version-pins"
 Requires-Dist: numba==0.60.0; extra == "ci-version-pins"
 Requires-Dist: numpy==2.0.2; extra == "ci-version-pins"
 Requires-Dist: scipy==1.15.3; extra == "ci-version-pins"
@@ -47,9 +47,9 @@ Requires-Dist: ThrustRTC==0.3.20; extra == "ci-version-pins"
 Requires-Dist: CURandRTC==0.1.7; extra == "ci-version-pins"
 Requires-Dist: jax==0.8.0; extra == "ci-version-pins"
 Requires-Dist: Pint==0.24.4; extra == "ci-version-pins"
-Requires-Dist: chempy==0.10.1; extra == "ci-version-pins"
+Requires-Dist: chempy==0.10.2; extra == "ci-version-pins"
 Requires-Dist: pyevtk==1.7.0; extra == "ci-version-pins"
-Requires-Dist: pyparsing==3.2.5; extra == "ci-version-pins"
+Requires-Dist: pyparsing==3.3.3; extra == "ci-version-pins"
 Dynamic: license-file
 
 # <img src="https://raw.githubusercontent.com/open-atmos/PySDM/main/docs/logos/pysdm_logo.svg" width=100 height=146 alt="pysdm logo">

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PySDM/attributes/impl/temperature_variation_option_attribute.py

@@ -7,7 +7,7 @@ class TemperatureVariationOptionAttribute:  # pylint: disable=too-few-public-met
 
     def __init__(self, particulator, neglect_temperature_variations: bool):
         if neglect_temperature_variations:
-            assert particulator.environment.mesh.dimension == 0
+            assert particulator.environment.mesh.n_dims == 0
         self.neglect_temperature_variations = neglect_temperature_variations
         self.initial_temperature = (
             particulator.Storage.from_ndarray(

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PySDM/attributes/numerics/cell_origin.py

@@ -12,5 +12,5 @@ class CellOrigin(CellAttribute):
             particulator,
             name="cell origin",
             dtype=int,
-            n_vector_components=particulator.mesh.dim,
+            n_vector_components=particulator.mesh.n_dims,
         )

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PySDM/attributes/numerics/position_in_cell.py

@@ -11,5 +11,5 @@ class PositionInCell(CellAttribute):
         super().__init__(
             particulator,
             name="position in cell",
-            n_vector_components=particulator.mesh.dim,
+            n_vector_components=particulator.mesh.n_dims,
         )

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PySDM/attributes/physics/critical_saturation.py

@@ -12,7 +12,7 @@ from PySDM.attributes.impl import (
 @register_attribute()
 class CriticalSaturation(DerivedAttribute, TemperatureVariationOptionAttribute):
     def __init__(self, particulator, neglect_temperature_variations=False):
-        assert particulator.mesh.dimension == 0
+        assert particulator.mesh.n_dims == 0
 
         self.v_crit = particulator.get_attribute("critical volume")
         self.v_dry = particulator.get_attribute("dry volume")

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PySDM/backends/__init__.py

@@ -84,10 +84,18 @@ else:
 _BACKEND_CACHE = {}
 
 
-def _cached_backend(formulae=None, backend_class=None, **kwargs):
-    key = backend_class.__name__ + ":" + str(formulae) + ":" + str(kwargs)
+def _cached_backend(formulae=None, n_dims=0, backend_class=None, **kwargs):
+    key = (
+        backend_class.__name__
+        + ":"
+        + str(n_dims)
+        + ":"
+        + str(formulae)
+        + ":"
+        + str(kwargs)
+    )
     if key not in _BACKEND_CACHE:
-        _BACKEND_CACHE[key] = backend_class(formulae=formulae, **kwargs)
+        _BACKEND_CACHE[key] = backend_class(formulae=formulae, n_dims=n_dims, **kwargs)
     return _BACKEND_CACHE[key]
 
 

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PySDM/backends/impl_common/backend_methods.py

@@ -17,3 +17,5 @@ class BackendMethods:
             self.default_jit_flags = {}
         if not hasattr(self, "block_until_ready"):
             self.block_until_ready = False
+        if not hasattr(self, "n_dims"):
+            self.n_dims = None

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PySDM/backends/impl_numba/methods/displacement_methods.py

@@ -107,9 +107,8 @@ class DisplacementMethods(BackendMethods):
     def calculate_displacement(
         self, *, dim, displacement, courant, cell_origin, position_in_cell, n_substeps
     ):
-        n_dims = len(courant.shape)
         scheme = self.formulae.particle_advection.displacement
-        if n_dims == 1:
+        if self.n_dims == 1:
             DisplacementMethods.calculate_displacement_body_1d(
                 dim,
                 scheme,
@@ -119,7 +118,7 @@ class DisplacementMethods(BackendMethods):
                 position_in_cell.data,
                 n_substeps,
             )
-        elif n_dims == 2:
+        elif self.n_dims == 2:
             DisplacementMethods.calculate_displacement_body_2d(
                 dim,
                 scheme,
@@ -129,7 +128,7 @@ class DisplacementMethods(BackendMethods):
                 position_in_cell.data,
                 n_substeps,
             )
-        elif n_dims == 3:
+        elif self.n_dims == 3:
             DisplacementMethods.calculate_displacement_body_3d(
                 dim,
                 scheme,

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PySDM/backends/jax.py

@@ -26,15 +26,17 @@ class Jax(
         self,
         formulae=None,
         *,
+        n_dims=0,
         double_precision=True,
         override_jit_flags=None,  # pylint: disable=unused-argument
         # TODO #1913: investigate if there are any jit/jax flags we can add configuration for
         block_until_ready=False,
     ):
         jax.config.update("jax_enable_x64", True)
-        if not double_precision:
+        if not double_precision or n_dims > 0:
             raise NotImplementedError()
 
+        self.n_dims = n_dims
         self.block_until_ready = (
             block_until_ready  # TODO #1913: implement switch in jit code
         )

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PySDM/backends/numba.py

@@ -40,13 +40,18 @@ class Numba(  # pylint: disable=too-many-ancestors,duplicate-code
     default_croupier = "local"
 
     def __init__(
-        self, formulae=None, *, double_precision=True, override_jit_flags=None
+        self,
+        formulae=None,
+        *,
+        n_dims=0,
+        double_precision=True,
+        override_jit_flags=None,
     ):
         if not double_precision:
             raise NotImplementedError()
         self.formulae = formulae or Formulae()
         self.formulae_flattened = self.formulae.flatten
-
+        self.n_dims = n_dims
         parallel_default = True
 
         if override_jit_flags is not None and "parallel" in override_jit_flags:

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PySDM/backends/thrust_rtc.py

@@ -47,10 +47,16 @@ class ThrustRTC(  # pylint: disable=duplicate-code,too-many-ancestors
     default_croupier = "global"
 
     def __init__(
-        self, formulae=None, *, double_precision=False, debug=False, verbose=False
+        self,
+        formulae=None,
+        *,
+        n_dims=0,
+        double_precision=False,
+        debug=False,
+        verbose=False,
     ):
         self.formulae = formulae or Formulae()
-
+        self.n_dims = n_dims
         self._conv_function = trtc.DVDouble if double_precision else trtc.DVFloat
         self._real_type = "double" if double_precision else "float"
         self._np_dtype = np.float64 if double_precision else np.float32

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PySDM/environments/box.py

@@ -16,6 +16,8 @@ class Box:
         self.particulator = None
         self._ambient_air = {}
         self.backend = backend
+        if backend is not None:
+            assert backend.n_dims == self.mesh.n_dims
 
     def __getitem__(self, item):
         return self._ambient_air[item]

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PySDM/environments/impl/moist.py

@@ -18,6 +18,7 @@ class Moist:  # pylint: disable=too-many-instance-attributes
         self.dt = dt
         self.mesh = mesh
         self.variables = variables
+        assert backend.n_dims == self.mesh.n_dims
 
         if backend.formulae.ventilation.__name__ != "Neglect":
             for var in ("air density", "air dynamic viscosity"):

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PySDM/environments/impl/register_environment.py

@@ -5,8 +5,8 @@ from copy import deepcopy
 
 
 def _instantiate(self, particulator):
-    """Creating a copy without backend as a workaround
-    for long execution times: see PR #1885"""  # to be addressed in TODO #1179
+    """Creating a deep copy to facilitate environment reuse in user code,
+    but exluding backend (and the formulae within) to avoid re-JIT-compilations"""
     backend = self.backend
     self.backend = None
     copy = deepcopy(self)

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PySDM/environments/parcel.py

@@ -39,7 +39,7 @@ class Parcel(Moist):  # pylint: disable=too-many-instance-attributes
         super().__init__(
             dt, Mesh.mesh_0d(), variables, mixed_phase=mixed_phase, backend=backend
         )
-
+        assert backend.n_dims == self.mesh.n_dims
         self.mass_of_dry_air = mass_of_dry_air
         self.w = w if callable(w) else lambda _: w
         self.delta_liquid_water_mixing_ratio = np.nan

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PySDM/exporters/vtk_exporter.py

@@ -89,7 +89,7 @@ class VTKExporter:
             }
         )
 
-        if particulator.mesh.dimension == 2:
+        if particulator.mesh.n_dims == 2:
             y = (
                 particulator.mesh.size[0]
                 / particulator.mesh.grid[0]
@@ -122,7 +122,7 @@ class VTKExporter:
                 print("Exporting Products to vtk, path: " + path)
             payload = {}
 
-            if particulator.mesh.dimension != 2:
+            if particulator.mesh.n_dims != 2:
                 raise NotImplementedError(
                     "Only 2 dimensions data is supported at the moment."
                 )

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PySDM/impl/mesh.py

@@ -21,14 +21,6 @@ class Mesh:
     def dz(self):
         return self.size[-1] / self.grid[-1]
 
-    @property
-    def dimension(self):
-        return self.n_dims
-
-    @property
-    def dim(self):
-        return self.n_dims
-
     @property
     def domain_bottom_surface_area(self):
         assert self.n_dims > 0

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PySDM/initialisation/sampling/spatial_sampling.py

@@ -2,9 +2,6 @@
 spatial sampling logic (i.e., physical x-y-z coordinates)
 """
 
-# TODO #305 QUASIRANDOM & GRID
-#  https://slayoo.github.io/workshop_2019/files/talk_Shima.pdf
-
 
 class Pseudorandom:  # pylint: disable=too-few-public-methods
     @staticmethod

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/PySDM/particulator.py

@@ -725,7 +725,7 @@ class Particulator:  # pylint: disable=too-many-public-methods,too-many-instance
                 **self.condensation_params,
             )
         attributes["multiplicity"] = int_caster(attributes["multiplicity"])
-        if self.mesh.dimension == 0:
+        if self.mesh.n_dims == 0:
             attributes["cell id"] = np.zeros_like(
                 attributes["multiplicity"], dtype=np.int64
             )

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/examples/PySDM_examples/Arabas_et_al_2015/example_benchmark.py

@@ -54,7 +54,7 @@ def main():
         for sd in n_sd:
             settings.n_sd_per_gridbox = sd
             storage = Storage()
-            simulation = Simulation(settings, storage, None, backend())
+            simulation = Simulation(settings, storage, None, backend(n_dims=2))
             simulation.reinit(products=[WallTime()])
             simulation.run()
             times[key].append(storage.load("wall time")[-1])

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/examples/PySDM_examples/Shipway_and_Hill_2012/simulation.py

@@ -56,7 +56,7 @@ class Simulation:
             thd_of_z=settings.thd,
             rhod_of_z=settings.rhod,
             z0=-settings.particle_reservoir_depth,
-            backend=backend(formulae=settings.formulae),
+            backend=backend(formulae=settings.formulae, n_dims=1),
             solvers=mpdata,
         )
 

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/examples/PySDM_examples/Ware_et_al_2025/fig_6_7_8.ipynb

@@ -145,7 +145,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "backend = CPU()\n",
+    "backend = CPU(n_dims=2)\n",
     "runs = {}"
    ]
   },

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/examples/PySDM_examples/utils/kinematic_2d/simulation.py

@@ -25,7 +25,7 @@ class Simulation:
         self.settings = settings
         self.storage = storage
         self.particulator = None
-        self.backend = backend or CPU(settings.formulae)
+        self.backend = backend or CPU(settings.formulae, n_dims=2)
         self.SpinUp = SpinUp
 
     @property

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/examples/pyproject.toml

@@ -48,7 +48,7 @@ CI_version_pins = [
     "ipywidgets==8.1",
     "seaborn==0.13.2",
     "numdifftools==0.11.1",
-    "vtk==9.6.2",
+    "vtk==9.7",
     "pyrcel==2.0.0",
     "pyvinecopulib==1.0.0",
     "cmcrameri==1.10",

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/pyproject.toml

@@ -64,7 +64,7 @@ nonunit-tests = [
 ]
 
 CI_version_pins = [
-  "PyPartMC==2.1.0",
+  "PyPartMC==2.1.2",
   "numba==0.60.0",
   "numpy==2.0.2",
   "scipy==1.15.3",
@@ -72,9 +72,9 @@ CI_version_pins = [
   "CURandRTC==0.1.7",
   "jax==0.8.0",
   "Pint==0.24.4",
-  "chempy==0.10.1",
+  "chempy==0.10.2",
   "pyevtk==1.7.0",
-  "pyparsing==3.2.5",
+  "pyparsing==3.3.3",
 ]
 
 [project.urls]

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/pysdm.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pysdm
-Version: 3.0.0rc41
+Version: 3.0.0rc42
 Summary: Pythonic particle-based (super-droplet) warm-rain/aqueous-chemistry cloud microphysics package with box, parcel & 1D/2D prescribed-flow examples in Python, Julia and Matlab
 Author-email: "https://github.com/open-atmos/PySDM/graphs/contributors" <sylwester.arabas@agh.edu.pl>
 License-Expression: GPL-3.0-only
@@ -39,7 +39,7 @@ Requires-Dist: pytest; extra == "nonunit-tests"
 Requires-Dist: PySDM-examples; extra == "nonunit-tests"
 Requires-Dist: PyPartMC; extra == "nonunit-tests"
 Provides-Extra: ci-version-pins
-Requires-Dist: PyPartMC==2.1.0; extra == "ci-version-pins"
+Requires-Dist: PyPartMC==2.1.2; extra == "ci-version-pins"
 Requires-Dist: numba==0.60.0; extra == "ci-version-pins"
 Requires-Dist: numpy==2.0.2; extra == "ci-version-pins"
 Requires-Dist: scipy==1.15.3; extra == "ci-version-pins"
@@ -47,9 +47,9 @@ Requires-Dist: ThrustRTC==0.3.20; extra == "ci-version-pins"
 Requires-Dist: CURandRTC==0.1.7; extra == "ci-version-pins"
 Requires-Dist: jax==0.8.0; extra == "ci-version-pins"
 Requires-Dist: Pint==0.24.4; extra == "ci-version-pins"
-Requires-Dist: chempy==0.10.1; extra == "ci-version-pins"
+Requires-Dist: chempy==0.10.2; extra == "ci-version-pins"
 Requires-Dist: pyevtk==1.7.0; extra == "ci-version-pins"
-Requires-Dist: pyparsing==3.2.5; extra == "ci-version-pins"
+Requires-Dist: pyparsing==3.3.3; extra == "ci-version-pins"
 Dynamic: license-file
 
 # <img src="https://raw.githubusercontent.com/open-atmos/PySDM/main/docs/logos/pysdm_logo.svg" width=100 height=146 alt="pysdm logo">

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/pysdm.egg-info/requires.txt

@@ -10,7 +10,7 @@ pyevtk
 pyparsing
 
 [CI_version_pins]
-PyPartMC==2.1.0
+PyPartMC==2.1.2
 numba==0.60.0
 numpy==2.0.2
 scipy==1.15.3
@@ -18,9 +18,9 @@ ThrustRTC==0.3.20
 CURandRTC==0.1.7
 jax==0.8.0
 Pint==0.24.4
-chempy==0.10.1
+chempy==0.10.2
 pyevtk==1.7.0
-pyparsing==3.2.5
+pyparsing==3.3.3
 
 [nonunit-tests]
 pytest

pysdm-3.0.0rc42/pysdm.egg-info/scm_version.json

@@ -0,0 +1,8 @@
+{
+  "tag": "3.0.0rc42",
+  "distance": 0,
+  "node": "gd1360ae719a1a243d2b8c55885023654601892dd",
+  "dirty": false,
+  "branch": "HEAD",
+  "node_date": "2026-10-08"
+}

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/tests/smoke_tests/kinematic_2d/arabas_et_al_2015/test_export.py

@@ -16,14 +16,16 @@ from PySDM import Formulae
 from PySDM.exporters import NetCDFExporter, VTKExporter
 
 
-def test_export(backend_instance, tmp_path):
+def test_export(backend_class, tmp_path):
     # Arrange
     settings = Settings()
     settings.simulation_time = settings.dt
     settings.output_interval = settings.dt
 
     storage = Storage()
-    simulator = Simulation(settings, storage, SpinUp=SpinUp, backend=backend_instance)
+    simulator = Simulation(
+        settings, storage, SpinUp=SpinUp, backend=backend_class(n_dims=2)
+    )
     _, temp_file = tempfile.mkstemp(dir=tmp_path, suffix=".nc")
     sut = NetCDFExporter(storage, settings, simulator, temp_file)
 

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/tests/smoke_tests/kinematic_2d/arabas_et_al_2015/test_initialisation.py

@@ -8,13 +8,15 @@ from PySDM.physics.constants import si
 
 
 # pylint: disable=too-many-locals
-def test_initialisation(backend_instance, plot=False):
+def test_initialisation(backend_class, plot=False):
     settings = Settings()
     settings.simulation_time = -1 * settings.dt
     settings.grid = (10, 5)
     settings.n_sd_per_gridbox = 5000
 
-    simulation = Simulation(settings, None, SpinUp=SpinUp, backend=backend_instance)
+    simulation = Simulation(
+        settings, None, SpinUp=SpinUp, backend=backend_class(n_dims=2)
+    )
 
     n_levels = settings.grid[1]
     n_cell = int(np.prod(np.array(settings.grid)))

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/tests/smoke_tests/kinematic_2d/arabas_et_al_2015/test_spin_up.py

@@ -18,7 +18,7 @@ from .dummy_storage import DummyStorage
         pytest.param(True, id="fastmath: True"),
     ),
 )
-def test_spin_up(backend_instance, fastmath, plot=False):
+def test_spin_up(backend_class, fastmath, plot=False):
     # Arrange
     settings = Settings(Formulae(fastmath=fastmath))
     settings.dt = 0.5 * si.second
@@ -27,7 +27,9 @@ def test_spin_up(backend_instance, fastmath, plot=False):
     settings.output_interval = 1 * settings.dt
 
     storage = DummyStorage()
-    simulation = Simulation(settings, storage, SpinUp=SpinUp, backend=backend_instance)
+    simulation = Simulation(
+        settings, storage, SpinUp=SpinUp, backend=backend_class(n_dims=2)
+    )
     simulation.reinit()
 
     # Act

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/tests/unit_tests/dynamics/displacement/displacement_settings.py

@@ -25,7 +25,7 @@ class DisplacementSettings:  # pylint: disable=too-few-public-methods,too-many-a
 
     def get_displacement(self, backend, scheme, adaptive=True):
         formulae = Formulae(particle_advection=scheme)
-        backend = backend(formulae, double_precision=True)
+        backend = backend(formulae, n_dims=len(self.grid), double_precision=True)
         environment = DummyEnvironment(
             timestep=self.dt,
             grid=self.grid,

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/tests/unit_tests/dynamics/displacement/test_courant_product.py

@@ -26,7 +26,11 @@ def test_courant_product(courant_field):
     particulator = Particulator(
         n_sd=n_sd,
         environment=Kinematic2D(
-            dt=1, grid=GRID, size=(100, 100), rhod_of=lambda x: x * 0 + 1, backend=CPU()
+            dt=1,
+            grid=GRID,
+            size=(100, 100),
+            rhod_of=lambda x: x * 0 + 1,
+            backend=CPU(n_dims=2),
         ),
         dynamics=(Displacement(),),
         attributes={

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/tests/unit_tests/environments/test_impl.py

@@ -16,8 +16,8 @@ class Env:  # pylint: disable=too-few-public-methods
     def __init__(self, backend):
         self.particulator = None
         self.backend = backend
-        self.mesh = namedtuple("MeshMock", ("grid", "dimension", "n_cell"))(
-            grid=(1, 1), dimension=0, n_cell=1
+        self.mesh = namedtuple("MeshMock", ("grid", "n_dims", "n_cell"))(
+            grid=(1, 1), n_dims=0, n_cell=1
         )
 
     def register(self, *, particulator):

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/tests/unit_tests/exporters/test_vtk_exporter.py

@@ -36,8 +36,8 @@ def test_vtk_exporter_copies_product_data(tmp_path):
             "b": prod,
         },
         dt=0,
-        mesh=namedtuple(typename="MockMesh", field_names=("dimension", "grid", "size"))(
-            dimension=2,
+        mesh=namedtuple(typename="MockMesh", field_names=("n_dims", "grid", "size"))(
+            n_dims=2,
             grid=grid,
             size=(1, 1),
         ),

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/tests/unit_tests/impl/test_mesh.py

@@ -68,13 +68,13 @@ class TestMesh:
 
         # assert
         assert cell_id.shape == (n_sd,)
-        assert cell_origin.shape == (mesh.dimension, n_sd)
-        assert position_in_cell.shape == (mesh.dimension, n_sd)
+        assert cell_origin.shape == (mesh.n_dims, n_sd)
+        assert position_in_cell.shape == (mesh.n_dims, n_sd)
 
         assert 0 <= min(cell_id) <= max(cell_id) < mesh.n_cell
         assert cell_id.dtype == np.int64
 
-        for dim in range(mesh.dimension):
+        for dim in range(mesh.n_dims):
             assert (
                 0
                 <= min(cell_origin[dim, :])
@@ -83,7 +83,7 @@ class TestMesh:
             )
         assert cell_origin.dtype == np.int64
 
-        for dim in range(mesh.dimension):
+        for dim in range(mesh.n_dims):
             assert (
                 0 <= min(position_in_cell[dim, :]) <= max(position_in_cell[dim, :]) < 1
             )

{pysdm-3.0.0rc41 → pysdm-3.0.0rc42}/tests/unit_tests/products/test_cooling_rate.py

@@ -106,7 +106,7 @@ class TestCoolingRate:
             mesh=Mesh(grid=(nz,), size=(z_max,)),
             thd_of_z=lambda z: signed_thd_lapse_rate * z + 300 * si.K,
             rhod_of_z=lambda z: 0 * z + constant_rhod,
-            backend=CPU(),
+            backend=CPU(n_dims=1),
             solvers=namedtuple(typename="_", field_names=("advectee",))(
                 advectee=namedtuple(typename="__", field_names=("ravel", "shape"))(
                     ravel=lambda: None, shape=(nz,)