Forward difference gradient has the same number of function calls as central difference
Nessuno ha ancora preso questa issue.
Valutazione
- Difficoltà
- 4/5
- Tempo stimato
- 3-5 giorni
- Idoneità per principianti
- 35/100
- Tipo di issue
- Bug
- Chiarezza
- Abbastanza chiara
- Stato di attività
- Ferma
- Stack tecnologico
- julia
- Ambito
- performance
Direzione di ricerca
Inizia in src/gradients.jl, al costruttore di GradientCache e nel codice di finite_difference_gradient! intorno alle righe 138 e 144-145. Traccia come viene fornito fx per le differenze in avanti e confronta il comportamento delle chiamate di funzione con quello delle differenze centrali. Il lavoro è completato quando le differenze in avanti possono riutilizzare f(x) tramite la cache e richiedono meno chiamate di funzione, e sono comprese le implicazioni sulle prestazioni dell'approccio scelto.
Scritto dal modello di indicizzazione a partire dal testo della issue.
Descrizione
Changes to gradients.jl made in fc5a08e removed any way to reduce the number of function calls made when using forward difference for gradients as compared to central difference. The finite_difference_gradient! functions assume that f(x) will be provided in the GradientCache for forward differences, but the only constructor for GradientCache forces fx = nothing.
Here are the lines returning the only cache the constructor will give you:
https://github.com/JuliaDiff/FiniteDiff.jl/blob/9945e7d7f567dc2eb6a95a714b1de9782d0ecea4/src/gradients.jl#L52-L53
I can think of three ways to fix it:
- There could be a constructor that accepts an input of
fxreturns the a cache with fx, although fc5a08e was trying to get rid of extra constructors. - Instead of using the type
Nothingfor the first parameter of the returned GradientCache, useUnion{Nothing,returntype}so the user can manually update it later. - Or, don't bother using the cache and change the following line
https://github.com/JuliaDiff/FiniteDiff.jl/blob/9945e7d7f567dc2eb6a95a714b1de9782d0ecea4/src/gradients.jl#L138
to_fx, c1, c2, c3 = cache.fx, cache.c1, cache.c2, cache.c3and then before this for loop:
https://github.com/JuliaDiff/FiniteDiff.jl/blob/9945e7d7f567dc2eb6a95a714b1de9782d0ecea4/src/gradients.jl#L144-L145
add something along the lines offx = _fx == Nothing ? f(x) : _fx
2 makes the most sense to me, but I'm not sure what the performance implications are.
- Lingua principale
- Julia
- Stelle
- 305
- Fork
- 42
- Merge medio
- 24m
- PR unite (30g)
- 2
Preparare l'ambiente
Questo progetto non fornisce container di sviluppo, Dockerfile né guida per i contributori, quindi l'ambiente è a tuo carico: parti dal suo README e consulta la nostra guida al primo contributo per i passaggi generali.
Come iniziare
- Leggi tutta la issue e poi la guida ai contributi del progetto.
- Commenta sulla issue per dire che te ne occupi tu — evita che due persone facciano lo stesso lavoro.
- Fai un fork del repository e lavora su un branch.
- Apri una pull request che faccia riferimento al numero della issue.
Altre issue di JuliaDiff/FiniteDiff.jl
-
tagbot-manual
Difficoltà 2/5 1-3 ore Idoneità per principianti 55/100
JuliaDiff/FiniteDiff.jl#218 ·
-
Feature request: vectorized finite difference evaluationForse già presa @ChrisRackauckas-Claude l’ha presa 191 giorni fa. Aperta
Difficoltà 5/5 Più di una settimana Idoneità per principianti 45/100
JuliaDiff/FiniteDiff.jl#210 · 1 commento · 1 reazione ·
-
Difficoltà 3/5 1-2 giorni Idoneità per principianti 45/100
JuliaDiff/FiniteDiff.jl#199 · 3 commenti ·
-
StridedArraysAperta
Difficoltà 4/5 3-5 giorni Idoneità per principianti 25/100
JuliaDiff/FiniteDiff.jl#158 · 1 commento ·
-
Difficoltà 5/5 Più di una settimana Idoneità per principianti 25/100
JuliaDiff/FiniteDiff.jl#157 · 1 commento ·
Tutte le issue di JuliaDiff/FiniteDiff.jl
Issue simili
-
Difficoltà 2/5 1-3 ore Idoneità per principianti 70/100
SciML/OrdinaryDiffEq.jl#4875 ·
I maintainer di solito rispondono entro 3 giorni
-
documentation
Difficoltà 2/5 1-3 ore Idoneità per principianti 62/100
I maintainer di solito rispondono entro 1 giorno
-
Difficoltà 2/5 1-3 ore Idoneità per principianti 66/100
SciML/NonlinearSolve.jl#1372 ·
I maintainer di solito rispondono entro 1 giorno
-
Truncated discrete ccdf and logccdf return incorrect values at the lower atomForse già presa Una pull request collegata a questa issue è aperta o già unita. Aperta
Difficoltà 2/5 1-3 ore Idoneità per principianti 82/100
JuliaStats/Distributions.jl#2095 ·
-
documentation 📜 microphysics💧 ❄️
Difficoltà 2/5 1-3 ore Idoneità per principianti 63/100
NumericalEarth/Breeze.jl#1073 ·
I maintainer di solito rispondono entro 1 giorno