essreduce 25.2.4__tar.gz → 25.2.6__tar.gz

{essreduce-25.2.4/src/essreduce.egg-info → essreduce-25.2.6}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: essreduce
-Version: 25.2.4
+Version: 25.2.6
 Summary: Common data reduction tools for the ESS facility
 Author: Scipp contributors
 License: BSD 3-Clause License

{essreduce-25.2.4 → essreduce-25.2.6}/src/ess/reduce/live/raw.py

@@ -15,14 +15,12 @@ options:
 - `'xy_plane'`: Project the data onto the x-y plane, i.e., perpendicular to the beam.
 - `'cylinder_mantle_z'`: Project the data onto the mantle of a cylinder aligned with the
    z-axis.
-- `LogicalView`: Not a projection in the traditional sense, but a way to select and
-  flatten dimensions of the data.
+- A callable, e.g., to select and flatten dimensions of the data.
 """
 
 from __future__ import annotations
 
 from collections.abc import Callable, Sequence
-from dataclasses import dataclass, field
 from math import ceil
 from typing import Literal, NewType
 
@@ -140,44 +138,6 @@ class Histogrammer:
         return self._hist(replicated, coords=self._coords) / self.replicas
 
 
-@dataclass
-class LogicalView:
-    """
-    Logical view of a multi-dimensional detector.
-
-    Instances can be used as a "projection" function for a detector view.
-
-    Parameters
-    ----------
-    fold:
-        Dimensions to fold. This is useful is the raw data has a single dimension that
-        corresponds to multiple dimensions in the logical view.
-    transpose:
-        Dimensions to transpose. This is useful for reordering dimensions.
-    select:
-        Dimensions with associated index to select from the data. This extracts a slice
-        of the data for each given dimension.
-    flatten:
-        Dimensions to flatten.
-    """
-
-    fold: dict[str, int] | None = None
-    transpose: tuple[str, ...] | None = None
-    select: dict[str, int] = field(default_factory=dict)
-    flatten: dict[str, list[str]] = field(default_factory=dict)
-
-    def __call__(self, da: sc.DataArray) -> sc.DataArray:
-        if self.fold is not None:
-            da = da.fold(da.dim, sizes=self.fold)
-        if self.transpose is not None:
-            da = da.transpose(self.transpose)
-        for dim, index in self.select.items():
-            da = da[dim, index]
-        for to, dims in self.flatten.items():
-            da = da.flatten(dims, to=to)
-        return da.copy()
-
-
 class Detector:
     def __init__(self, detector_number: sc.Variable):
         self._data = sc.DataArray(
@@ -295,7 +255,7 @@ class RollingDetectorView(Detector):
         else:
             indices = sc.ones(sizes=self.data.sizes, dtype='int32', unit=None)
             indices = sc.cumsum(indices, mode='exclusive')
-            if isinstance(self._projection, LogicalView):
+            if self._projection is not None:
                 indices = self._projection(indices)
         return roi.ROIFilter(indices=indices, norm=norm)
 
@@ -357,17 +317,23 @@ class RollingDetectorView(Detector):
         )
 
     @staticmethod
-    def from_detector_and_logical_view(
-        detector: CalibratedDetector[SampleRun],
-        window: RollingDetectorViewWindow,
-        projection: LogicalView,
-    ) -> RollingDetectorView:
-        """Helper for constructing via a Sciline workflow."""
-        return RollingDetectorView(
-            detector_number=detector.coords['detector_number'],
-            window=window,
-            projection=projection,
-        )
+    def from_detector_with_projection(
+        projection: Callable[[sc.DataArray], sc.DataArray] | None,
+    ) -> Callable[
+        [CalibratedDetector[SampleRun], RollingDetectorViewWindow], RollingDetectorView
+    ]:
+        def factory(
+            detector: CalibratedDetector[SampleRun],
+            window: RollingDetectorViewWindow,
+        ) -> RollingDetectorView:
+            """Helper for constructing via a Sciline workflow."""
+            return RollingDetectorView(
+                detector_number=detector.coords['detector_number'],
+                window=window,
+                projection=projection,
+            )
+
+        return factory
 
     @staticmethod
     def from_nexus(
@@ -375,7 +341,9 @@ class RollingDetectorView(Detector):
         *,
         detector_name: str,
         window: int,
-        projection: Literal['xy_plane', 'cylinder_mantle_z'] | LogicalView,
+        projection: Literal['xy_plane', 'cylinder_mantle_z']
+        | Callable[[sc.DataArray], sc.DataArray]
+        | None = None,
         resolution: dict[str, int] | None = None,
         pixel_noise: Literal['cylindrical'] | sc.Variable | None = None,
     ) -> RollingDetectorView:
@@ -396,10 +364,12 @@ class RollingDetectorView(Detector):
             Size of the rolling window.
         projection:
             Projection to use for the detector data. This can be a string selecting a
-            predefined projection or a LogicalView instance.
+            predefined projection or a function that takes a DataArray and returns a
+            DataArray. The predefined projections are 'xy_plane' and
+            'cylinder_mantle_z'.
         resolution:
-            Resolution to use for histogramming the detector data. Not required when the
-            projection is a LogicalView.
+            Resolution to use for histogramming the detector data. Only required for
+            'xy_plane' and 'cylinder_mantle_z' projections.
         pixel_noise:
             Noise to add to the pixel positions. This can be a scalar value to add
             Gaussian noise to the pixel positions or the string 'cylindrical' to add
@@ -413,13 +383,7 @@ class RollingDetectorView(Detector):
             noise_replica_count = 16
         wf = GenericNeXusWorkflow(run_types=[SampleRun], monitor_types=[])
         wf[RollingDetectorViewWindow] = window
-        if isinstance(projection, LogicalView):
-            wf[LogicalView] = projection
-            wf[NeXusTransformation[snx.NXdetector, SampleRun]] = NeXusTransformation[
-                snx.NXdetector, SampleRun
-            ](sc.scalar(1))
-            wf.insert(RollingDetectorView.from_detector_and_logical_view)
-        elif projection == 'cylinder_mantle_z':
+        if projection == 'cylinder_mantle_z':
             wf.insert(make_cylinder_mantle_coords)
             wf.insert(RollingDetectorView.from_detector_and_histogrammer)
             wf[DetectorViewResolution] = resolution
@@ -428,7 +392,12 @@ class RollingDetectorView(Detector):
             wf.insert(RollingDetectorView.from_detector_and_histogrammer)
             wf[DetectorViewResolution] = resolution
         else:
-            raise ValueError(f"Invalid {projection=}.")
+            wf[NeXusTransformation[snx.NXdetector, SampleRun]] = NeXusTransformation[
+                snx.NXdetector, SampleRun
+            ](sc.scalar(1))
+            wf.insert(
+                RollingDetectorView.from_detector_with_projection(projection=projection)
+            )
         if isinstance(pixel_noise, sc.Variable):
             wf.insert(gaussian_position_noise)
             wf[PositionNoiseSigma] = pixel_noise

{essreduce-25.2.4 → essreduce-25.2.6}/src/ess/reduce/nexus/_nexus_loader.py

@@ -180,10 +180,15 @@ def _attempt_to_open_without_locking(
     # HDF5 tracks file locking flags internally within a single process.
     # If the same file is opened multiple times, we can get a flag mismatch.
     # We can try opening without locking, maybe this matches the original flags.
-    if "file locking flag values don't match" in err.args[0]:
-        return True
-    if "file locking 'ignore disabled locks' flag values don't match" in err.args[0]:
-        return True
+    error_message = err.args[0]
+    if isinstance(error_message, str):
+        if "file locking flag values don't match" in error_message:
+            return True
+        if (
+            "file locking 'ignore disabled locks' flag values don't match"
+            in error_message
+        ):
+            return True
     return False
 
 

{essreduce-25.2.4 → essreduce-25.2.6}/src/ess/reduce/streaming.py

@@ -68,7 +68,18 @@ class Accumulator(ABC, Generic[T]):
     def _do_push(self, value: T) -> None: ...
 
     @property
-    @abstractmethod
+    def is_empty(self) -> bool:
+        """
+        Check if the accumulator is empty.
+
+        Returns
+        -------
+        :
+            True if the accumulator is empty, False otherwise.
+        """
+        return False
+
+    @property
     def value(self) -> T:
         """
         Get the accumulated value.
@@ -77,6 +88,24 @@ class Accumulator(ABC, Generic[T]):
         -------
         :
             Accumulated value.
+
+        Raises
+        ------
+        ValueError
+            If the accumulator is empty.
+        """
+        if self.is_empty:
+            raise ValueError("Cannot get value from empty accumulator")
+        return self._get_value()
+
+    @abstractmethod
+    def _get_value(self) -> T:
+        """Return the accumulated value, assuming it exists."""
+
+    @abstractmethod
+    def clear(self) -> None:
+        """
+        Clear the accumulator, resetting it to its initial state.
         """
 
 
@@ -92,7 +121,10 @@ class EternalAccumulator(Accumulator[T]):
         self._value: T | None = None
 
     @property
-    def value(self) -> T:
+    def is_empty(self) -> bool:
+        return self._value is None
+
+    def _get_value(self) -> T:
         return deepcopy(self._value)
 
     def _do_push(self, value: T) -> None:
@@ -101,6 +133,10 @@ class EternalAccumulator(Accumulator[T]):
         else:
             self._value += value
 
+    def clear(self) -> None:
+        """Clear the accumulated value."""
+        self._value = None
+
 
 class RollingAccumulator(Accumulator[T]):
     """
@@ -121,7 +157,10 @@ class RollingAccumulator(Accumulator[T]):
         self._values: list[T] = []
 
     @property
-    def value(self) -> T:
+    def is_empty(self) -> bool:
+        return len(self._values) == 0
+
+    def _get_value(self) -> T:
         # Naive and potentially slow implementation if values and/or window are large!
         return sc.reduce(self._values).sum()
 
@@ -130,6 +169,68 @@ class RollingAccumulator(Accumulator[T]):
         if len(self._values) > self._window:
             self._values.pop(0)
 
+    def clear(self) -> None:
+        """Clear the accumulated values."""
+        self._values = []
+
+
+class MinAccumulator(Accumulator):
+    """Keeps the minimum value seen so far.
+
+    Only supports scalar values.
+    """
+
+    def __init__(self, **kwargs: Any) -> None:
+        super().__init__(**kwargs)
+        self._cur_min: sc.Variable | None = None
+
+    def _do_push(self, value: sc.Variable) -> None:
+        if self._cur_min is None:
+            self._cur_min = value
+        else:
+            self._cur_min = min(self._cur_min, value)
+
+    @property
+    def is_empty(self) -> bool:
+        """Check if the accumulator has collected a minimum value."""
+        return self._cur_min is None
+
+    def _get_value(self) -> Any:
+        return self._cur_min
+
+    def clear(self) -> None:
+        """Clear the accumulated minimum value."""
+        self._cur_min = None
+
+
+class MaxAccumulator(Accumulator):
+    """Keeps the maximum value seen so far.
+
+    Only supports scalar values.
+    """
+
+    def __init__(self, **kwargs: Any) -> None:
+        super().__init__(**kwargs)
+        self._cur_max: sc.Variable | None = None
+
+    @property
+    def is_empty(self) -> bool:
+        """Check if the accumulator has collected a maximum value."""
+        return self._cur_max is None
+
+    def _do_push(self, value: sc.Variable) -> None:
+        if self._cur_max is None:
+            self._cur_max = value
+        else:
+            self._cur_max = max(self._cur_max, value)
+
+    def _get_value(self) -> sc.Variable | None:
+        return self._cur_max
+
+    def clear(self) -> None:
+        """Clear the accumulated maximum value."""
+        self._cur_max = None
+
 
 class StreamProcessor:
     """
@@ -299,6 +400,16 @@ class StreamProcessor:
             self._finalize_workflow[key] = self._accumulators[key].value
         return self._finalize_workflow.compute(self._target_keys)
 
+    def clear(self) -> None:
+        """
+        Clear all accumulators, resetting them to their initial state.
+
+        This is useful for restarting a streaming computation without
+        creating a new StreamProcessor instance.
+        """
+        for accumulator in self._accumulators.values():
+            accumulator.clear()
+
 
 def _find_descendants(
     workflow: sciline.Pipeline, keys: tuple[sciline.typing.Key, ...]

{essreduce-25.2.4 → essreduce-25.2.6/src/essreduce.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: essreduce
-Version: 25.2.4
+Version: 25.2.6
 Summary: Common data reduction tools for the ESS facility
 Author: Scipp contributors
 License: BSD 3-Clause License

{essreduce-25.2.4 → essreduce-25.2.6}/tests/live/raw_test.py

@@ -354,13 +354,12 @@ def test_ROIFilter_from_trivial_RollingDetectorView() -> None:
     assert sc.identical(scale, sc.ones(dims=['detector_number'], shape=[2]))
 
 
-def test_ROIFilter_from_RollingDetectorView_with_LogicalView() -> None:
-    logical_view = raw.LogicalView(select={'z': 0})
+def test_ROIFilter_from_RollingDetectorView_with_custom_projection() -> None:
     detector_number = sc.array(
         dims=['x', 'y', 'z'], values=[[[1, 2], [3, 4]], [[5, 6], [7, 8]]], unit=None
     )
     view = raw.RollingDetectorView(
-        detector_number=detector_number, window=2, projection=logical_view
+        detector_number=detector_number, window=2, projection=lambda da: da['z', 0]
     )
     roi_filter = view.make_roi_filter()
     data = detector_number.copy()

{essreduce-25.2.4 → essreduce-25.2.6}/tests/live/roi_test.py

@@ -3,7 +3,7 @@
 import pytest
 import scipp as sc
 
-from ess.reduce.live import raw, roi
+from ess.reduce.live import roi
 
 
 @pytest.fixture
@@ -146,13 +146,12 @@ def test_apply_selection_fails_with_out_of_bounds_index():
         roi.apply_selection(data, selection=selection)
 
 
-@pytest.fixture
-def logical_view():
-    return raw.LogicalView(fold={'x': 3, 'y': 4, 'z': 2}, select={'z': 0})
+def logical_view(da: sc.DataArray) -> sc.DataArray:
+    return da.fold(da.dim, sizes={'x': 3, 'y': 4, 'z': 2})['z', 0]
 
 
 @pytest.fixture
-def roi_filter(logical_view: raw.LogicalView):
+def roi_filter() -> roi.ROIFilter:
     indices = sc.ones(sizes={'detector_number': 24}, dtype='int32', unit=None)
     indices = sc.cumsum(indices, mode='exclusive')
     return roi.ROIFilter(logical_view(indices))

{essreduce-25.2.4 → essreduce-25.2.6}/tests/nexus/nexus_loader_test.py

@@ -1,6 +1,7 @@
 # SPDX-License-Identifier: BSD-3-Clause
 # Copyright (c) 2024 Scipp contributors (https://github.com/scipp)
-
+import os
+import sys
 from contextlib import contextmanager
 from io import BytesIO
 from pathlib import Path
@@ -661,17 +662,34 @@ def test_open_nexus_file_multiple_times(tmp_path: Path, locks: tuple[Any, Any])
             assert f1.name == f2.name
 
 
+def _in_conda_env():
+    return 'CONDA_PREFIX' in os.environ
+
+
+def _test_open_nexus_file_with_mismatched_locking(
+    tmp_path: Path, locks: tuple[Any, Any]
+) -> None:
+    from ess.reduce.nexus._nexus_loader import _open_nexus_file
+
+    path = FilePath(tmp_path / "file.nxs")
+    with snx.File(path, "w"):
+        pass
+
+    with _open_nexus_file(path, locking=locks[0]):
+        with pytest.raises(OSError, match="flag values don't match"):
+            _ = _open_nexus_file(path, locking=locks[1])
+
+
+@pytest.mark.skipif(
+    sys.platform in ("darwin", "win32")
+    or (sys.platform == "linux" and _in_conda_env()),
+    reason="HDF5 has different file locking flags on MacOS, Windows and Linux(conda)",
+)
 @pytest.mark.parametrize(
     "locks",
     [
-        (True, False),
         (True, None),
-        (False, True),
-        (False, None),
         (None, True),
-        (None, False),
-        # On a read-only filesystem, this would work:
-        (NoLockingIfNeeded, False),
         # This could be supported, but it could cause problems because the first
         # user expects the file to be locked.
         (True, NoLockingIfNeeded),
@@ -679,15 +697,31 @@ def test_open_nexus_file_multiple_times(tmp_path: Path, locks: tuple[Any, Any])
         (NoLockingIfNeeded, True),
     ],
 )
-def test_open_nexus_file_with_mismatched_locking(
+def test_open_nexus_file_with_mismatched_locking_pypi_linux(
     tmp_path: Path, locks: tuple[Any, Any]
 ) -> None:
-    from ess.reduce.nexus._nexus_loader import _open_nexus_file
+    _test_open_nexus_file_with_mismatched_locking(tmp_path, locks)
 
-    path = FilePath(tmp_path / "file.nxs")
-    with snx.File(path, "w"):
-        pass
 
-    with _open_nexus_file(path, locking=locks[0]):
-        with pytest.raises(OSError, match="flag values don't match"):
-            _ = _open_nexus_file(path, locking=locks[1])
+@pytest.mark.parametrize(
+    "locks",
+    [
+        (True, False),
+        (False, True),
+        (False, None),
+        (None, False),
+        # On a read-only filesystem, this would work:
+        (NoLockingIfNeeded, False),
+    ],
+)
+def test_open_nexus_file_with_mismatched_locking_all(
+    tmp_path: Path, locks: tuple[Any, Any]
+) -> None:
+    _test_open_nexus_file_with_mismatched_locking(tmp_path, locks)
+
+
+def test_open_nonexisting_file_raises_filenotfounderror():
+    from ess.reduce.nexus._nexus_loader import _open_nexus_file
+
+    with pytest.raises(FileNotFoundError):
+        _open_nexus_file(nexus.types.FilePath(Path("doesnotexist.hdf")))

{essreduce-25.2.4 → essreduce-25.2.6}/tests/streaming_test.py

@@ -45,6 +45,18 @@ def test_eternal_accumulator_does_not_modify_pushed_values() -> None:
     assert sc.identical(var, original)
 
 
+def test_eternal_accumulator_clear() -> None:
+    accum = streaming.EternalAccumulator()
+    var = sc.linspace(dim='x', start=0, stop=1, num=10)
+    for i in range(10):
+        accum.push(var[i].copy())
+    assert sc.identical(accum.value, sc.sum(var))
+    accum.clear()
+    assert accum.is_empty
+    with pytest.raises(ValueError, match="Cannot get value from empty accumulator"):
+        _ = accum.value
+
+
 def test_rolling_accumulator_sums_over_window() -> None:
     accum = streaming.RollingAccumulator(window=3)
     var = sc.linspace(dim='x', start=0, stop=1, num=10)
@@ -94,6 +106,76 @@ def test_rolling_accumulator_does_not_modify_pushed_values() -> None:
     assert sc.identical(var, original)
 
 
+def test_rolling_accumulator_clear() -> None:
+    accum = streaming.RollingAccumulator(window=3)
+    var = sc.linspace(dim='x', start=0, stop=1, num=10)
+    for i in range(5):
+        accum.push(var[i].copy())
+    assert sc.identical(accum.value, var[2:5].sum())
+    accum.clear()
+    assert accum.is_empty
+    with pytest.raises(ValueError, match="Cannot get value from empty accumulator"):
+        _ = accum.value
+
+
+def test_eternal_accumulator_is_empty() -> None:
+    accum = streaming.EternalAccumulator()
+    assert accum.is_empty
+    with pytest.raises(ValueError, match="Cannot get value from empty accumulator"):
+        _ = accum.value
+
+    var = sc.linspace(dim='x', start=0, stop=1, num=10)
+    accum.push(var[0].copy())
+    assert not accum.is_empty
+    assert sc.identical(accum.value, var[0])
+
+    accum.clear()
+    assert accum.is_empty
+    with pytest.raises(ValueError, match="Cannot get value from empty accumulator"):
+        _ = accum.value
+
+
+def test_rolling_accumulator_is_empty() -> None:
+    accum = streaming.RollingAccumulator(window=3)
+    assert accum.is_empty
+    with pytest.raises(ValueError, match="Cannot get value from empty accumulator"):
+        _ = accum.value
+
+    var = sc.linspace(dim='x', start=0, stop=1, num=10)
+    accum.push(var[0].copy())
+    assert not accum.is_empty
+    assert sc.identical(accum.value, var[0])
+
+    accum.clear()
+    assert accum.is_empty
+    with pytest.raises(ValueError, match="Cannot get value from empty accumulator"):
+        _ = accum.value
+
+
+def test_min_accumulator() -> None:
+    accum = streaming.MinAccumulator()
+    var = sc.array(dims=['x'], values=[1.0, 2.0, 3.0, 2.0, 1.0])
+    for scalar_var in var:
+        accum.push(scalar_var)
+    assert sc.identical(accum.value, sc.min(var))
+    accum.clear()
+    assert accum.is_empty
+    with pytest.raises(ValueError, match="Cannot get value from empty accumulator"):
+        _ = accum.value
+
+
+def test_max_accumulator() -> None:
+    accum = streaming.MaxAccumulator()
+    var = sc.array(dims=['x'], values=[1.0, 2.0, 3.0, 2.0, 1.0])
+    for scalar_var in var:
+        accum.push(scalar_var)
+    assert sc.identical(accum.value, sc.max(var))
+    accum.clear()
+    assert accum.is_empty
+    with pytest.raises(ValueError, match="Cannot get value from empty accumulator"):
+        _ = accum.value
+
+
 DynamicA = NewType('DynamicA', float)
 DynamicB = NewType('DynamicB', float)
 DynamicC = NewType('DynamicC', float)
@@ -161,9 +243,16 @@ def test_StreamProcessor_uses_custom_accumulator() -> None:
             pass
 
         @property
-        def value(self) -> sc.Variable:
+        def is_empty(self) -> bool:
+            return False
+
+        def _get_value(self) -> sc.Variable:
             return sc.scalar(42)
 
+        def clear(self) -> None:
+            # Nothing to clear
+            pass
+
     base_workflow = sciline.Pipeline(
         (make_static_a, make_accum_a, make_accum_b, make_target)
     )
@@ -370,10 +459,9 @@ def test_StreamProcessor_raises_given_partial_update_for_accumulator() -> None:
         accumulators=(Target, AccumC),  # Target depends on both A and B
     )
     # We can update either (A, B) and/or C...
-    result = streaming_wf.add_chunk({DynamicA: sc.scalar(1), DynamicB: sc.scalar(4)})
-    assert sc.identical(result[Target], sc.scalar(2 * 1.0 / 4.0))
-    assert result[AccumC] is None
-    result = streaming_wf.add_chunk({DynamicC: sc.scalar(11)})
+    streaming_wf.accumulate({DynamicA: sc.scalar(1), DynamicB: sc.scalar(4)})
+    streaming_wf.accumulate({DynamicC: sc.scalar(11)})
+    result = streaming_wf.finalize()
     assert sc.identical(result[Target], sc.scalar(2 * 1.0 / 4.0))
     assert sc.identical(result[AccumC], sc.scalar(11))
     result = streaming_wf.add_chunk({DynamicA: sc.scalar(2), DynamicB: sc.scalar(5)})
@@ -428,3 +516,35 @@ def test_StreamProcessor_raises_when_trying_to_update_non_dynamic_key() -> None:
         match=r'Got non-dynamic keys: {tests.streaming_test.Target}',
     ):
         result = streaming_wf.add_chunk({Target: sc.scalar(2)})
+
+
+def test_StreamProcessor_clear() -> None:
+    base_workflow = sciline.Pipeline(
+        (make_static_a, make_accum_a, make_accum_b, make_target)
+    )
+
+    # Reset call counter to ensure we can track it properly
+    make_static_a.call_count = 0
+
+    streaming_wf = streaming.StreamProcessor(
+        base_workflow=base_workflow,
+        dynamic_keys=(DynamicA, DynamicB),
+        target_keys=(Target,),
+        accumulators=(AccumA, AccumB),
+    )
+    # Add some data
+    result = streaming_wf.add_chunk({DynamicA: sc.scalar(1), DynamicB: sc.scalar(4)})
+    assert sc.identical(result[Target], sc.scalar(2 * 1.0 / 4.0))
+    result = streaming_wf.add_chunk({DynamicA: sc.scalar(2), DynamicB: sc.scalar(5)})
+    assert sc.identical(result[Target], sc.scalar(2 * 3.0 / 9.0))
+
+    # Make sure static_a was called exactly once
+    assert make_static_a.call_count == 1
+
+    # Clear and verify we get back to initial state
+    streaming_wf.clear()
+    result = streaming_wf.add_chunk({DynamicA: sc.scalar(1), DynamicB: sc.scalar(4)})
+    assert sc.identical(result[Target], sc.scalar(2 * 1.0 / 4.0))
+
+    # Static values should be preserved after clear, so call_count remains 1
+    assert make_static_a.call_count == 1