regdiff 3.6.1__tar.gz → 3.6.2__tar.gz

{regdiff-3.6.1 → regdiff-3.6.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.3
 Name: regdiff
-Version: 3.6.1
+Version: 3.6.2
 Summary: Regularized differentiation of time-series as numpy ndarrays
 Requires-Dist: numpy>=2.5.4
 Requires-Dist: scipy>=1.18.1

{regdiff-3.6.1 → regdiff-3.6.2}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "regdiff"
-version = "3.6.1"
+version = "3.6.2"
 description = "Regularized differentiation of time-series as numpy ndarrays"
 requires-python = ">=3.14"
 dependencies = [

{regdiff-3.6.1 → regdiff-3.6.2}/pyproject.toml.orig

@@ -1,6 +1,6 @@
 [project]
 name = "regdiff"
-version = "3.6.1"
+version = "3.6.2"
 description = "Regularized differentiation of time-series as numpy ndarrays"
 requires-python = ">=3.14"
 dependencies = [

regdiff-3.6.2/src/regdiff/__init__.py

@@ -0,0 +1 @@
+from .regdiff import matrix_anti_diff, matrix_diff, reg_diff  # noqa: F401

{regdiff-3.6.1 → regdiff-3.6.2}/src/regdiff/regdiff.py

@@ -4,16 +4,6 @@ import scipy as sp
 
 # NOTE - all matrix/vector operations are called using np.matvel/np.subtract, etc. Only scalars should multiply ndarrays.
 
-
-def finite_diff(x):
-    """take the finite difference of the input vector. Performs no validation.
-    :param np.ndarray(x): data to take the finite difference of -- 1-dimensional array of length n
-
-    :return: diff(x) -- 1 dimensional array of length n-1
-    """
-    return np.subtract(x[1:], x[:-1])
-
-
 def matrix_diff(n, dx=1.0):
     """return the matrix representations of the first-order central finite first difference, first-order central second
     finite difference, and simple finite difference for a specified size and delta x increment.
@@ -149,11 +139,11 @@ def reg_diff(
     if diag:
         u_cost = [
             npl.norm(np.subtract(sp.integrate.cumulative_trapezoid(u), data_comp))
-            + alpha * npl.norm(finite_diff(u), ord=1)
+            + alpha * npl.norm(np.diff(u), ord=1)
         ]
     prev_cost = None
     for iter in range(max_iter):
-        c = np.reciprocal(np.sqrt(np.power(finite_diff(u), 2) + ep))
+        c = np.reciprocal(np.sqrt(np.power(np.diff(u), 2) + ep))
         En = np.diag(c)
         Ln = np.matmul(np.matmul(d_tp, En), d)
         gn = np.add(np.subtract(np.matvec(a_tp_a, u), a_tp_b), alpha * np.matvec(Ln, u))
@@ -162,7 +152,7 @@ def reg_diff(
         u = np.subtract(u, s)
         cost_n = npl.norm(
             np.subtract(sp.integrate.cumulative_trapezoid(u), data_comp)
-        ) + alpha * npl.norm(finite_diff(u), ord=1)
+        ) + alpha * npl.norm(np.diff(u), ord=1)
         if diag:
             print(
                 f"u, iteration {iter + 1}: cost: {cost_n}, relative change: {npl.norm(s) / npl.norm(u)}"
@@ -196,11 +186,11 @@ def reg_diff(
                         data_comp,
                     )
                 )
-                + alpha2 * npl.norm(finite_diff(v), ord=1)
+                + alpha2 * npl.norm(np.diff(v), ord=1)
             ]
         prev_cost = None
         for iter in range(max_iter):
-            c = np.reciprocal(np.sqrt(np.power(finite_diff(v), 2) + ep))
+            c = np.reciprocal(np.sqrt(np.power(np.diff(v), 2) + ep))
             En = np.diag(c)
             Ln = np.matmul(np.matmul(d_tp, En), d)
             gn = np.add(
@@ -216,7 +206,7 @@ def reg_diff(
                     ),
                     data_comp,
                 )
-            ) + alpha2 * npl.norm(finite_diff(v), ord=1)
+            ) + alpha2 * npl.norm(np.diff(v), ord=1)
             if diag:
                 print(
                     f"v, iteration {iter + 1}: cost: {cost_n}, relative change: {npl.norm(s) / npl.norm(v)}"

regdiff-3.6.1/src/regdiff/__init__.py

@@ -1 +0,0 @@
-from .regdiff import finite_diff, matrix_anti_diff, matrix_diff, reg_diff  # noqa: F401