Project.toml

@@ -1,7 +1,7 @@
 name = "STFT"
 uuid = "58bb99bf-048b-48b7-93e7-1cbf3ee61509"
 authors = ["Szymon M. Woźniak <s@zymon.org>"]
-version = "1.3.0"
+version = "1.1.0"
 
 [deps]
 FFTW = "7a1cc6ca-52ef-59f5-83cd-3a7055c09341"

src/STFT.jl

@@ -5,21 +5,9 @@ using FFTW
 
 export stft, istft
 
-_fft(x::AbstractArray{<:Real}, d) = rfft(x, d)
-_fft(x::AbstractArray{<:Complex}, d) = fft(x, d)
+_fft(x::AbstractMatrix{<:Real}, d) = rfft(x, d)
+_fft(x::AbstractMatrix{<:Complex}, d) = fft(x, d)
 
-function malloc_stft(
-    x::M,
-    w::V,
-    L::I = zero(I),
-    N::I = length(w);
-) where {T<:Number, I<:Integer, V<:AbstractVector{T}, M<:AbstractMatrix{T}}
-    X, K = size(x)      # Length of the signal in samples
-    W = length(w)       # Length of the window in samples
-    S = (X-L) ÷ (W - L) # Number of segments
-    N = N < W ? W : N   # DFT size
-    zeros(T, (N, S, K)) # Allocate container for signal segments
-end
 
 
 doc_analysis = """
@@ -107,31 +95,24 @@ function analysis() end
 stft(x, w, L=0, N=length(w)) = analysis(x, w, L, N)
 
 function analysis(
-    x::M,
-    w::V,
-    L::I = zero(I),
-    N::I = length(w);
-) where {T<:Number, I<:Integer, V<:AbstractVector{T}, M<:AbstractMatrix{T}}
-    sc = malloc_stft(x, w, L, N)
-    N, S, K = sc |> size
-    W = w |> length
+    x::AbstractVector{T},
+    w::AbstractVector{T},
+    L::Integer = 0,
+    N::Integer = length(w);
+)::AbstractMatrix{<:Complex} where {T<:Number}
+    X = length(x)       # Length of the signal in samples
+    W = length(w)       # Length of the window in samples
+    H = W - L           # Hop
+    S = (X-L) ÷ H       # Number of segments
+    N = N < W ? W : N   # DFT size
+    sc = zeros(T, N, S) # Allocate container for signal segments
 
-    @turbo for k ∈ 1:K, s ∈ 1:S, n ∈ 1:W
-        sc[n, s, k] = w[n] * x[(s-1)*(W-L)+n, k]
+    @turbo for s ∈ 1:S, n ∈ 1:W # Slice the signal
+        sc[n, s] = w[n] * x[(s-1)*H+n]
     end
     _fft(sc, 1) # Convert segments to frequency-domain
 end
 
-function analysis(
-    x::V,
-    w::V,
-    L::I = zero(I),
-    N::I = length(w);
-)::Matrix{T |> complex} where {T<:Number, I<:Integer, V<:AbstractVector{T}}
-    xx = @view x[:,:]
-    @view analysis(xx, w, L, N)[:, :, 1]
-end
-
 function analysis(
     xs::AbstractArray{<:AbstractVector{T}},
     w::AbstractVector{T},
@@ -249,32 +230,14 @@ function synthesis(
 
     xs = irfft(X, N, 1) # Convert segments to time-domain
 
-    @turbo for s ∈ 0:(S-1)
-        ss = s*H # Segment shift
+    @turbo for s ∈ 1:S
+        ss = (s-1)*H # Segment start
         for n = 1:W
-            xn[ss+n] += xs[n, s+1] * w[n]
+            xn[ss+n] += xs[n, s] * w[n]
             xd[ss+n] += w²[n]
         end
     end
     xn ./ xd # Normalize
 end
 
-"""
-Real-value signal STFT with constant input size.
-"""
-function rSTFTm(A, w, L, N=length(w))
-    mem = STFT.malloc_stft(A, w, L, N)
-    N, S, K = mem |> size
-    W = w |> length
-    P = plan_rfft(mem, 1)
-    function f(x)
-        @turbo for k ∈ 1:K, s ∈ 1:S, n ∈ 1:W
-            mem[n, s, k] = w[n] * x[(s-1)*(W-L)+n, k]
-        end
-        P * mem
-    end
-    return f
-end
-
-
 end # module