Skip to content

Commit 59b6ea8

Browse files
committed
Fix radix 3
1 parent 3a7d2bb commit 59b6ea8

File tree

1 file changed

+3
-3
lines changed

1 file changed

+3
-3
lines changed

src/algos.jl

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -188,14 +188,14 @@ function fft_pow3!(out::AbstractVector{T}, in::AbstractVector{U}, N::Int, start_
188188
end
189189

190190
# Size of subproblem
191-
Nprime = N/3
191+
Nprime = N ÷ 3
192192

193193
ds = direction_sign(d)
194194

195195
# Dividing into subproblems
196196
fft_pow3!(out, in, Nprime, start_out, stride_out, start_in, stride_in*3, d, plus120, minus120)
197-
fft_pow3!(out, in, Nprime, start_out + N_prime*stride_out, stride_out, start_in + stride_in, stride_in*3, d, plus120, minus120)
198-
fft_pow3!(out, in, Nprime, start_out + 2*N_prime*stride_out, stride_out, start_in + 2*stride_in, stride_in*3, d, plus120, minus120)
197+
fft_pow3!(out, in, Nprime, start_out + Nprime*stride_out, stride_out, start_in + stride_in, stride_in*3, d, plus120, minus120)
198+
fft_pow3!(out, in, Nprime, start_out + 2*Nprime*stride_out, stride_out, start_in + 2*stride_in, stride_in*3, d, plus120, minus120)
199199

200200
w1 = convert(T, cispi(ds*2/N))
201201
w2 = convert(T, cispi(ds*4/N))

0 commit comments

Comments
 (0)