@@ -840,7 +840,9 @@ function _imfilter_fft!(r::AbstractCPU{FFT},
840
840
out
841
841
end
842
842
843
- function planned_fft (A:: AbstractArray{T,N} , kernel:: Tuple{AbstractArray,Vararg{AbstractArray}} , border:: BorderSpecAny = Pad (:replicate )) where {T,N}
843
+ function planned_fft (A:: AbstractArray{T,N} ,
844
+ kernel:: Tuple{AbstractArray,Vararg{AbstractArray}} ,
845
+ border:: BorderSpecAny = Pad (:replicate )) where {T,N}
844
846
bord = border (kernel, A, Algorithm. FFT ())
845
847
_A = padarray (T, A, bord)
846
848
p1 = plan_rfft (_A)
@@ -853,12 +855,6 @@ function planned_fft(A::AbstractArray{T,N}, kernel::Tuple{AbstractArray,Vararg{A
853
855
return Algorithm. FFT (p1, p2, p3)
854
856
end
855
857
856
- function kernel_plan_rfft (A:: AbstractArray{T,N} , kernel:: Tuple{AbstractArray,Vararg{AbstractArray}} ) where {T,N}
857
- kern = samedims (A, kernelconv (kernel... ))
858
- krn = FFTView (zeros (eltype (kern), map (length, axes (A))))
859
- return plan_rfft (krn), krn
860
- end
861
-
862
858
function filtfft (A, krn, plan_A:: FFTW.rFFTWPlan , plan_krn:: FFTW.rFFTWPlan , plan_B:: FFTW.AbstractFFTs.ScaledPlan )
863
859
B = plan_A * A
864
860
B .*= conj! (plan_krn * krn)
0 commit comments