🔥 remove function: propagate_statistical_error & derivative

arpes/plotting/spin.py

@@ -16,7 +16,6 @@
 from arpes.analysis.statistics import mean_and_deviation
 from arpes.bootstrap import bootstrap
 from arpes.provenance import save_plot_provenance
-from arpes.utilities.math import polarization, propagate_statistical_error
 
 from .tof import scatter_with_std
 from .utils import label_for_dim, path_for_plot, polarization_colorbar, savefig
@@ -35,9 +34,6 @@
 )
 
 
-test_polarization = propagate_statistical_error(polarization)
-
-
 @save_plot_provenance
 def spin_colored_spectrum(
     spin_dr: xr.Dataset,

arpes/utilities/math.py

@@ -2,7 +2,7 @@
 from __future__ import annotations
 
 import itertools
-from collections.abc import Callable, Iterable
+from collections.abc import Iterable
 from typing import TYPE_CHECKING
 
 import numpy as np
@@ -16,44 +16,11 @@
     from numpy.typing import NDArray
 
 
-def derivative(f: Callable[..., float], arg_idx: int = 0) -> float:
-    """Defines a simple midpoint derivative."""
-
-    def d(*args: Incomplete):
-        args = list(args)
-        ref_arg = args[arg_idx]
-        d = ref_arg / 100
-        args[arg_idx] = ref_arg + d
-        high = f(*args)
-        args[arg_idx] = ref_arg - d
-        low = f(*args)
-        return (high - low) / (2 * d)
-
-    return d
-
-
 def polarization(up: NDArray[np.float_], down: NDArray[np.float_]) -> NDArray[np.float_]:
     """The equivalent normalized difference for a two component signal."""
     return (up - down) / (up + down)
 
 
-def propagate_statistical_error(f):
-    """To compute a function which propagates statistical error.
-
-    It Uses numerical derivatives and sampling.
-    """
-
-    def compute_propagated_error(*args):
-        running_sum = 0
-        for i, arg in enumerate(args):
-            df_darg_i = derivative(f, i)
-            running_sum += df_darg_i(*args) ** 2 * arg
-
-        return np.sqrt(running_sum)
-
-    return compute_propagated_error
-
-
 def shift_by(
     arr: NDArray[np.float_],
     value: xr.DataArray | NDArray[np.float_],