From 5ea0991c57da443bdffa58039528b03895992ca3 Mon Sep 17 00:00:00 2001
From: Kamal Choudhary <writetokamal.1989@gmail.com>
Date: Mon, 23 Sep 2024 11:12:00 -0400
Subject: [PATCH] Update utils.py

---
 jarvis/core/utils.py | 31 ++++++++++++++++++++++++++++++-
 1 file changed, 30 insertions(+), 1 deletion(-)

diff --git a/jarvis/core/utils.py b/jarvis/core/utils.py
index 7fde5013..69e0298b 100644
--- a/jarvis/core/utils.py
+++ b/jarvis/core/utils.py
@@ -1,7 +1,8 @@
 """Set of useful utility functions."""
 
 from collections import OrderedDict
-from collections import defaultdict
+from scipy import sparse
+from scipy.sparse.linalg import spsolve
 import random
 import numpy as np
 import math
@@ -334,6 +335,34 @@ def cos_formula(a, b, c):
     return np.arccos(res)
 
 
+def baseline_als(y, lam, p, niter=10):
+    """
+    Adaptive Least Squares fitting for baseline correction
+    
+    Parameters:
+    y: array_like
+        Input signal
+    lam: float
+        Lambda (smoothness)
+    p: float
+        Asymmetry
+    niter: int, optional
+        Number of iterations
+    
+    Returns:
+    array_like
+        The estimated baseline
+    """
+    L = len(y)
+    D = sparse.diags([1,-2,1],[0,-1,-2], shape=(L,L-2))
+    w = np.ones(L)
+    for i in range(niter):
+        W = sparse.spdiags(w, 0, L, L)
+        Z = W + lam * D.dot(D.transpose())
+        z = spsolve(Z, w*y)
+        w = p * (y > z) + (1-p) * (y < z)
+    return z
+
 # def is_xml_valid(xsd="jarvisdft.xsd", xml="JVASP-1002.xml"):
 #   """Check if XML is valid."""
 #   xml_file = etree.parse(xml)