Skip to content

Attempt to fix FFT #95

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Draft
jagoosw wants to merge 2 commits into
base: master
Choose a base branch
from
Draft

Attempt to fix FFT #95

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
bfe3ab4
attempted to fix, doesn't run (seg faults)
jagoosw Sep 25, 2023
c6c4b3c
Merge branch 'master' into jsw/fft
jagoosw Nov 27, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
74 changes: 35 additions & 39 deletions src/DSP.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -372,42 +372,38 @@ function plan_destroy(setup::DFTSetup)
end


# const FFT_FORWARD = +1
# const FFT_INVERSE = -1
# const SIGNAL_STRIDE = 1
#
# typealias vDSP_Stride Clong
# typealias vDSP_Length Culong
# typealias FFTDirection Cint
#
# immutable DSPDoubleSplitComplex
# realp::Ptr{Float64}
# imagp::Ptr{Float64}
# end
#
# DSPDoubleSplitComplex(realp::Vector{Float64},imagp::Vector{Float64})=DSPDoubleSplitComplex(pointer(realp),pointer(imagp))
#
#
# function plan_fft(n,radix::Integer=2)
# @assert isinteger(log2(n))
# logn=round(Int,log2(n))
# ccall(("vDSP_create_fftsetupD",libacc),Ptr{Cvoid},(Cuint,Cint),logn,radix)
# end
#
# function newfft(r::Vector{Complex128},plan)
# n=length(r)
# @assert isinteger(log2(n))
# logn=round(Int,log2(n))
#
# realp,imagp=real(r),imag(r) # keep references to avoid garbage collection
# vals=DSPDoubleSplitComplex(realp,imagp)
# retr,reti=Array{Float64}(undef,n),Array{Float64}(undef,n) # keep references to avoid garbage collection
# ret=DSPDoubleSplitComplex(retr,reti)
#
#
# ccall(("vDSP_fft_zopD",libacc),Cvoid,
# (Ptr{Cvoid},DSPDoubleSplitComplex,vDSP_Stride, DSPDoubleSplitComplex,vDSP_Stride, vDSP_Length,FFTDirection),
# plan, vals, SIGNAL_STRIDE, ret, SIGNAL_STRIDE,logn, FFT_FORWARD)
#
# retr+im*reti
# End
const FFT_FORWARD = 1
const FFT_INVERSE = -1
const SIGNAL_STRIDE = 1

struct DSPDoubleSplitComplex
realp::Ptr{Float64}
imagp::Ptr{Float64}
end

DSPDoubleSplitComplex(realp::Vector{Float64}, imagp::Vector{Float64}) = DSPDoubleSplitComplex(pointer(realp), pointer(imagp))

function plan_fft(n, radix::Integer = 2)
logn = log2(n)
@assert isinteger(logn)
logn = round(Int, logn)
ccall(("vDSP_create_fftsetupD", libacc), Ptr{Cvoid}, (Cuint, Cint), logn, radix)
end

function newfft(r, plan)
n = length(r)
@assert isinteger(log2(n))
logn = round(Int, log2(n))

realp, imagp = real(r), imag(r) # keep references to avoid garbage collection
vals = DSPDoubleSplitComplex(realp, imagp)
retr, reti = Array{Float64}(undef, n), Array{Float64}(undef, n) # keep references to avoid garbage collection
ret = DSPDoubleSplitComplex(retr, reti)


ccall(("vDSP_fft_zopD",libacc), Cvoid,
(Ptr{Cvoid}, DSPDoubleSplitComplex, Cint, DSPDoubleSplitComplex, Cint, Cint, Cint),
plan, vals, SIGNAL_STRIDE, ret, SIGNAL_STRIDE, logn, FFT_FORWARD)

return retr + im*reti
end