Turbo doesn't work?

Author: dannys4

src/NSFFT.jl

@@ -1,7 +1,7 @@
 module NSFFT
 
 using Primes, StaticArrays, LoopVectorization
-export pow2FFT!, myfft_sv
+export pow2FFT!, DFT!
 
 include("algos.jl")
 

src/algos.jl

@@ -48,4 +48,29 @@ function pow2FFT!(out::AbstractVector{T}, in::AbstractVector{T}, ::Val{FFT_BACKW
         out[j+m] = out_j - wj*out[j+m]
         wj *= w1
     end
+end
+
+function DFT!(out::AbstractVector{T}, in::AbstractVector{T}, ::Val{FFT_FORWARD}) where {T<:Complex}
+    N = length(out)
+    inc = 2*π/N
+    wn² = wn = w = T(cos(inc), -sin(inc));
+    wn_1 = T(1., 0.);
+
+    tmp = in[1];
+    out .= tmp;
+    tmp = sum(in)
+    out[1] = tmp;
+
+    wk = wn²;
+    @turbo for d in 1:N
+        for k in (d+1):N
+            wk *= wn
+            out[d] = in[k]*wk + out[d]
+            out[k] = in[d]*wk + out[k]
+        end
+        wn_1 = wn
+        wn *= w
+        wn² *= (wn*wn_1)
+        wk = wn²
+    end
 end