Error with nested calls to TaylorDiff
Chưa có ai nhận issue này.
Đánh giá
- Độ khó
- 4/5
- Thời gian dự kiến
- 3-5 ngày
- Mức phù hợp với người mới
- 38/100
Hướng nghiên cứu
Tái hiện ví dụ đạo hàm lồng nhau từ issue bằng Julia 1.11.5 và TaylorDiff v0.3.3. Bắt đầu với derivative.jl và các định nghĩa phép nhân trong primitive.jl, đặc biệt là các method được báo cáo trong phần mơ hồ. Hoàn thành khi lời gọi TaylorDiff lồng nhau không còn phát sinh MethodError mơ hồ cho * và ví dụ có thể tính gradient được yêu cầu.
Do mô hình lập chỉ mục viết ra từ nội dung của issue.
Mô tả
Hi, I have the following nested AD problem that I'd like to get to work with TaylorDiff. Starting with a vector-valued function f(params, x) for scalar x, take a high-order derivative with respect to x and evaluate for a specific value of x. Then apply some reduction function to the result to obtain a scalar-valued function g(params). Finally I want to evaluate the gradient d g / d params. Example:
# Arbitrary function:
f(params, x) = [params[1] * x^3 + params[2], params[2] * sin(x - params[1]), sqrt(x + params[2])]
function g(params)
closure(x) = f(params, x)
some_x = 0.7
d3f_dx3 = TaylorDiff.derivative(closure, some_x, Val(3))
return sum(d3f_dx3)
end
some_params = [1.3, 2.1]
@show g(some_params) # Fine, gives 6.095380076578732
TaylorDiff.derivative(g, some_params, [1.0, 0.0], Val(1)) # First element of the gradient
Results, using julia 1.11.5 and TaylorDiff v0.3.3:
ERROR: MethodError: *(::TaylorScalar{Float64, 1}, ::TaylorScalar{Float64, 3}) is ambiguous.
Candidates:
*(a::TaylorScalar, b::Number)
@ TaylorDiff ~/.julia/packages/TaylorDiff/qw5aY/src/primitive.jl:119
*(a::Number, b::TaylorScalar)
@ TaylorDiff ~/.julia/packages/TaylorDiff/qw5aY/src/primitive.jl:114
Possible fix, define
*(::TaylorScalar, ::TaylorScalar)
Stacktrace:
[1] f(params::Vector{TaylorScalar{Float64, 1}}, x::TaylorScalar{Float64, 3})
@ Main ./REPL[4]:1
[2] (::var"#closure#1"{Vector{TaylorScalar{Float64, 1}}})(x::TaylorScalar{Float64, 3})
@ Main ./REPL[5]:2
[3] derivatives
@ ~/.julia/packages/TaylorDiff/qw5aY/src/derivative.jl:41 [inlined]
[4] derivative
@ ~/.julia/packages/TaylorDiff/qw5aY/src/derivative.jl:16 [inlined]
[5] g(params_in::Vector{TaylorScalar{Float64, 1}})
@ Main ./REPL[5]:4
[6] derivatives
@ ~/.julia/packages/TaylorDiff/qw5aY/src/derivative.jl:41 [inlined]
[7] derivative(f::Function, x::Vector{Float64}, l::Vector{Float64}, p::Val{1})
@ TaylorDiff ~/.julia/packages/TaylorDiff/qw5aY/src/derivative.jl:17
[8] top-level scope
@ REPL[8]:1
Any idea how this could be made to work?
While Zygote-over-TaylorDiff does work for this problem, @btime shows it is much faster to use ForwardDiff-over-ForwardDiff (probably due to the overhead of reverse mode), so I imagine TaylorDiff-over-TaylorDiff (or ForwardDiff-over-TaylorDiff) might be even faster due to the high-order inner derivative. Thanks.
- Ngôn ngữ chính
- Julia
- Star
- 83
- Fork
- 12
- Chỉ số merge pull request
- Không có pull request nào được merge trong 30 ngày
Chuẩn bị môi trường
Khởi chạy dev container của dự án ngay trên trình duyệt, bằng tài khoản GitHub của bạn.
- Không có Dockerfile hay tệp Docker Compose
- Không có mẫu pull request
- Không có hướng dẫn đóng góp
Bắt đầu từ đâu
- Đọc hết issue, rồi đọc hướng dẫn đóng góp của dự án.
- Bình luận trên issue rằng bạn sẽ nhận — tránh hai người làm cùng một việc.
- Fork repository và làm thay đổi trên một nhánh.
- Mở pull request có tham chiếu số hiệu của issue.
Issue khác của JuliaDiff/TaylorDiff.jl
-
Ambiguity with forwarddiff ?Đang mở
Độ khó 4/5 3-5 ngày Mức phù hợp với người mới 35/100
JuliaDiff/TaylorDiff.jl#103 · 6 bình luận ·
-
Độ khó 4/5 3-5 ngày Mức phù hợp với người mới 38/100
JuliaDiff/TaylorDiff.jl#102 · 3 bình luận ·
-
Độ khó 5/5 Hơn một tuần Mức phù hợp với người mới 25/100
JuliaDiff/TaylorDiff.jl#95 · 1 bình luận ·
-
Độ khó 4/5 3-5 ngày Mức phù hợp với người mới 35/100
JuliaDiff/TaylorDiff.jl#93 ·
-
Độ khó 4/5 3-5 ngày Mức phù hợp với người mới 25/100
JuliaDiff/TaylorDiff.jl#92 · 1 bình luận ·
Tất cả issue của JuliaDiff/TaylorDiff.jl
Issue tương tự
-
Độ khó 1/5 Dưới một giờ Mức phù hợp với người mới 85/100
-
Độ khó 2/5 1-3 giờ Mức phù hợp với người mới 83/100
JuliaGenAI/HuggingFaceDatasets.jl#96 ·
Maintainer thường phản hồi trong vòng 1 ngày
-
Độ khó 2/5 1-3 giờ Mức phù hợp với người mới 78/100
-
Độ khó 2/5 1-3 giờ Mức phù hợp với người mới 83/100
Maintainer thường phản hồi trong vòng 1 ngày
-
enhancement
Độ khó 2/5 1-3 giờ Mức phù hợp với người mới 76/100