Hacktoberfest 2026: le issue che i maintainer hanno segnato per ottobre, aperte e adatte ai principianti. Sfoglia le issue Hacktoberfest

cholesky decomposition

Aperta
#414 1 commento 0 reazioni 0 assegnatari Vedi su GitHub

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

Direzione di ricerca

Inizia con la rrule per cholesky in cui A <: LinearAlgebra.HermOrSym in src/rulesets/LinearAlgebra/factorization.jl intorno alla riga 452, quindi riproduci gli esempi di Zygote.jacobian dell’issue. Confronta il percorso Hermitian con i percorsi Matrix e cholesky semplice, inclusa la issue collegata di ChainRulesCore mul!. Il lavoro è completato quando il Jacobian di Hermitian cholesky non produce più zeri.

Scritto dal modello di indicizzazione a partire dal testo della issue.

Descrizione

I'm getting incorrect results when working with the rrule for cholesky where A <: LinearAlgebra.HermOrSym

Passing the input matrix through Matrix fixes the issue. The mul! fix relates to this issue.

using Zygote, ChainRules,LinearAlgebra
# Example matrix
A = [2. -1. 0.0; -1. 2. -1.; 0. -1. 2. ]

import LinearAlgebra.mul!
LinearAlgebra.mul!(C, ::ChainRulesCore.ZeroTangent, ::Any, ::Any, b) = C *=b

# produces zeros
Zygote.jacobian(a -> cholesky(Hermitian(a)).L , A)[1]

# both of the following produce expected result
Zygote.jacobian(a -> cholesky(Matrix(Hermitian(a))).L , A)[1]
Zygote.jacobian(a -> cholesky(a).L , A)[1]
Lingua principale
Julia
Stelle
475
Fork
98
Metriche di merge delle PR
Nessuna PR unita negli ultimi 30g

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

  1. Leggi tutta la issue e poi la guida ai contributi del progetto.
  2. Commenta sulla issue per dire che te ne occupi tu — evita che due persone facciano lo stesso lavoro.
  3. Fai un fork del repository e lavora su un branch.
  4. Apri una pull request che faccia riferimento al numero della issue.

Altre issue di JuliaDiff/ChainRules.jl

Tutte le issue di JuliaDiff/ChainRules.jl

Issue simili

Altre issue su Julia

Ricevi le nuove issue nella tua casella

Un breve riepilogo di issue GitHub adatte ai principianti.