JuliaMath / JuliaMath/FFTW.jl

`plan_rfft` eltype should be `Complex`

Open
#256 0 comments 0 reactions 0 assignees View on GitHub
Dominant language
Julia
Stars
300
Forks
61
PR merge metrics
No merged PRs in 30d

Description

This caught me by surprise:
```julia
julia> x = randn(5); P = plan_rfft(x);

julia> eltype(P)
Float64

julia> P * x
3-element Vector{ComplexF64}:
-0.3514539021978911 + 0.0im
-0.8023047651675712 + 1.243731722698685im
2.8759555200391196 - 2.679171588261927im
```
If a "plan" is in some sense analogous to a matrix the `eltype` should be type of the output, in this case `ComplexF64`.

Contributor guide

No contributing guide indexed for this repository

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.