PySDM 2.68__tar.gz → 2.70__tar.gz

{PySDM-2.68 → PySDM-2.70}/.codecov.yml

@@ -2,4 +2,4 @@ coverage:
   status:
     project:
       default:
-        threshold: 0.01%
+        threshold: 0.1%

{PySDM-2.68 → PySDM-2.70}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: PySDM
-Version: 2.68
+Version: 2.70
 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

{PySDM-2.68 → PySDM-2.70}/PySDM/builder.py

@@ -61,9 +61,9 @@ class Builder:
     def _register_product(self, product, buffer):
         if product.name in self.particulator.products:
             raise ValueError(f'product name "{product.name}" already registered')
-        product.set_buffer(buffer)
-        product.register(self)
-        self.particulator.products[product.name] = product
+        self.particulator.products[product.name] = product.instantiate(
+            builder=self, buffer=buffer
+        )
 
     def _resolve_attribute(self, attr_name):
         if attr_name not in self.req_attr:

{PySDM-2.68 → PySDM-2.70}/PySDM/formulae.py

@@ -119,7 +119,7 @@ class Formulae:  # pylint: disable=too-few-public-methods,too-many-instance-attr
             **constants_defaults
         )
         self.constants = constants
-        self.seed = seed or physics.constants.default_random_seed
+        self.seed = seed if seed is not None else physics.constants.default_random_seed
         self.fastmath = fastmath
         self.handle_all_breakups = handle_all_breakups
         dimensional_analysis = physics.impl.flag.DIMENSIONAL_ANALYSIS

{PySDM-2.68 → PySDM-2.70}/PySDM/impl/particle_attributes.py

@@ -22,7 +22,6 @@ class ParticleAttributes:  # pylint: disable=too-many-instance-attributes
         attributes: Dict[str, Attribute],
     ):
         self.__valid_n_sd = particulator.n_sd
-        self.healthy = True
         self.__healthy_memory = particulator.Storage.from_ndarray(np.full((1,), 1))
         self.__idx = idx
 
@@ -40,6 +39,14 @@ class ParticleAttributes:  # pylint: disable=too-many-instance-attributes
         self.__sorted = False
         self.__attributes = attributes
 
+    @property
+    def healthy(self) -> bool:
+        return bool(self.__healthy_memory[0])
+
+    @healthy.setter
+    def healthy(self, value: bool):
+        self.__healthy_memory[:] = value
+
     @property
     def cell_start(self):
         if not self.__sorted:
@@ -63,7 +70,6 @@ class ParticleAttributes:  # pylint: disable=too-many-instance-attributes
             self.__idx.remove_zero_n_or_flagged(self["multiplicity"])
             self.__valid_n_sd = self.__idx.length
             self.healthy = True
-            self.__healthy_memory[:] = 1
             self.__sorted = False
 
     def cut_working_length(self, length):

{PySDM-2.68 → PySDM-2.70}/PySDM/particulator.py

@@ -202,9 +202,6 @@ class Particulator:  # pylint: disable=too-many-public-methods,too-many-instance
                 coalescence_rate=coalescence_rate,
                 is_first_in_pair=is_first_in_pair,
             )
-        self.attributes.healthy = bool(
-            self.attributes._ParticleAttributes__healthy_memory
-        )
         self.attributes.sanitize()
         self.attributes.mark_updated("multiplicity")
         for key in self.attributes.get_extensive_attribute_keys():
@@ -410,9 +407,6 @@ class Particulator:  # pylint: disable=too-many-public-methods,too-many-instance
             precipitation_counting_level_index=precipitation_counting_level_index,
             displacement=displacement,
         )
-        self.attributes.healthy = bool(
-            self.attributes._ParticleAttributes__healthy_memory
-        )
         self.attributes.sanitize()
         return res
 
@@ -425,9 +419,6 @@ class Particulator:  # pylint: disable=too-many-public-methods,too-many-instance
             healthy=self.attributes._ParticleAttributes__healthy_memory,
             domain_top_level_index=self.mesh.grid[-1],
         )
-        self.attributes.healthy = bool(
-            self.attributes._ParticleAttributes__healthy_memory
-        )
         self.attributes.sanitize()
 
     def calculate_displacement(

{PySDM-2.68 → PySDM-2.70}/PySDM/products/ambient_thermodynamics/ambient_dry_air_density.py

@@ -2,9 +2,10 @@
 ambient dry-air density
 """
 
-from PySDM.products.impl.moist_environment_product import MoistEnvironmentProduct
+from PySDM.products.impl import MoistEnvironmentProduct, register_product
 
 
+@register_product()
 class AmbientDryAirDensity(MoistEnvironmentProduct):
     def __init__(self, name="rhod", unit="kg/m^3", var=None):
         super().__init__(name=name, unit=unit, var=var)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/ambient_thermodynamics/ambient_dry_air_potential_temperature.py

@@ -2,9 +2,10 @@
 ambient dry-air potential temperature (computed using dry air partial pressure)
 """
 
-from PySDM.products.impl.moist_environment_product import MoistEnvironmentProduct
+from PySDM.products.impl import MoistEnvironmentProduct, register_product
 
 
+@register_product()
 class AmbientDryAirPotentialTemperature(MoistEnvironmentProduct):
     def __init__(self, unit="K", name=None, var=None):
         super().__init__(unit=unit, name=name, var=var)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/ambient_thermodynamics/ambient_pressure.py

@@ -2,9 +2,10 @@
 ambient pressure
 """
 
-from PySDM.products.impl.moist_environment_product import MoistEnvironmentProduct
+from PySDM.products.impl import MoistEnvironmentProduct, register_product
 
 
+@register_product()
 class AmbientPressure(MoistEnvironmentProduct):
     def __init__(self, name=None, unit="Pa", var=None):
         super().__init__(name=name, unit=unit, var=var)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/ambient_thermodynamics/ambient_relative_humidity.py

@@ -2,9 +2,10 @@
 ambient relative humidity (wrt water or ice)
 """
 
-from PySDM.products.impl.moist_environment_product import MoistEnvironmentProduct
+from PySDM.products.impl import MoistEnvironmentProduct, register_product
 
 
+@register_product()
 class AmbientRelativeHumidity(MoistEnvironmentProduct):
     def __init__(self, name=None, unit="dimensionless", var=None, ice=False):
         super().__init__(name=name, unit=unit, var=var)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/ambient_thermodynamics/ambient_temperature.py

@@ -2,9 +2,10 @@
 ambient temperature
 """
 
-from PySDM.products.impl.moist_environment_product import MoistEnvironmentProduct
+from PySDM.products.impl import MoistEnvironmentProduct, register_product
 
 
+@register_product()
 class AmbientTemperature(MoistEnvironmentProduct):
     def __init__(self, name=None, unit="K", var=None):
         super().__init__(name=name, unit=unit, var=var)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/ambient_thermodynamics/ambient_water_vapour_mixing_ratio.py

@@ -2,9 +2,10 @@
 ambient water vapour mixing ratio (mass of water vapour per mass of dry air)
 """
 
-from PySDM.products.impl.moist_environment_product import MoistEnvironmentProduct
+from PySDM.products.impl import MoistEnvironmentProduct, register_product
 
 
+@register_product()
 class AmbientWaterVapourMixingRatio(MoistEnvironmentProduct):
     def __init__(self, name=None, unit="dimensionless", var=None):
         super().__init__(unit=unit, name=name, var=var)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/aqueous_chemistry/acidity.py

@@ -5,9 +5,10 @@ with weighting either by number or volume
 
 import numpy as np
 
-from PySDM.products.impl.moment_product import MomentProduct
+from PySDM.products.impl import MomentProduct, register_product
 
 
+@register_product()
 class Acidity(MomentProduct):
     def __init__(
         self,

{PySDM-2.68 → PySDM-2.70}/PySDM/products/aqueous_chemistry/aqueous_mass_spectrum.py

@@ -8,9 +8,10 @@ from chempy import Substance
 
 from PySDM.dynamics.impl.chemistry_utils import AQUEOUS_COMPOUNDS
 from PySDM.physics.constants import si
-from PySDM.products.impl.spectrum_moment_product import SpectrumMomentProduct
+from PySDM.products.impl import SpectrumMomentProduct, register_product
 
 
+@register_product()
 class AqueousMassSpectrum(SpectrumMomentProduct):
     def __init__(
         self, *, key, dry_radius_bins_edges, specific=False, name=None, unit="kg/m^3"
@@ -58,6 +59,7 @@ class AqueousMassSpectrum(SpectrumMomentProduct):
         return vals
 
 
+@register_product()
 class SpecificAqueousMassSpectrum(AqueousMassSpectrum):
     def __init__(self, key, dry_radius_bins_edges, name=None, unit="dimensionless"):
         super().__init__(

{PySDM-2.68 → PySDM-2.70}/PySDM/products/aqueous_chemistry/aqueous_mole_fraction.py

@@ -2,11 +2,12 @@
 mole fractions of aqueous-chemistry relevant compounds
 """
 
-from PySDM.products.impl.moment_product import MomentProduct
+from PySDM.products.impl import MomentProduct, register_product
 
 DUMMY_SPECIFIC_GRAVITY = 44
 
 
+@register_product()
 class AqueousMoleFraction(MomentProduct):
     def __init__(self, key, unit="dimensionless", name=None):
         super().__init__(unit=unit, name=name)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/aqueous_chemistry/gaseous_mole_fraction.py

@@ -3,9 +3,10 @@ mole fractions of gaseous compounds relevant for aqueous chemistry
 """
 
 from PySDM.dynamics.impl.chemistry_utils import GASEOUS_COMPOUNDS
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class GaseousMoleFraction(Product):
     def __init__(self, key, unit="dimensionless", name=None):
         super().__init__(name=name, unit=unit)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/aqueous_chemistry/total_dry_mass_mixing_ratio.py

@@ -4,9 +4,10 @@ dry aerosol mass summed over all particles in a grid box per mass of dry air
 
 import numpy as np
 
-from PySDM.products.impl.moment_product import MomentProduct
+from PySDM.products.impl import MomentProduct, register_product
 
 
+@register_product()
 class TotalDryMassMixingRatio(MomentProduct):
     def __init__(self, density, name=None, unit="kg/kg"):
         super().__init__(unit=unit, name=name)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/collision/collision_rates.py

@@ -4,9 +4,10 @@ Rates of collision events and their deficit wrt expected values in case mutiplic
  or due to unrepresentable breakups from integer overflow
 """
 
-from PySDM.products.impl.rate_product import RateProduct
+from PySDM.products.impl import RateProduct, register_product
 
 
+@register_product()
 class CollisionRateDeficitPerGridbox(RateProduct):
     def __init__(self, name=None, unit="s^-1 kg^-1"):
         super().__init__(
@@ -14,6 +15,7 @@ class CollisionRateDeficitPerGridbox(RateProduct):
         )
 
 
+@register_product()
 class CollisionRatePerGridbox(RateProduct):
     def __init__(self, name=None, unit="s^-1 kg^-1"):
         super().__init__(
@@ -21,6 +23,7 @@ class CollisionRatePerGridbox(RateProduct):
         )
 
 
+@register_product()
 class CoalescenceRatePerGridbox(RateProduct):
     def __init__(self, name=None, unit="s^-1 kg^-1"):
         super().__init__(
@@ -28,6 +31,7 @@ class CoalescenceRatePerGridbox(RateProduct):
         )
 
 
+@register_product()
 class BreakupRatePerGridbox(RateProduct):
     def __init__(self, name=None, unit="s^-1 kg^-1"):
         super().__init__(
@@ -35,6 +39,7 @@ class BreakupRatePerGridbox(RateProduct):
         )
 
 
+@register_product()
 class BreakupRateDeficitPerGridbox(RateProduct):
     def __init__(self, name=None, unit="s^-1 kg^-1"):
         super().__init__(

{PySDM-2.68 → PySDM-2.70}/PySDM/products/collision/collision_timestep_mean.py

@@ -7,9 +7,10 @@ import numba
 import numpy as np
 
 from PySDM.backends.impl_numba.conf import JIT_FLAGS
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class CollisionTimestepMean(Product):
     def __init__(self, unit="s", name=None):
         super().__init__(unit=unit, name=name)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/collision/collision_timestep_min.py

@@ -3,9 +3,10 @@ Minimal collision timestep used when adaptive timestepping is enabled in the
  `PySDM.dynamics.collisions.collision.Collision` dynamic (fetching a value resets the counter)
 """
 
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class CollisionTimestepMin(Product):
     def __init__(self, unit="s", name=None):
         super().__init__(unit=unit, name=name)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/condensation/activable_fraction.py

@@ -3,9 +3,10 @@ fraction of particles with critical supersaturation lower than a given supersatu
  (passed as keyword argument while calling `get()`)
 """
 
-from PySDM.products.impl.moment_product import MomentProduct
+from PySDM.products.impl import MomentProduct, register_product
 
 
+@register_product()
 class ActivableFraction(MomentProduct):
     def __init__(self, unit="dimensionless", name=None):
         super().__init__(name=name, unit=unit)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/condensation/condensation_timestep.py

@@ -4,9 +4,10 @@ minimum and maximum condensation timestep (fetching a value resets the counter)
 
 import numpy as np
 
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class _CondensationTimestep(Product):
     def __init__(self, name, unit, extremum, reset_value):
         super().__init__(
@@ -37,11 +38,13 @@ class _CondensationTimestep(Product):
         return self.buffer
 
 
+@register_product()
 class CondensationTimestepMin(_CondensationTimestep):
     def __init__(self, name=None, unit="s"):
         super().__init__(name=name, unit=unit, extremum=np.minimum, reset_value=np.inf)
 
 
+@register_product()
 class CondensationTimestepMax(_CondensationTimestep):
     def __init__(self, name=None, unit="s"):
         super().__init__(name=name, unit=unit, extremum=np.maximum, reset_value=-np.inf)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/condensation/event_rates.py

@@ -5,9 +5,10 @@ rates of activation, deactivation and ripening events (take into account substep
 
 import numpy as np
 
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class EventRate(Product):
     def __init__(self, what, name=None, unit=None):
         super().__init__(name=name, unit=unit)
@@ -42,16 +43,19 @@ class EventRate(Product):
         return self.buffer
 
 
+@register_product()
 class RipeningRate(EventRate):
     def __init__(self, name=None, unit="s^-1 kg^-1"):
         super().__init__("ripening", name=name, unit=unit)
 
 
+@register_product()
 class ActivatingRate(EventRate):
     def __init__(self, name=None, unit="s^-1 kg^-1"):
         super().__init__("activating", name=name, unit=unit)
 
 
+@register_product()
 class DeactivatingRate(EventRate):
     def __init__(self, name=None, unit="s^-1 kg^-1"):
         super().__init__("deactivating", name=name, unit=unit)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/condensation/peak_supersaturation.py

@@ -7,9 +7,10 @@ highest supersaturation encountered while solving for condensation/evaporation (
 
 import numpy as np
 
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class PeakSupersaturation(Product):
     def __init__(self, unit="dimensionless", name=None):
         super().__init__(unit=unit, name=name)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/displacement/averaged_terminal_velocity.py

@@ -4,9 +4,10 @@ average terminal with weighting either by number or volume
 
 import numpy as np
 
-from PySDM.products.impl.moment_product import MomentProduct
+from PySDM.products.impl import MomentProduct, register_product
 
 
+@register_product()
 class AveragedTerminalVelocity(MomentProduct):
     def __init__(
         self, *, radius_range=(0, np.inf), weighting="volume", unit="m/s", name=None

{PySDM-2.68 → PySDM-2.70}/PySDM/products/displacement/flow_velocity_component.py

@@ -5,9 +5,10 @@ absolute values of Courant number on all edges of a cell)
 
 import numpy as np
 
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class FlowVelocityComponent(Product):
     def __init__(self, component: int, name=None, unit="m/s"):
         super().__init__(unit=unit, name=name)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/displacement/max_courant_number.py

@@ -5,9 +5,10 @@ absolute values of Courant number on all edges of a cell)
 
 import numpy as np
 
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class MaxCourantNumber(Product):
     def __init__(self, name=None, unit="dimensionless"):
         super().__init__(unit=unit, name=name)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/displacement/surface_precipitation.py

@@ -2,9 +2,10 @@
 water volume flux derived from sizes of particles crossing bottom domain boundary
 """
 
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class SurfacePrecipitation(Product):
     def __init__(self, name=None, unit="m/s"):
         super().__init__(unit=unit, name=name)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/freezing/cooling_rate.py

@@ -2,9 +2,10 @@
 number-averaged cooling rate
 """
 
-from PySDM.products.impl import MomentProduct
+from PySDM.products.impl import MomentProduct, register_product
 
 
+@register_product()
 class CoolingRate(MomentProduct):
     def __init__(self, unit="K/s", name=None):
         super().__init__(unit=unit, name=name)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/freezing/freezable_specific_concentration.py

@@ -4,9 +4,10 @@ freezing-temperature binned specific concentration of particles
 
 import numpy as np
 
-from PySDM.products.impl.spectrum_moment_product import SpectrumMomentProduct
+from PySDM.products.impl import SpectrumMomentProduct, register_product
 
 
+@register_product()
 class FreezableSpecificConcentration(SpectrumMomentProduct):
     def __init__(self, temperature_bins_edges, name=None, unit="kg^-1 K^-1"):
         super().__init__(name=name, unit=unit, attr_unit="K")

{PySDM-2.68 → PySDM-2.70}/PySDM/products/freezing/frozen_particle_concentration.py

@@ -4,9 +4,10 @@ concentration of frozen particles (unactivated, activated or both)
 
 import numpy as np
 
-from PySDM.products.impl.concentration_product import ConcentrationProduct
+from PySDM.products.impl import ConcentrationProduct, register_product
 
 
+@register_product()
 class FrozenParticleConcentration(ConcentrationProduct):
     def __init__(
         self,
@@ -39,6 +40,7 @@ class FrozenParticleConcentration(ConcentrationProduct):
         return super()._impl(**kwargs)
 
 
+@register_product()
 class FrozenParticleSpecificConcentration(FrozenParticleConcentration):
     def __init__(self, *, count_unactivated, count_activated, unit="kg^-1", name=None):
         super().__init__(

{PySDM-2.68 → PySDM-2.70}/PySDM/products/freezing/ice_nuclei_concentration.py

@@ -4,9 +4,10 @@ immersed ice nucleus concentration (both within frozen and unfrozen particles)
 
 import numpy as np
 
-from PySDM.products.impl.concentration_product import ConcentrationProduct
+from PySDM.products.impl import ConcentrationProduct, register_product
 
 
+@register_product()
 class IceNucleiConcentration(ConcentrationProduct):
     def __init__(self, unit="m^-3", name=None, __specific=False, stp=False):
         super().__init__(unit=unit, name=name, specific=__specific, stp=stp)
@@ -34,6 +35,7 @@ class IceNucleiConcentration(ConcentrationProduct):
         return super()._impl(**kwargs)
 
 
+@register_product()
 class SpecificIceNucleiConcentration(IceNucleiConcentration):
     def __init__(self, unit="kg^-1", name=None, __specific=True):
         super().__init__(unit=unit, name=name)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/freezing/total_unfrozen_immersed_surface_area.py

@@ -4,9 +4,10 @@ mean immersed surface area per particle for unfrozen particles
 
 import numpy as np
 
-from PySDM.products.impl.moment_product import MomentProduct
+from PySDM.products.impl import MomentProduct, register_product
 
 
+@register_product()
 class TotalUnfrozenImmersedSurfaceArea(MomentProduct):
     def __init__(self, unit="m^2", name=None):
         super().__init__(unit=unit, name=name)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/housekeeping/dynamic_wall_time.py

@@ -2,9 +2,10 @@
 wall-time for a given dynamic (fetching a value resets the counter)
 """
 
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class DynamicWallTime(Product):
     def __init__(self, dynamic, name=None, unit="s"):
         super().__init__(name=name, unit=unit)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/housekeeping/super_droplet_count_per_gridbox.py

@@ -5,9 +5,10 @@ super-droplet count per gridbox (dimensionless)
 import numba
 
 from PySDM.backends.impl_numba.conf import JIT_FLAGS
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class SuperDropletCountPerGridbox(Product):
     def __init__(self, unit="dimensionless", name=None):
         super().__init__(unit=unit, name=name)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/housekeeping/time.py

@@ -2,9 +2,10 @@
 physical time (in dt increments)
 """
 
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class Time(Product):
     def __init__(self, name=None, unit="s"):
         super().__init__(name=name, unit=unit)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/housekeeping/timers.py

@@ -5,7 +5,7 @@ CPU- and wall-time counters (fetching a value resets the counter)
 import time
 from abc import abstractmethod
 
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
 class _Timer(Product):
@@ -33,6 +33,7 @@ class _Timer(Product):
         raise NotImplementedError()
 
 
+@register_product()
 class CPUTime(_Timer):
     def __init__(self, name="CPU Time", unit="s"):
         super().__init__(unit=unit, name=name)
@@ -42,6 +43,7 @@ class CPUTime(_Timer):
         return time.process_time()
 
 
+@register_product()
 class WallTime(_Timer):
     def __init__(self, name=None, unit="s"):
         super().__init__(unit=unit, name=name)

{PySDM-2.68 → PySDM-2.70}/PySDM/products/impl/__init__.py

@@ -5,3 +5,6 @@ from PySDM.products.impl.moment_product import MomentProduct
 from PySDM.products.impl.product import Product
 from PySDM.products.impl.rate_product import RateProduct
 from PySDM.products.impl.spectrum_moment_product import SpectrumMomentProduct
+from PySDM.products.impl.concentration_product import ConcentrationProduct
+from PySDM.products.impl.activation_filtered_product import ActivationFilteredProduct
+from PySDM.products.impl.register_product import register_product

{PySDM-2.68 → PySDM-2.70}/PySDM/products/impl/activation_filtered_product.py

@@ -5,7 +5,7 @@ common base class for products filtering droplets based on their activation stat
 import numpy as np
 
 
-class _ActivationFilteredProduct:
+class ActivationFilteredProduct:
     def __init__(
         self,
         *,

{PySDM-2.68 → PySDM-2.70}/PySDM/products/impl/product.py

@@ -60,7 +60,11 @@ class Product:
 
         default_unit_arg = init.parameters["unit"].default
 
-        if default_unit_arg is None or str(default_unit_arg).strip() == "":
+        if (
+            default_unit_arg is inspect._empty
+            or default_unit_arg is None
+            or str(default_unit_arg).strip() == ""
+        ):
             raise AssertionError(
                 f"unit parameter of {type(self).__name__}.__init__"
                 f" is expected to have a non-empty default value"