Can't use empty vectors in `th_cheb2leg` or `th_leg2cheb`
まだ誰も着手していません。
評価
調査の方向性
src/toeplitzhankel.jl から開始し、特にスタックトレースで特定された plan_th_leg2cheb! と plan_th_cheb2leg! のパスを確認してください。Float64[] を使って両方の呼び出しを再現し、長さ 0 の次元がどのように処理されるかを追跡してください。th_leg2cheb(Float64[]) と th_cheb2leg(Float64[]) が、報告された BoundsError や無効な次元エラーなしに完了すれば完了です。
索引モデルが issue の本文から書いたものです。
説明
This originally came up in https://github.com/JuliaApproximation/ClassicalOrthogonalPolynomials.jl/issues/196
julia> FastTransforms.th_leg2cheb(Float64[])
ERROR: BoundsError: attempt to access 0×0 Matrix{Float64} at index [1, 1:0]
Stacktrace:
[1] throw_boundserror(A::Matrix{Float64}, I::Tuple{Int64, Base.Slice{Base.OneTo{Int64}}})
@ Base .\abstractarray.jl:737
[2] checkbounds
@ .\abstractarray.jl:702 [inlined]
[3] view
@ .\subarray.jl:184 [inlined]
[4] maybeview
@ .\views.jl:148 [inlined]
[5] dotview
@ .\broadcast.jl:1244 [inlined]
[6] _leg2chebTH_TLC(::Type{Float64}, mn::Tuple{Int64}, d::Int64)
@ FastTransforms C:\Users\danjv\.julia\packages\FastTransforms\S7huq\src\toeplitzhankel.jl:181
[7] _broadcast_getindex_evalf
@ .\broadcast.jl:709 [inlined]
[8] _broadcast_getindex
@ .\broadcast.jl:692 [inlined]
[9] #31
@ .\broadcast.jl:1118 [inlined]
[10] ntuple
@ .\ntuple.jl:48 [inlined]
[11] copy
@ .\broadcast.jl:1118 [inlined]
[12] materialize
@ .\broadcast.jl:903 [inlined]
[13] plan_th_leg2cheb!
@ C:\Users\danjv\.julia\packages\FastTransforms\S7huq\src\toeplitzhankel.jl:203 [inlined]
[14] plan_th_leg2cheb!
@ C:\Users\danjv\.julia\packages\FastTransforms\S7huq\src\toeplitzhankel.jl:651 [inlined]
[15] th_leg2cheb(::Vector{Float64})
@ FastTransforms C:\Users\danjv\.julia\packages\FastTransforms\S7huq\src\toeplitzhankel.jl:652
[16] top-level scope
@ REPL[20]:1
julia> FastTransforms.th_cheb2leg(Float64[])
ERROR: ArgumentError: invalid Array dimensions
Stacktrace:
[1] Array
@ .\boot.jl:477 [inlined]
[2] Array
@ .\boot.jl:486 [inlined]
[3] zeros
@ .\array.jl:636 [inlined]
[4] zeros
@ .\array.jl:632 [inlined]
[5] _cheb2legTH_TLC(::Type{Float64}, mn::Tuple{Int64}, d::Int64)
@ FastTransforms C:\Users\danjv\.julia\packages\FastTransforms\S7huq\src\toeplitzhankel.jl:218
[6] _broadcast_getindex_evalf
@ .\broadcast.jl:709 [inlined]
[7] _broadcast_getindex
@ .\broadcast.jl:692 [inlined]
[8] #31
@ .\broadcast.jl:1118 [inlined]
[9] ntuple
@ .\ntuple.jl:48 [inlined]
[10] copy
@ .\broadcast.jl:1118 [inlined]
[11] materialize
@ .\broadcast.jl:903 [inlined]
[12] plan_th_cheb2leg!
@ C:\Users\danjv\.julia\packages\FastTransforms\S7huq\src\toeplitzhankel.jl:235 [inlined]
[13] plan_th_cheb2leg!
@ C:\Users\danjv\.julia\packages\FastTransforms\S7huq\src\toeplitzhankel.jl:651 [inlined]
[14] th_cheb2leg(::Vector{Float64})
@ FastTransforms C:\Users\danjv\.julia\packages\FastTransforms\S7huq\src\toeplitzhankel.jl:652
[15] top-level scope
@ REPL[21]:1
- 主要言語
- Julia
- スター
- 282
- フォーク
- 27
- 平均マージ
- 53分
- マージ済み PR(30日)
- 1
環境構築
このプロジェクトには開発コンテナ、Dockerfile、コントリビューションガイドがありません。まず README を読み、一般的な手順ははじめてのコントリビューションガイドを参照してください。
はじめの一歩
- issue を最後まで読み、次にプロジェクトのコントリビューションガイドを読みます。
- 着手することを issue にコメントします — 二人が同じ作業をするのを防げます。
- リポジトリをフォークし、ブランチを切って変更します。
- issue 番号を参照したプルリクエストを送ります。
JuliaApproximation/FastTransforms.jl のほかの issue
-
Loading FastTransforms.jl can make FFTs via FFTW.jl 100x slower due to threading conflicts対応中かも @dlfivefifty が 96 日前に担当しました。 オープン
JuliaApproximation/FastTransforms.jl#267 · リアクション 2 件 · 担当者 2 名 ·
-
難易度 3/5 1〜2日 初心者へのやさしさ 38/100
JuliaApproximation/FastTransforms.jl#266 · コメント 1 件 ·
-
難易度 4/5 3〜5日 初心者へのやさしさ 42/100
JuliaApproximation/FastTransforms.jl#263 · コメント 2 件 ·
-
難易度 4/5 3〜5日 初心者へのやさしさ 35/100
JuliaApproximation/FastTransforms.jl#257 · コメント 1 件 ·
-
Allocating lmul!オープン
難易度 4/5 3〜5日 初心者へのやさしさ 25/100
JuliaApproximation/FastTransforms.jl#253 · コメント 7 件 ·
JuliaApproximation/FastTransforms.jl の issue をすべて見る
似ている issue
-
難易度 2/5 1時間未満 初心者へのやさしさ 78/100
メンテナーはふだん 1 日以内に返信
-
難易度 2/5 1〜3時間 初心者へのやさしさ 62/100
CliMA/ClimaAtmos.jl#4891 · コメント 1 件 ·
メンテナーはふだん 1 日以内に返信
-
tagbot-manual
難易度 1/5 1時間未満 初心者へのやさしさ 75/100
tshort/StaticCompiler.jl#191 ·
-
難易度 2/5 1〜3時間 初心者へのやさしさ 68/100
SciML/SciMLSensitivity.jl#1703 · コメント 1 件 ·
メンテナーはふだん 1 日以内に返信
-
難易度 2/5 1〜3時間 初心者へのやさしさ 68/100
JuliaGPU/GPUCompiler.jl#1008 ·
メンテナーはふだん 1 日以内に返信