Merge pull request #73 from butala/master

Use scipy.signal.correlate for efficiency
cokelaer · Dec 8, 2024 · 822b3df · 822b3df
2 parents 9afb2d7 + 0c81dda
commit 822b3df
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/spectrum/correlation.py b/src/spectrum/correlation.py
@@ -21,6 +21,7 @@
 
 """
 import numpy as np
+import scipy.signal
 
 
 __all__ = ['CORRELATION', 'xcorr']
@@ -208,7 +209,7 @@ def xcorr(x, y=None, maxlags=None, norm='biased'):
         assert maxlags <= N, 'maxlags must be less than data length'
         lags = np.arange(N-maxlags-1, N+maxlags)
 
-    res = np.correlate(x, y, mode='full')
+    res = scipy.signal.correlate(x, y, mode='full')
 
     if norm == 'biased':
         Nf = float(N)