papyri
This is a pre-release version (2.0.0.dev0+git20260907.41eeb59).Go to latest (1.18.0)
{ } Raw JSON

bundles / scipy 2.0.0.dev0+git20260907.41eeb59 / scipy / signal / _signaltools / oaconvolve

function

scipy.signal._signaltools:oaconvolve

source: build-install/usr/lib/python3.14/site-packages/scipy/signal/_signaltools.py:851

Signature

def oaconvolve(in1, in2, mode = full, axes = None)

Summary

Convolve two N-dimensional arrays using the overlap-add method.

Extended Summary

Convolve in1 and in2 using the overlap-add method, with the output size determined by the mode argument.

This is generally much faster than convolve for large arrays (n > ~500), and generally much faster than fftconvolve when one array is much larger than the other, but can be slower when only a few output values are needed or when the arrays are very similar in shape, and can only output float arrays (int or object array inputs will be cast to float).

Parameters

in1 : array_like

First input.

in2 : array_like

Second input. Should have the same number of dimensions as in1.

mode : str {'full', 'valid', 'same'}, optional

A string indicating the size of the output:

full

The output is the full discrete linear convolution of the inputs. (Default)

valid

The output consists only of those elements that do not rely on the zero-padding. In 'valid' mode, either in1 or in2 must be at least as large as the other in every dimension.

same

Output size will be N, the same size as in1, centered with respect to the 'full' output. Boundary effects may be visible.

axes : int or array_like of ints or None, optional

Axes over which to compute the convolution. The default is over all axes.

Returns

out : array

An N-dimensional array containing a subset of the discrete linear convolution of in1 with in2.

The output size along each axis depends on mode. For input sizes N and M along a given axis:

  • fullN + M - 1

  • sameN (same as the size of in1)

  • validmax(N, M) - min(N, M) + 1

Notes

Examples

Convolve a 100,000 sample signal with a 512-sample filter. 
import numpy as np
from scipy import signal
rng = np.random.default_rng()
sig = rng.standard_normal(100000)
filt = signal.firwin(512, 0.01)
fsig = signal.oaconvolve(sig, filt)
✓
import matplotlib.pyplot as plt
fig, (ax_orig, ax_mag) = plt.subplots(2, 1)
✓
ax_orig.plot(sig)
ax_orig.set_title('White noise')
ax_mag.plot(fsig)
ax_mag.set_title('Filtered noise')
✗
Text(0.5, 1.0, 'Filtered noise')

fig.tight_layout()
fig.show()
✓

See also

convolve

Uses the direct convolution or FFT convolution algorithm depending on which is faster.

fftconvolve

An implementation of convolution using FFT.

Aliases

  • scipy.signal.oaconvolve

Referenced by

This package