diff --git a/src/dspbase.jl b/src/dspbase.jl index 4c985925..e7db6d3f 100644 --- a/src/dspbase.jl +++ b/src/dspbase.jl @@ -797,7 +797,7 @@ function conv(u::AbstractVector{T}, v::Transpose{T,<:AbstractVector}, A::Abstrac m = length(u) + size(A, 1) - 1 n = length(v) + size(A, 2) - 1 B = zeros(T, m, n) - B[1:size(A, 1), 1:size(A, 2)] = A + B[CartesianIndices(A)] = A u, v = fft.(_zeropad.((u, transpose(v)), (m, n))) p = plan_fft(B) C = inv(p) * ((p * B) .*= u .* transpose(v))