From 8b729604ea1419ef086e5573637663deb4e33df4 Mon Sep 17 00:00:00 2001
From: Blaise deB Frederick <blaise.frederick@gmail.com>
Date: Tue, 9 Apr 2024 13:50:34 -0400
Subject: [PATCH] Proactive fixes of some things that will break in numpy 2.0.0

---
 capcalc/fit.py | 16 ++++++++--------
 1 file changed, 8 insertions(+), 8 deletions(-)

diff --git a/capcalc/fit.py b/capcalc/fit.py
index fc92579..5307f35 100644
--- a/capcalc/fit.py
+++ b/capcalc/fit.py
@@ -1555,7 +1555,7 @@ def mlregress(x, y, intercept=True):
         xc = x
         beta = np.ones(p)
 
-    solution = np.linalg.lstsq(np.mat(xc).T, np.mat(y).T, rcond=-1)
+    solution = np.linalg.lstsq(np.asmatrix(xc).T, np.asmatrix(y).T, rcond=-1)
 
     # Computation of the coefficient of determination.
     Rx = np.atleast_2d(np.corrcoef(x, rowvar=1))
@@ -1737,7 +1737,7 @@ def peakdetect(y_axis, x_axis=None, lookahead=200, delta=0.0):
 
     # maxima and minima candidates are temporarily stored in
     # mx and mn respectively
-    mn, mx = np.Inf, -np.Inf
+    mn, mx = np.inf, -np.inf
 
     # Only detect peak if there is 'lookahead' amount of points after it
     for index, (x, y) in enumerate(zip(x_axis[:-lookahead], y_axis[:-lookahead])):
@@ -1749,15 +1749,15 @@ def peakdetect(y_axis, x_axis=None, lookahead=200, delta=0.0):
             mnpos = x
 
         ####look for max####
-        if y < mx - delta and mx != np.Inf:
+        if y < mx - delta and mx != np.inf:
             # Maxima peak candidate found
             # look ahead in signal to ensure that this is a peak and not jitter
             if y_axis[index : index + lookahead].max() < mx:
                 max_peaks.append([mxpos, mx])
                 dump.append(True)
                 # set algorithm to only find minima now
-                mx = np.Inf
-                mn = np.Inf
+                mx = np.inf
+                mn = np.inf
                 if index + lookahead >= length:
                     # end is within lookahead no more peaks can be found
                     break
@@ -1767,15 +1767,15 @@ def peakdetect(y_axis, x_axis=None, lookahead=200, delta=0.0):
                 #    mxpos = x_axis[np.where(y_axis[index:index+lookahead]==mx)]
 
         ####look for min####
-        if y > mn + delta and mn != -np.Inf:
+        if y > mn + delta and mn != -np.inf:
             # Minima peak candidate found
             # look ahead in signal to ensure that this is a peak and not jitter
             if y_axis[index : index + lookahead].min() > mn:
                 min_peaks.append([mnpos, mn])
                 dump.append(False)
                 # set algorithm to only find maxima now
-                mn = -np.Inf
-                mx = -np.Inf
+                mn = -np.inf
+                mx = -np.inf
                 if index + lookahead >= length:
                     # end is within lookahead no more peaks can be found
                     break