ForwardDiff and Zygote cannot automatically differentiate (AD) function from C^n to R that uses FFT

Specific Domains Numerics
23

did you ever reach an acceptable solution for this?