Error when computing jacobian for broadcasted function
还没有人认领这个 Issue。
评估
- 难度
- 3/5
- 预计耗时
- 1-2 天
- 新手友好度
- 35/100
- Issue 类型
- 缺陷
- 描述清晰度
- 基本清楚
- 活跃度
- 停滞
- 技术栈
- julia
调研方向
使用 Julia 1.10.0-beta2 以及列出的 AbstractDifferentiation 和 Diffractor 版本重现三个 jacobian 调用。从 src/stage1/broadcast.jl:27 以及 src/tangent.jl:351 和 :375 处的 createinstance/rebundle 路径开始,然后比较正常工作和失败的 broadcast 案例。逐元素 broadcast 的 sin 表达式能够计算其 jacobian,且不再出现报告的 TangentBundle MethodError,即表示完成。
由索引模型根据 Issue 内容生成。
描述
Simple broadcasting errors for me. See below.
I am on Julia 1.10.0-beta2+0.x64.linux.gnu
within the following environment:
Status ~/.julia/environments/v1.10/Project.toml
[c29ec348] AbstractDifferentiation v0.5.2
[9f5e2b26] Diffractor v0.2.1
[31c24e10] Distributions v0.25.100
This shows the issue (verified that I get the same problem with gradient from AbstractDifferentiation.jl):
using Diffractor: DiffractorForwardBackend
using AbstractDifferentiation: derivative, jacobian
jacobian(DiffractorForwardBackend(), θ -> sin(θ[1]) + sin(θ[2]), [0.3, 0.8]) # works
jacobian(DiffractorForwardBackend(), θ -> sum(sin, θ), [0.3, 0.8]) # works
jacobian(DiffractorForwardBackend(), θ -> sum(sin.(θ)), [0.3, 0.8]) # errors
last command errors with (truncated error):
ERROR: MethodError: no method matching Diffractor.TangentBundle{1, Float64, Diffractor.TaylorTangent{Tuple{Float64}}}(::Float64, ::Tuple{Float64})
Stacktrace:
[1] createinstance(::Type{Diffractor.TangentBundle{1, Float64, Diffractor.TaylorTangent{Tuple{Float64}}}}, ::Float64, ::Vararg{Float64})
@ Diffractor ~/.julia/packages/Diffractor/vzpvE/src/tangent.jl:351
[2] _getindex
@ Diffractor ~/.julia/packages/StructArrays/dNQpc/src/structarray.jl:353 [inlined]
[3] getindex
@ Diffractor ~/.julia/packages/StructArrays/dNQpc/src/structarray.jl:348 [inlined]
[4] iterate
@ Diffractor ./abstractarray.jl:1212 [inlined]
[5] iterate
@ Diffractor ./abstractarray.jl:1210 [inlined]
[6] iterate
@ Diffractor ./generator.jl:44 [inlined]
[7] _collect
@ Diffractor ./array.jl:852 [inlined]
[8] collect_similar
@ Diffractor ./array.jl:761 [inlined]
[9] map
@ Diffractor ./abstractarray.jl:3273 [inlined]
[10] rebundle(A::StructArrays.StructVector{Diffractor.TangentBundle{1, Float64, Diffractor.TaylorTangent{Tuple{Float64}}}, Tuple{Vector{Float64}, Vector{Float64}}, Int64})
@ Diffractor ~/.julia/packages/Diffractor/vzpvE/src/tangent.jl:375
[11] (::Diffractor.∂☆{1})(zc::ZeroBundle{1, typeof(copy)}, bc::Diffractor.CompositeBundle{1, Base.Broadcast.Broadcasted{…}, Tuple{…}})
@ Diffractor ~/.julia/packages/Diffractor/vzpvE/src/stage1/broadcast.jl:27
[12] materialize
@ Diffractor ./broadcast.jl:903 [inlined]
[13] (::Diffractor.∂☆internal{1})(::ZeroBundle{1, Type{Base.Broadcast.Broadcasted}}, ::Vararg{Diffractor.AbstractTangentBundle{1}})
@ Diffractor ~/.julia/packages/Diffractor/vzpvE/src/stage1/forward.jl:121 [inlined]
[14] ∂☆
@ Diffractor ~/.julia/packages/Diffractor/vzpvE/src/stage1/forward.jl:160 [inlined]
[15] #61
@ Diffractor ./REPL[45]:1 [inlined]
[16] (::Diffractor.∂☆recurse{1})(::ZeroBundle{1, var"#61#62"}, ::Diffractor.TangentBundle{1, Vector{Float64}, Diffractor.TaylorTangent{Tuple{Vector{Float64}}}})
@ Diffractor ~/.julia/packages/Diffractor/vzpvE/src/stage1/recurse_fwd.jl:0
[17] (::Diffractor.∂☆internal{1})(::ZeroBundle{1, var"#61#62"}, ::Vararg{Diffractor.AbstractTangentBundle{1}})
@ Diffractor ~/.julia/packages/Diffractor/vzpvE/src/stage1/forward.jl:121
[18] (::Diffractor.∂☆{1})(::ZeroBundle{1, var"#61#62"}, ::Vararg{Diffractor.AbstractTangentBundle{1}})
@ Diffractor ~/.julia/packages/Diffractor/vzpvE/src/stage1/forward.jl:160
[19] (::Diffractor.var"#pushforward#359"{var"#61#62", Tuple{Vector{Float64}}})(vs::Tuple{Vector{Float64}})
@ Diffractor ~/.julia/packages/Diffractor/vzpvE/src/AbstractDifferentiation.jl:13
[20] (::Diffractor.var"#358#362"{Diffractor.var"#pushforward#359"{var"#61#62", Tuple{Vector{Float64}}}})(cols::Vector{Float64})
@ Diffractor ~/.julia/packages/AbstractDifferentiation/eEkWP/src/AbstractDifferentiation.jl:526
[21] mapslices(f::Diffractor.var"#358#362"{Diffractor.var"#pushforward#359"{var"#61#62", Tuple{Vector{Float64}}}}, A::Matrix{Float64}; dims::Int64)
@ Base ./abstractarray.jl:3179
[22] jacobian(b::DiffractorForwardBackend, f::Function, args::Vector{Float64})
@ Diffractor ~/.julia/packages/AbstractDifferentiation/eEkWP/src/AbstractDifferentiation.jl:524
[23] top-level scope
@ REPL[45]:1
Some type information was truncated. Use show(err) to see complete types.
- 主要语言
- Julia
- 星标
- 453
- 派生
- 33
- PR 合并指标
- 30 天内没有已合并 PR
环境准备
这个项目没有提供开发容器、Dockerfile 或贡献指南,环境需要你自己搭建:先看它的 README,通用步骤见我们的新手贡献指南。
从这里开始
- 先读完整个 Issue,再读项目的贡献指南。
- 在 Issue 下留言说明你要接手 —— 这能避免两个人做同样的事。
- Fork 仓库,在一个分支上完成修改。
- 提交 Pull Request,并在描述里引用这个 Issue 编号。
JuliaDiff/Diffractor.jl 的其他 Issue
-
难度 2/5 1-3 小时 新手友好度 25/100
JuliaDiff/Diffractor.jl#302 · 2 条评论 ·
-
难度 3/5 1-2 天 新手友好度 52/100
JuliaDiff/Diffractor.jl#294 ·
-
bug
难度 4/5 3-5 天 新手友好度 35/100
JuliaDiff/Diffractor.jl#293 ·
-
bug
难度 4/5 3-5 天 新手友好度 35/100
JuliaDiff/Diffractor.jl#290 · 8 条评论 ·
-
难度 3/5 1-2 天 新手友好度 35/100
JuliaDiff/Diffractor.jl#280 · 1 条评论 ·
查看 JuliaDiff/Diffractor.jl 的全部 Issue
相似的 Issue
-
难度 2/5 1-3 小时 新手友好度 83/100
JuliaGenAI/HuggingFaceDatasets.jl#96 ·
维护者通常 1 天内回复
-
难度 2/5 1-3 小时 新手友好度 66/100
JuliaData/DBFTables.jl#50 ·
-
code-review
难度 2/5 1-3 小时 新手友好度 72/100
-
难度 2/5 1-3 小时 新手友好度 72/100
JuliaNLSolvers/Optim.jl#1296 ·
-
OptimizationCUTEst jobs fail: gfortran not installed on the runner可能已有人在做 @ChrisRackauckas-Claude 于 1 天前认领。 未关闭
难度 2/5 1-3 小时 新手友好度 68/100
SciML/SciMLBenchmarks.jl#1946 · 1 条评论 ·
维护者通常 1 天内回复