perf(dsp): use conv for simple filt as it is much faster on complex s…

…ignals
org-arl · Jul 24, 2024 · a9a968e · a9a968e
1 parent b5ac0f6
commit a9a968e
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/src/dsp.jl b/src/dsp.jl
@@ -410,7 +410,11 @@ function pll(x::AbstractVecOrMat, fc=0.0, bandwidth=1e-3; fs=framerate(x))
   signal(y, fs)
 end
 
-sfilt(f::AbstractVector{<:Number}, x::AbstractVector) = signal(filt(f, samples(x)), framerate(x))
+function sfilt(f::AbstractVector{<:Number}, x::AbstractVector)
+  x̄ = samples(x)
+  signal(conv(f, x̄)[eachindex(x̄)], framerate(x))
+end
+
 sfilt(f::AbstractVector{<:Number}, x::AbstractVector, si) = signal(filt(f, samples(x), si), framerate(x))
 sfilt(b::AbstractVector{<:Number}, a::Union{Number,AbstractVector}, x::AbstractVector) = signal(filt(b, a, samples(x)), framerate(x))
 sfilt(b::AbstractVector{<:Number}, a::Union{Number,AbstractVector}, x::AbstractVector, si) = signal(filt(b, a, samples(x), si), framerate(x))