Working Reals

Author: dannys4

src/algos.jl

@@ -127,8 +127,8 @@ function fft_dft!(out::AbstractVector{Complex{T}}, in::AbstractVector{T}, N::Int
         wk *= w
         wkn = wk
     end
-    @inbounds @turbo for i in 1:halfN-1
-        out[start_out + stride_out*(N-i)] = conj(out[start_out + stride_out*i])
+    @inbounds @turbo for k in halfN+1:N-1
+        out[start_out + stride_out*k] = conj(out[start_out + stride_out*(N-k)])
     end
 end
 

test/real_forward.jl

@@ -1,6 +1,6 @@
 using FFTA, Test
 test_nums = [8, 11, 15, 100]
-@testset verbose = true "forward" begin 
+@testset verbose = true " forward" begin 
     for N in test_nums
         x = ones(Float64, N)
         y = fft(x)

test/runtests.jl

@@ -12,9 +12,17 @@ end
     @testset verbose = true "Complex" begin
         include("complex_forward.jl")
         include("complex_backward.jl")
+        x = rand(ComplexF64, 100)
+        y = fft(x)
+        x2 = bfft(y)/length(x)
+        @test x ≈ x2 atol=1e-12
     end
     @testset verbose = true "Real" begin
         include("real_forward.jl")
         include("real_backward.jl")
+        x = rand(Float64, 100)
+        y = fft(x)
+        x2 = bfft(y)/length(x)
+        @test x ≈ x2 atol=1e-12
     end
 end