Incorrect Jacobian calculation
还没有人认领这个 Issue。
评估
调研方向
首先运行 issue 中的 Julia 复现程序,并将 DiffractorForwardBackend 的 Jacobian 与 AD.ForwardDiffBackend 以及 GTPSA.jl 的交叉检查结果进行比较。跟踪测试函数的 Diffractor Jacobian 路径,并确定其最初的对角线元素为何不同;当 Diffractor 结果与交叉检查得到的 Jacobian 一致时,即表示完成。
由索引模型根据 Issue 内容生成。
描述
using ForwardDiff
using Diffractor: DiffractorForwardBackend
import AbstractDifferentiation as AD
function test(z0)
z = Vector{eltype(z0)}(undef, length(z0))
L = 0.5
k1 = 0.36
z[1] = cos(sqrt(k1)*L)*z0[1]+1/sqrt(k1)*sin(sqrt(k1)*L)*z0[2]
z[2] = -sqrt(k1)*sin(sqrt(k1)*L)*z0[1]+cos(sqrt(k1)*L)*z0[2]
return z
end
m(z) = test([z[1], z[2]])
Then:
julia> j = AD.jacobian(DiffractorForwardBackend(), m, zeros(2)) |> only
2×2 Matrix{Float64}:
0.29552 0.492534
-0.177312 0.29552
julia> j = AD.jacobian(AD.ForwardDiffBackend(), m, zeros(2)) |> only
2×2 Matrix{Float64}:
0.955336 0.492534
-0.177312 0.955336
ForwardDiff.jl's Jacobian was also cross-checked with GTPSA.jl
- 主要语言
- 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#290 · 8 条评论 ·
-
难度 3/5 1-2 天 新手友好度 35/100
JuliaDiff/Diffractor.jl#280 · 1 条评论 ·
-
难度 4/5 3-5 天 新手友好度 45/100
JuliaDiff/Diffractor.jl#278 · 1 条评论 · 1 个 reaction ·
查看 JuliaDiff/Diffractor.jl 的全部 Issue
相似的 Issue
-
enhancement good first issue help wanted
难度 2/5 半天 新手友好度 82/100
-
authority-vechtstromen ready-for-agent
难度 1/5 1 小时以内 新手友好度 75/100
Deltares/Ribasim-NL#890 ·
维护者通常 1 天内回复
-
难度 1/5 1 小时以内 新手友好度 88/100
维护者通常 1 天内回复
-
难度 2/5 1-3 小时 新手友好度 78/100
JuliaLinearAlgebra/GenericLinearAlgebra.jl#183 · 1 条评论 ·
-
难度 2/5 1-3 小时 新手友好度 72/100
JuliaGPU/KernelAbstractions.jl#897 ·
维护者通常 1 天内回复