Incorrect Gradient Computation when Nesting Zygote over TaylorDiff Description
Nadie ha tomado este issue todavía.
Evaluación
- Dificultad
- 4/5
- Tiempo estimado
- 3-5 días
- Aptitud para principiantes
- 25/100
Línea de trabajo
Empieza ejecutando el ejemplo mínimo con las versiones indicadas de Julia, TaylorDiff y Zygote, centrándote en TaylorDiff.derivative anidado dentro de Zygote.gradient. Sigue la interacción entre esos puntos de entrada y añade una prueba de regresión para el ejemplo. Se considera completado cuando el gradiente calculado ya no es todo ceros y coincide con el resultado esperado distinto de cero.
Escrito por el modelo de indexación a partir del texto del issue.
Descripción
Bug Report: Incorrect Gradient Computation when Nesting Zygote over TaylorDiff
Description
When attempting to compute the gradient of a function that internally uses TaylorDiff for differentiation, Zygote returns incorrect results (all zeros).
Minimal Working Example
The following mwe demonstrates the issue:
import TaylorDiff
import Zygote
# Define a function
f(x) = sum(x .^ 2)
# Define function to return unit vectors
unit_vectors(x, i) = [j == i ? one(eltype(x)) : zero(eltype(x)) for j in 1:length(x)]
# Define input
x = Float64.(collect(1:5))
# Compute unit vectors for each direction to compute the gradient
e_vectors = [unit_vectors(x, i) for i in 1:length(x)]
# Define a function to compute gradient using TaylorDiff
∇f_taylor(f, x) = TaylorDiff.derivative.(Ref(f), Ref(x), e_vectors, Ref(Val(1)))
# Evaluate gradient
∇f_taylor(f, x) # Returns 5-element Vector{Float64}: 2.0 4.0 6.0 8.0 10.0 (correct)
# Define a simple function of the gradient with TaylorDiff
g_taylor(x) = sum(∇f_taylor(f, x))
# Evaluate gradient
g_taylor(x) # returns 30.0 (Correct)
# Compute gradient using Zygote
Zygote.gradient(g_taylor, x) # returns ([0.0, 0.0, 0.0, 0.0, 0.0],) (incorrect)
The function g_taylor(x) computes the sum of the gradient of f(x), which should return a nonzero result. Applying Zygote.gradient to g_taylor(x) should correctly compute the derivative.
Instead of returning the correct gradient, Zygote.gradient(g_taylor, x) returns all zeros.
Additional Information
Julia Version: 1.11.3
TaylorDiff Version: v0.3.1
Zygote Version: v0.6.75
- Lenguaje dominante
- Julia
- Estrellas
- 83
- Forks
- 12
- Métricas de merge de PR
- Sin PR fusionados en 30 d
Preparar el entorno
Inicia el contenedor de desarrollo del proyecto en tu navegador, con tu propia cuenta de GitHub.
- Sin Dockerfile ni archivo de Docker Compose
- Sin plantilla de pull request
- Sin guía de contribución
Primeros pasos
- Lee el issue completo y luego la guía de contribución del proyecto.
- Comenta en el issue que vas a ocuparte — evita que dos personas hagan lo mismo.
- Haz un fork del repositorio y trabaja en una rama.
- Abre un pull request que haga referencia al número del issue.
Más de JuliaDiff/TaylorDiff.jl
-
Ambiguity with forwarddiff ?Abierto
Dificultad 4/5 3-5 días Aptitud para principiantes 35/100
JuliaDiff/TaylorDiff.jl#103 · 6 comentarios ·
-
Dificultad 4/5 3-5 días Aptitud para principiantes 38/100
JuliaDiff/TaylorDiff.jl#102 · 3 comentarios ·
-
Dificultad 4/5 3-5 días Aptitud para principiantes 38/100
JuliaDiff/TaylorDiff.jl#99 · 3 comentarios ·
-
Dificultad 5/5 Más de una semana Aptitud para principiantes 25/100
JuliaDiff/TaylorDiff.jl#95 · 1 comentario ·
-
Dificultad 4/5 3-5 días Aptitud para principiantes 35/100
JuliaDiff/TaylorDiff.jl#93 ·
Todos los issues de JuliaDiff/TaylorDiff.jl
Issues similares
-
Add more snippets to R and JuliaAbiertoenhancement good first issue help wanted
Dificultad 2/5 Medio día Aptitud para principiantes 82/100
-
Dificultad 2/5 1-3 horas Aptitud para principiantes 78/100
SciML/SciMLBase.jl#1643 ·
Los mantenedores suelen responder en 2 días
-
authority-vechtstromen ready-for-agent
Dificultad 1/5 Menos de una hora Aptitud para principiantes 75/100
Deltares/Ribasim-NL#890 ·
Los mantenedores suelen responder en 1 día
-
Dificultad 1/5 Menos de una hora Aptitud para principiantes 88/100
Los mantenedores suelen responder en 1 día
-
Precompilation fails on Julia 1.14 nightly: HessenbergQ lmul! methods overwrite LinearAlgebra'sAbierto
Dificultad 2/5 1-3 horas Aptitud para principiantes 78/100
JuliaLinearAlgebra/GenericLinearAlgebra.jl#183 · 1 comentario ·