FFTW planned 2D real ldiv! unexpectedly overwrites input
Open
- Dominant language
- Julia
- Stars
- 300
- Forks
- 61
- PR merge metrics
- No merged PRs in 30d
Description
The following seems unexpected:
```julia
using FFTW, LinearAlgebra
y = rand(4,4)
p = plan_rfft(y)
x = p * y
x_orig = copy(x)
ldiv!(y,p,x)
x == x_orig # unexpectedly returns false
```
This only happens with FFTW and not MKL, only with 2D and real FFTs, and only with `ldiv!` not `mul!`.
Contributor guide
No contributing guide indexed for this repository
Assessment
This issue has not been assessed yet.