Incorrect substitution of multivariate derivatives
还没有人认领这个 Issue。
评估
调研方向
先运行 issue 正文中的 Julia 复现，重点关注 dww 和 u0 附近的 diff 及替换调用。跟踪替换后如何计算多元导数；当列出的所有导数（包括总阶数为六的情况）在自身代入后都变为 0 时，工作即完成。
由索引模型根据 Issue 内容生成。
描述
using SymPy
@vars(u, w)
@symfuns(x)
dww = diff(x(u,w), (u, 1), (w, 2))(w=>0)
u0 = (
u=>0,
diff(x(u,w), (u, 2), (w, 4))(0,0)=>0,
)
diff(dww, (u, 3))(u0...)
for N in 1:6
t = diff(x(u, w), (u, N), (w, 1))
println(t, ':', t(0,0)(t(0,0)=>0))
end
gives the output
Derivative(x(u, w), u, w):0
Derivative(x(u, w), (u, 2), w):0
Derivative(x(u, w), (u, 3), w):0
Derivative(x(u, w), (u, 4), w):0
Derivative(x(u, w), (u, 5), w):Subs(Derivative(x(u, w), (u, 4), (w, 2)), (u, w), (0, 0))
Derivative(x(u, w), (u, 6), w):0
I expect them all to be 0, but for some reason the evaluation of a (u, 5), (w, 1) derivative and other multivariate derivatives of total order 6 (not shown here) do not self-substitute correctly after evaluating the line diff(dww, (u, 3))(u0...).
- 主要语言
- Julia
- 星标
- 282
- 派生
- 59
- PR 合并指标
- 30 天内没有已合并 PR
环境准备
这个项目没有提供开发容器、Dockerfile 或贡献指南，环境需要你自己搭建：先看它的 README，通用步骤见我们的新手贡献指南。
从这里开始
- 先读完整个 Issue，再读项目的贡献指南。
- 在 Issue 下留言说明你要接手 —— 这能避免两个人做同样的事。
- Fork 仓库，在一个分支上完成修改。
- 提交 Pull Request，并在描述里引用这个 Issue 编号。
JuliaPy/SymPy.jl 的其他 Issue
-
难度 4/5 3-5 天 新手友好度 35/100
-
难度 2/5 1-3 小时 新手友好度 45/100
-
难度 2/5 1-3 小时 新手友好度 35/100
-
难度 4/5 3-5 天 新手友好度 30/100
-
难度 4/5 3-5 天 新手友好度 35/100
相似的 Issue
-
Chains resumed from `initial_state` take `num_warmup + 1` warm-up steps可能已有人在做 @thevolatilebit 今天认领。 未关闭
难度 2/5 1-3 小时 新手友好度 80/100
TuringLang/AbstractMCMC.jl#220 ·
-
found-by-agent
难度 2/5 1-3 小时 新手友好度 68/100
exanauts/SparseDirectSolver.jl#92 ·
维护者通常 1 天内回复
-
`inv` of a dense matrix fails for arrays whose `parent` is not an array of the same kind可能已有人在做 @devmotion 今天认领。 未关闭
难度 2/5 1-3 小时 新手友好度 76/100
JuliaLang/LinearAlgebra.jl#1740 ·
维护者通常 2 天内回复
-
难度 2/5 1-3 小时 新手友好度 78/100
维护者通常 1 天内回复
-
难度 2/5 1-3 小时 新手友好度 90/100
NumericalEarth/Breeze.jl#1051 ·
维护者通常 1 天内回复