PySDM 2.73__py3-none-any.whl → 2.76__py3-none-any.whl

PySDM/backends/impl_numba/methods/__init__.py

@@ -12,3 +12,4 @@ from .moments_methods import MomentsMethods
 from .pair_methods import PairMethods
 from .physics_methods import PhysicsMethods
 from .terminal_velocity_methods import TerminalVelocityMethods
+from .seeding_methods import SeedingMethods

PySDM/backends/impl_numba/methods/condensation_methods.py

@@ -445,7 +445,7 @@ class CondensationMethods(BackendMethods):
                 v_drop = formulae.particle_shape_and_density__mass_to_volume(
                     attributes.water_mass[drop]
                 )
-                if v_drop < 0:
+                if v_drop <= 0:
                     continue
                 x_old = formulae.condensation_coordinate__x(v_drop)
                 r_old = formulae.trivia__radius(v_drop)

PySDM/backends/impl_numba/methods/seeding_methods.py

@@ -0,0 +1,68 @@
+""" CPU implementation of backend methods for particle injections """
+
+from functools import cached_property
+
+import numba
+
+from PySDM.backends.impl_common.backend_methods import BackendMethods
+
+
+class SeedingMethods(BackendMethods):  # pylint: disable=too-few-public-methods
+    @cached_property
+    def _seeding(self):
+        @numba.njit(**{**self.default_jit_flags, "parallel": False})
+        def body(  # pylint: disable=too-many-arguments
+            idx,
+            multiplicity,
+            extensive_attributes,
+            seeded_particle_index,
+            seeded_particle_multiplicity,
+            seeded_particle_extensive_attributes,
+            number_of_super_particles_to_inject: int,
+        ):
+            number_of_super_particles_already_injected = 0
+            # TODO #1387 start enumerating from the end of valid particle set
+            for i, mult in enumerate(multiplicity):
+                if (
+                    number_of_super_particles_to_inject
+                    == number_of_super_particles_already_injected
+                ):
+                    break
+                if mult == 0:
+                    idx[i] = -1
+                    s = seeded_particle_index[
+                        number_of_super_particles_already_injected
+                    ]
+                    number_of_super_particles_already_injected += 1
+                    multiplicity[i] = seeded_particle_multiplicity[s]
+                    for a in range(len(extensive_attributes)):
+                        extensive_attributes[a, i] = (
+                            seeded_particle_extensive_attributes[a, s]
+                        )
+            assert (
+                number_of_super_particles_to_inject
+                == number_of_super_particles_already_injected
+            )
+
+        return body
+
+    def seeding(
+        self,
+        *,
+        idx,
+        multiplicity,
+        extensive_attributes,
+        seeded_particle_index,
+        seeded_particle_multiplicity,
+        seeded_particle_extensive_attributes,
+        number_of_super_particles_to_inject: int,
+    ):
+        self._seeding(
+            idx=idx.data,
+            multiplicity=multiplicity.data,
+            extensive_attributes=extensive_attributes.data,
+            seeded_particle_index=seeded_particle_index.data,
+            seeded_particle_multiplicity=seeded_particle_multiplicity.data,
+            seeded_particle_extensive_attributes=seeded_particle_extensive_attributes.data,
+            number_of_super_particles_to_inject=number_of_super_particles_to_inject,
+        )

PySDM/backends/numba.py

@@ -28,6 +28,7 @@ class Numba(  # pylint: disable=too-many-ancestors,duplicate-code
     methods.DisplacementMethods,
     methods.TerminalVelocityMethods,
     methods.IsotopeMethods,
+    methods.SeedingMethods,
 ):
     Storage = ImportedStorage
     Random = ImportedRandom
@@ -75,3 +76,4 @@ class Numba(  # pylint: disable=too-many-ancestors,duplicate-code
         methods.DisplacementMethods.__init__(self)
         methods.TerminalVelocityMethods.__init__(self)
         methods.IsotopeMethods.__init__(self)
+        methods.SeedingMethods.__init__(self)

PySDM/builder.py

@@ -150,4 +150,8 @@ class Builder:
         for key in self.particulator.dynamics:
             self.particulator.timers[key] = WallTimer()
 
+        if (attributes["multiplicity"] == 0).any():
+            self.particulator.attributes.healthy = False
+            self.particulator.attributes.sanitize()
+
         return self.particulator

PySDM/dynamics/__init__.py

@@ -15,3 +15,4 @@ from PySDM.dynamics.displacement import Displacement
 from PySDM.dynamics.eulerian_advection import EulerianAdvection
 from PySDM.dynamics.freezing import Freezing
 from PySDM.dynamics.relaxed_velocity import RelaxedVelocity
+from PySDM.dynamics.seeding import Seeding

PySDM/dynamics/seeding.py

@@ -0,0 +1,94 @@
+""" particle injection handling, requires initalising a simulation with
+enough particles flagged with NaN multiplicity (translated to zeros
+at multiplicity discretisation """
+
+from collections.abc import Sized
+
+import numpy as np
+
+from PySDM.dynamics.impl import register_dynamic
+from PySDM.initialisation import discretise_multiplicities
+
+
+@register_dynamic()
+class Seeding:
+    def __init__(
+        self,
+        *,
+        super_droplet_injection_rate: callable,
+        seeded_particle_extensive_attributes: dict,
+        seeded_particle_multiplicity: Sized,
+    ):
+        for attr in seeded_particle_extensive_attributes.values():
+            assert len(seeded_particle_multiplicity) == len(attr)
+        self.particulator = None
+        self.super_droplet_injection_rate = super_droplet_injection_rate
+        self.seeded_particle_extensive_attributes = seeded_particle_extensive_attributes
+        self.seeded_particle_multiplicity = seeded_particle_multiplicity
+        self.rnd = None
+        self.u01 = None
+        self.index = None
+
+    def register(self, builder):
+        self.particulator = builder.particulator
+
+    def post_register_setup_when_attributes_are_known(self):
+        if tuple(self.particulator.attributes.get_extensive_attribute_keys()) != tuple(
+            self.seeded_particle_extensive_attributes.keys()
+        ):
+            raise ValueError(
+                f"extensive attributes ({self.seeded_particle_extensive_attributes.keys()})"
+                " do not match those used in particulator"
+                f" ({self.particulator.attributes.get_extensive_attribute_keys()})"
+            )
+
+        self.index = self.particulator.Index.identity_index(
+            len(self.seeded_particle_multiplicity)
+        )
+        if len(self.seeded_particle_multiplicity) > 1:
+            self.rnd = self.particulator.Random(
+                len(self.seeded_particle_multiplicity), self.particulator.formulae.seed
+            )
+            self.u01 = self.particulator.Storage.empty(
+                len(self.seeded_particle_multiplicity), dtype=float
+            )
+        self.seeded_particle_multiplicity = (
+            self.particulator.IndexedStorage.from_ndarray(
+                self.index,
+                discretise_multiplicities(
+                    np.asarray(self.seeded_particle_multiplicity)
+                ),
+            )
+        )
+        self.seeded_particle_extensive_attributes = (
+            self.particulator.IndexedStorage.from_ndarray(
+                self.index,
+                np.asarray(list(self.seeded_particle_extensive_attributes.values())),
+            )
+        )
+
+    def __call__(self):
+        if self.particulator.n_steps == 0:
+            self.post_register_setup_when_attributes_are_known()
+
+        time = self.particulator.n_steps * self.particulator.dt
+        number_of_super_particles_to_inject = self.super_droplet_injection_rate(time)
+
+        if number_of_super_particles_to_inject > 0:
+            assert number_of_super_particles_to_inject <= len(
+                self.seeded_particle_multiplicity
+            )
+
+            if self.rnd is not None:
+                self.u01.urand(self.rnd)
+                # TODO #1387 make shuffle smarter
+                # e.g. don't need to shuffle if only one type of seed particle
+                # or if the number of super particles to inject
+                # is equal to the number of possible seeds
+                self.index.shuffle(self.u01)
+            self.particulator.seeding(
+                seeded_particle_index=self.index,
+                number_of_super_particles_to_inject=number_of_super_particles_to_inject,
+                seeded_particle_multiplicity=self.seeded_particle_multiplicity,
+                seeded_particle_extensive_attributes=self.seeded_particle_extensive_attributes,
+            )

PySDM/impl/particle_attributes.py

@@ -118,3 +118,8 @@ class ParticleAttributes:  # pylint: disable=too-many-instance-attributes
 
     def has_attribute(self, attr):
         return attr in self.__attributes
+
+    def reset_idx(self):
+        self.__valid_n_sd = self.__idx.shape[0]
+        self.__idx.reset_index()
+        self.healthy = False

PySDM/initialisation/discretise_multiplicities.py

@@ -6,11 +6,22 @@ import numpy as np
 
 
 def discretise_multiplicities(values_arg):
-    values_int = values_arg.round().astype(np.int64)
+    """any NaN values in the input array are ignored and flagged
+    with zero multiplicities in the output array"""
+
+    values_int = np.where(np.isnan(values_arg), 0, values_arg).round().astype(np.int64)
 
     if np.issubdtype(values_arg.dtype, np.floating):
+        if np.isnan(values_arg).all():
+            return values_int
+
+        if not np.logical_or(values_int > 0, np.isnan(values_arg)).all():
+            raise ValueError(
+                f"int-casting resulted in multiplicity of zero (min(y_float)={min(values_arg)})"
+            )
+
         percent_diff = 100 * abs(
-            1 - np.sum(values_arg) / np.sum(values_int.astype(float))
+            1 - np.nansum(values_arg) / np.sum(values_int.astype(float))
         )
         if percent_diff > 1:
             raise ValueError(
@@ -18,9 +29,4 @@ def discretise_multiplicities(values_arg):
                 f" due to casting multiplicities to ints"
             )
 
-        if not (values_int > 0).all():
-            raise ValueError(
-                f"int-casting resulted in multiplicity of zero (min(y_float)={min(values_arg)})"
-            )
-
     return values_int

PySDM/initialisation/init_fall_momenta.py

@@ -5,7 +5,6 @@ of droplets
 
 import numpy as np
 
-from PySDM.backends import CPU
 from PySDM.dynamics.terminal_velocity import GunnKinzer1949
 from PySDM.formulae import Formulae
 from PySDM.particulator import Particulator
@@ -31,6 +30,8 @@ def init_fall_momenta(
     if zero:
         return np.zeros_like(water_mass)
 
+    from PySDM.backends import CPU  # pylint: disable=import-outside-toplevel
+
     particulator = Particulator(0, CPU(Formulae()))  # TODO #1155
 
     approximation = terminal_velocity_approx(particulator=particulator)

PySDM/particulator.py

@@ -13,7 +13,7 @@ from PySDM.impl.particle_attributes import ParticleAttributes
 
 
 class Particulator:  # pylint: disable=too-many-public-methods,too-many-instance-attributes
-    def __init__(self, n_sd, backend: BackendMethods):
+    def __init__(self, n_sd, backend):
         assert isinstance(backend, BackendMethods)
         self.__n_sd = n_sd
 
@@ -438,3 +438,44 @@ class Particulator:  # pylint: disable=too-many-public-methods,too-many-instance
         self.backend.isotopic_fractionation()
         for isotope in heavy_isotopes:
             self.attributes.mark_updated(f"moles_{isotope}")
+
+    def seeding(
+        self,
+        *,
+        seeded_particle_index,
+        seeded_particle_multiplicity,
+        seeded_particle_extensive_attributes,
+        number_of_super_particles_to_inject,
+    ):
+        n_null = self.n_sd - self.attributes.super_droplet_count
+        if n_null == 0:
+            raise ValueError(
+                "No available seeds to inject. Please provide particles with nan filled attributes."
+            )
+
+        if number_of_super_particles_to_inject > n_null:
+            raise ValueError(
+                "Trying to inject more super particles than space available."
+            )
+
+        if number_of_super_particles_to_inject > len(seeded_particle_multiplicity):
+            raise ValueError(
+                "Trying to inject multiple super particles with the same attributes. \
+                Instead increase multiplicity of injected particles."
+            )
+
+        self.backend.seeding(
+            idx=self.attributes._ParticleAttributes__idx,
+            multiplicity=self.attributes["multiplicity"],
+            extensive_attributes=self.attributes.get_extensive_attribute_storage(),
+            seeded_particle_index=seeded_particle_index,
+            seeded_particle_multiplicity=seeded_particle_multiplicity,
+            seeded_particle_extensive_attributes=seeded_particle_extensive_attributes,
+            number_of_super_particles_to_inject=number_of_super_particles_to_inject,
+        )
+        self.attributes.reset_idx()
+        self.attributes.sanitize()
+
+        self.attributes.mark_updated("multiplicity")
+        for key in self.attributes.get_extensive_attribute_keys():
+            self.attributes.mark_updated(key)

{PySDM-2.73.dist-info → PySDM-2.76.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: PySDM
-Version: 2.73
+Version: 2.76
 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
 Home-page: https://github.com/open-atmos/PySDM
 Author: https://github.com/open-atmos/PySDM/graphs/contributors
@@ -25,7 +25,7 @@ Description-Content-Type: text/markdown
 Requires-Dist: ThrustRTC==0.3.20
 Requires-Dist: CURandRTC>=0.1.2
 Requires-Dist: numba>=0.51.2
-Requires-Dist: numpy<2.0.0
+Requires-Dist: numpy
 Requires-Dist: Pint
 Requires-Dist: chempy
 Requires-Dist: scipy
@@ -39,6 +39,8 @@ Requires-Dist: pyrcel; extra == "tests"
 Requires-Dist: jupyter-core<5.0.0; extra == "tests"
 Requires-Dist: ipywidgets!=8.0.3; extra == "tests"
 
+# <img src="https://raw.githubusercontent.com/open-atmos/PySDM/main/.github/pysdm_logo.png" width=100 height=146 alt="pysdm logo">
+
 # PySDM
 
 [![Python 3](https://img.shields.io/static/v1?label=Python&logo=Python&color=3776AB&message=3)](https://www.python.org/)

{PySDM-2.73.dist-info → PySDM-2.76.dist-info}/RECORD

@@ -1,7 +1,7 @@
 PySDM/__init__.py,sha256=nrBhoz8kmvowBpYj_ijAM85wIqlcCxZg5Y-I1LLpeNo,1738
-PySDM/builder.py,sha256=NpUwIZLXhepjB2BO7NRt6L42vmURXclLnw82iSv671g,5772
+PySDM/builder.py,sha256=6hgxvCFmDe9sJJDT9CgNf9c9ByLvOSvC9cXcc164Yac,5934
 PySDM/formulae.py,sha256=KZmUm6JpVPt5Y8lc7BYc_EWY8ZSMDktFjVHnFENp6YM,13761
-PySDM/particulator.py,sha256=yIRFNMTBDKelLlhD-uejqUL_-QUYu2Y9ShjE5wPzzZw,16447
+PySDM/particulator.py,sha256=pPI2yWD9RsKkiQPKODRtiaBo_cgsMMJH6HwmSJEbj94,18120
 PySDM/attributes/__init__.py,sha256=cJVc2aGo7ZID8fTo_AcxSi0XhnxML-r9yc4L4RYJaE4,171
 PySDM/attributes/chemistry/__init__.py,sha256=sv9kFUXgcZjbjmrQKK8sgdgqM05MPxAw73Euj0O3n40,223
 PySDM/attributes/chemistry/acidity.py,sha256=XGRhQEi7XH9WtBmwzsaMphzuq3ZYGhqClo4PY1d79bA,1642
@@ -47,7 +47,7 @@ PySDM/attributes/physics/terminal_velocity.py,sha256=zu0U0g-qGcx5HJqhTFJHuLyp7lm
 PySDM/attributes/physics/volume.py,sha256=25eG5HEXeRMJV4iKywtZ0JTU5PiF2Ibp3X-WHC6BB6c,587
 PySDM/attributes/physics/water_mass.py,sha256=fzmCjIWLWgghWy_Pa4KZoCd2XgdUv3tb-qB9FFmCUF8,392
 PySDM/backends/__init__.py,sha256=6dPU5Y30rDovvkaTehryG9gA-tJk5sEOrmlxpQAP4Ek,2235
-PySDM/backends/numba.py,sha256=bBtZPmlfoyV1vT2_yfgW_h39IJMC_tOqKlr5nFMKBkw,2708
+PySDM/backends/numba.py,sha256=WaaMBnEpJF0DqjapZfuQu2-DNmJOcfqfss4km1jgRRA,2782
 PySDM/backends/thrust_rtc.py,sha256=5T2joHmLA5y-HPbLfw4m3wf9-BFuk-1YsmWJj62LPNU,2643
 PySDM/backends/impl_common/__init__.py,sha256=hP2LHBlnzcaOh8MMNdm2fbChkfRakwbKn0PaKhLdyP8,37
 PySDM/backends/impl_common/backend_methods.py,sha256=2TaqFIccgQfsNB0X_H4kTgUdgSUu3_H6lDyLUhN4Z2E,560
@@ -66,10 +66,10 @@ PySDM/backends/impl_numba/storage.py,sha256=2ht3sL2FRQcFs-U7373Y2EJ1rHujMZ00rbRp
 PySDM/backends/impl_numba/storage_impl.py,sha256=s1zd6NlpQIQgIlTLvk_GucYUseBouWEGDgvdVHm3LfQ,1957
 PySDM/backends/impl_numba/toms748.py,sha256=pcI0Fhocma_9jbfkW3oDemWikDhG1pOWF3iamn2R2eU,6393
 PySDM/backends/impl_numba/warnings.py,sha256=MiYn_5sqWYV__dkVVQu88X_XcYW1N5vqDZtyE5c8FmA,649
-PySDM/backends/impl_numba/methods/__init__.py,sha256=T9Ex5Ra_G5_Xr-YTlt3pGPrlTBU8ocxHfQr4GVCv8Vg,624
+PySDM/backends/impl_numba/methods/__init__.py,sha256=esVAXeSXLDLqrue--vFac0zMcUsL6bUc1q6Q8SnyeQo,668
 PySDM/backends/impl_numba/methods/chemistry_methods.py,sha256=qAz3VTZAPQeTuti11w6Wv1ZZUXSHa4jxkkihukgLx0k,15508
 PySDM/backends/impl_numba/methods/collisions_methods.py,sha256=PkCnV3apuC16mSxel733kkWsgXRTlQf_WvbkLzRwmzY,25941
-PySDM/backends/impl_numba/methods/condensation_methods.py,sha256=Uu1jhRVtqWXzIibV5K6ud6bSH7sra2ylBUPIjb59ZXk,24985
+PySDM/backends/impl_numba/methods/condensation_methods.py,sha256=Gu4c6MN17WKtSSesgvhO2SGP6xLPQItJ1qjpPUxtSXI,24986
 PySDM/backends/impl_numba/methods/displacement_methods.py,sha256=pgsf98FQNf-sdHx4aiSLwn8m8JKk3eG0J_DlsDC8hLo,7673
 PySDM/backends/impl_numba/methods/fragmentation_methods.py,sha256=LzaBL7DxwGoprOMz7t-vRrpnKLl-lq_mpTbAyoV_VxA,17127
 PySDM/backends/impl_numba/methods/freezing_methods.py,sha256=-ZKoKzIsMMuk6eL3QqNaZ9kd8QnBr3AEgcGCCoDUrQM,5019
@@ -78,6 +78,7 @@ PySDM/backends/impl_numba/methods/isotope_methods.py,sha256=aFdezr4MqneCxYJVXDT9
 PySDM/backends/impl_numba/methods/moments_methods.py,sha256=7UzggKu186vqouByHNcgE03QkV0OwJn4f8RpWtVfFKM,5512
 PySDM/backends/impl_numba/methods/pair_methods.py,sha256=NeRji-X1NBUFDLGsGyGYFUWkFTKX9JIZp_Bny30AVJ0,6259
 PySDM/backends/impl_numba/methods/physics_methods.py,sha256=aM3LaLsx06nH8_l9CDYoqMOiCRKt4D6MAl322kCYgRM,7053
+PySDM/backends/impl_numba/methods/seeding_methods.py,sha256=vZDmQ7XEVX8CivdMT3tmmHpOs9zSACYxL1mip-B_iSw,2571
 PySDM/backends/impl_numba/methods/terminal_velocity_methods.py,sha256=KO-RfmseucWEXG-i2CGNPpqyFYXzmOExNQdw4JxrGJo,2117
 PySDM/backends/impl_numba/test_helpers/__init__.py,sha256=8_cw4HYZE6jatsdpVWJHbbS2HwXWisbMpjFsos6JBqk,94
 PySDM/backends/impl_numba/test_helpers/scipy_ode_condensation_solver.py,sha256=MRnpnhEgaZjBdCMTEVHmOJ_Z2V2bLoLnXmjNDPYKdYM,10250
@@ -103,7 +104,7 @@ PySDM/backends/impl_thrust_rtc/test_helpers/__init__.py,sha256=c-QZ2cI70jcZ8Yvb1
 PySDM/backends/impl_thrust_rtc/test_helpers/cpp2python.py,sha256=akXMMjRPTXLhvTPmifFkPcbIdBjsAzSMby7QWDiur_s,6078
 PySDM/backends/impl_thrust_rtc/test_helpers/fake_thrust_rtc.py,sha256=LMzL3KtLtVJRfjV4PxbJz7syo46EFNPkJ2HGZm3p6AQ,7962
 PySDM/backends/impl_thrust_rtc/test_helpers/flag.py,sha256=jVDnnq0rN_2jXlgbXx6nDY6m9jncV2CwjdL3-tKZ9iY,170
-PySDM/dynamics/__init__.py,sha256=wl-gnX3q9pGHjD5javTNVSigxzOdeVqcZytECqMnEFw,721
+PySDM/dynamics/__init__.py,sha256=COm5PdJoSCxgwWiwMFFaXBlxLf0wYEUfw9jQOGSHcrM,764
 PySDM/dynamics/ambient_thermodynamics.py,sha256=GWRzhyVZLNZjneHJ7aFk7LK1UBaUs48XfxzqkfX8jj8,353
 PySDM/dynamics/aqueous_chemistry.py,sha256=zk4wwf-TPkHn5BVitR4UcpcU3aH6sUuEpGOviRcL9DI,4823
 PySDM/dynamics/condensation.py,sha256=n1Mm0yqBE-cZHigth-PLgekLIOYQyxXvX81a_NnpqR0,4361
@@ -112,6 +113,7 @@ PySDM/dynamics/eulerian_advection.py,sha256=NOPO6h0V9ApjeGNxxE0_86AzBUY27DL8Sinr
 PySDM/dynamics/freezing.py,sha256=iXiHZJOmK3RarTI7Vm5VF7fO4OosR9hd_Dz-9-R8dzA,3710
 PySDM/dynamics/isotopic_fractionation.py,sha256=xJoEJzuA-93Jyww_BQxduZqDNRrsUD95bLzvhgI6Q1k,1270
 PySDM/dynamics/relaxed_velocity.py,sha256=SCLVcqJIsFvzQ5LDZ7YxTSFAuojQf-sbiU-_tY_wteo,3131
+PySDM/dynamics/seeding.py,sha256=OBofNtF9M5-j7BoAAAAayvAlATx2CIZu56StVIOO76s,3842
 PySDM/dynamics/collisions/__init__.py,sha256=gutS6RlEXSp9S7z6zgUKq-K4hXLrHXTJFQkYCi2KibY,554
 PySDM/dynamics/collisions/collision.py,sha256=Wl-hSj3BamAZLY7fg2bQ7_hMAoucMuvEmZqvR3oAWno,13359
 PySDM/dynamics/collisions/breakup_efficiencies/__init__.py,sha256=PWWgvr3wLjlOKjiK_S9-Bd9bcGFIIwMlbwJtMriDc1Q,59
@@ -174,13 +176,13 @@ PySDM/impl/arakawa_c.py,sha256=ACzN8k4lVapTP5ra2r3foUO-olBJf5ZcFTpH5I4SCEE,331
 PySDM/impl/camel_case.py,sha256=2ivX9gjsLDPwYvGWL0omPqSBp9LXyg-rPLC_pJ62pTU,367
 PySDM/impl/mesh.py,sha256=h64y8vKLR2sqKcs7BMSrtlz2S6JIq15KJMW4KQfKynI,2528
 PySDM/impl/null_physics_class.py,sha256=-AEqs_yrPyDkgOlkXtG4nMtIeyp6tShiyVNvQqQexkI,128
-PySDM/impl/particle_attributes.py,sha256=g0LbhZ3cM8kLWOG07acFGTAmZ2qNuzYBV2nWBFWBENo,3566
+PySDM/impl/particle_attributes.py,sha256=97VK7EMyupc_FouzZGcD7TTS9bakWc4vIbjfewuMGlo,3702
 PySDM/impl/particle_attributes_factory.py,sha256=ZKkkbkFIoEcw7Yj7axlmrsqhy-_Hb9TusMnu1XtuH3w,4687
 PySDM/impl/wall_timer.py,sha256=cDApBqoTS_kUfnBzzQtP-XTJxBMBmyqv4m9mODaA2CE,465
 PySDM/initialisation/__init__.py,sha256=Li4IboyzZBIwoLbNKceZq9nDUgJpShQr8QyibDZM_5Q,351
-PySDM/initialisation/discretise_multiplicities.py,sha256=D_8EsoXwSLsgIGUyOuLDDVuaQpZAQbtEDtJfcXK90Bw,779
+PySDM/initialisation/discretise_multiplicities.py,sha256=rRUJUGQYCUeLbVoUonrCvK3_w8HcJz8bIP8FfBzUBlI,1040
 PySDM/initialisation/equilibrate_wet_radii.py,sha256=mQSMWChcF3YnfH9AIYAzYgjZXmXxZ3NaS0yNVCk5nuY,3774
-PySDM/initialisation/init_fall_momenta.py,sha256=LFMTVm5DaO8B2_txi4yj9jvHGpgNX3JxYXzdaCai_VE,1427
+PySDM/initialisation/init_fall_momenta.py,sha256=RlZSFtMgIRXAiimOfzY3Iu9GlUTxaIn0wClRuBIcrbA,1475
 PySDM/initialisation/aerosol_composition/__init__.py,sha256=IsaPAyhxPXHaXY2nu9cyYdtjWnwxAGKZejvpRsv3Dyc,206
 PySDM/initialisation/aerosol_composition/dry_aerosol.py,sha256=FxFqPyR7Fhg1mi5lYfim4Z4PsnWrwIKGsVRiWa2HW6s,4372
 PySDM/initialisation/impl/__init__.py,sha256=vejuF9ebIGXb7l9zEiBBbNjs1Qct1KZuZRPLh8w3SJE,125
@@ -381,8 +383,8 @@ PySDM/products/size_spectral/size_standard_deviation.py,sha256=A7fSj_keQi_bHW41P
 PySDM/products/size_spectral/total_particle_concentration.py,sha256=K7sO2jgUpYOwWvBK3jXp9KUxLRXk_Ve6ooOdY5q5jd0,470
 PySDM/products/size_spectral/total_particle_specific_concentration.py,sha256=isLZt5j2vJGr-FrzYFcUNPHZqMKhOtl51za1Hogh6hQ,476
 PySDM/products/size_spectral/water_mixing_ratio.py,sha256=PP54hfMy4SPDQDJgV4RNzVP03lmL1THZPa2KXrdgoQ4,1324
-PySDM-2.73.dist-info/LICENSE,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
-PySDM-2.73.dist-info/METADATA,sha256=wEzcfKjLwxOYHbjaYkvoAaLf4q83agwXc-VM18naydw,45910
-PySDM-2.73.dist-info/WHEEL,sha256=eOLhNAGa2EW3wWl_TU484h7q1UNgy0JXjjoqKoxAAQc,92
-PySDM-2.73.dist-info/top_level.txt,sha256=6NYaDAhUnaYMD7ebD_dQ7HlhiJS7NMuKEXb0ck7Wcsg,6
-PySDM-2.73.dist-info/RECORD,,
+PySDM-2.76.dist-info/LICENSE,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
+PySDM-2.76.dist-info/METADATA,sha256=t9WrHc4GBaUdP8_MmzJ8rDQZAGgAOfACj5UdCy4lmaQ,46036
+PySDM-2.76.dist-info/WHEEL,sha256=eOLhNAGa2EW3wWl_TU484h7q1UNgy0JXjjoqKoxAAQc,92
+PySDM-2.76.dist-info/top_level.txt,sha256=6NYaDAhUnaYMD7ebD_dQ7HlhiJS7NMuKEXb0ck7Wcsg,6
+PySDM-2.76.dist-info/RECORD,,