`WendlandKernel`: Compact support -> sparse kernel matrix
Nessuno ha ancora preso questa issue.
Valutazione
- Difficoltà
- 5/5
- Tempo stimato
- Più di una settimana
- Idoneità per principianti
- 25/100
- Tipo di issue
- Funzionalità
- Chiarezza
- Abbastanza chiara
- Stato di attività
- Ferma
- Stack tecnologico
- julia
- Ambito
- machine-learning
Direzione di ricerca
Inizia esaminando l’implementazione proposta di WendlandKernel in KernelFunctions.jl e le astrazioni kernel esistenti, quindi analizza come vengono rappresentate le matrici kernel sparse. Verifica l’interazione con AbstractGPs.jl e determina se il supporto per sparse Cholesky debba rientrare in quel progetto o in una modifica separata; il lavoro sarà considerato completato quando ci sarà accordo su questo ambito e sarà disponibile un’implementazione funzionante del kernel.
Scritto dal modello di indicizzazione a partire dal testo della issue.
Descrizione
Hey 👋
I think Wendland kernels are cool, and I would like to contribute a PR. Let me know what you think :)
What's this about?
Wendland kernels have compact support. If two points are "too far away" (quantified by the lengthscale of the kernel), their covariance under a Wendland kernel is zero. This results in sparse kernel matrices, which can be leveraged to save memory and compute. In particular, Julia has built-in support for highly optimised sparse Cholesky decompositions.
Within their support, the Wendland functions are defined by a rational polynomial, the coefficients of which can be computed in closed form.
At the same time, Wendland kernels have nice theoretical properties. In particular, much like the Matérn kernels, their smoothness is controllable directly through a smoothness parameter.
For more details, and in particular for the definition of Wendland functions, refer to Chapter 9 of Scattered Data Approximation by Holger Wendland.
What do I propose?
I would like to add the WendlandKernel to KernelFunctions.jl. I already have an implementation locally for Wendland kernels with arbitrary space dimension $d$ and smoothness parameter $k$. It produces sparse kernel matrices.
Difficulties
Currently, AbstractGPs.jl does not work with sparse Cholesky factorizations. I added some extensions locally, and I am not sure if these should go to KernelFunctions.jl or if I should make a separate PR for AbstractGPs.jl.
- Lingua principale
- Julia
- Stelle
- 276
- Fork
- 42
- Merge medio
- 1g 13h
- PR unite (30g)
- 1
Preparare l'ambiente
- Nessun Dockerfile né file Docker Compose
- Ha un modello di pull request
- Leggi la guida per i contributori
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 JuliaGaussianProcesses/KernelFunctions.jl
-
Difficoltà 2/5 1-3 ore Idoneità per principianti 75/100
JuliaGaussianProcesses/KernelFunctions.jl#609 · 1 commento ·
-
Future plans?Aperta
Difficoltà 5/5 Più di una settimana Idoneità per principianti 25/100
JuliaGaussianProcesses/KernelFunctions.jl#585 · 8 commenti · 3 reazioni ·
-
Difficoltà 4/5 3-5 giorni Idoneità per principianti 30/100
-
Difficoltà 3/5 1-2 giorni Idoneità per principianti 48/100
JuliaGaussianProcesses/KernelFunctions.jl#571 · 1 reazione ·
-
Difficoltà 2/5 1-3 ore Idoneità per principianti 25/100
JuliaGaussianProcesses/KernelFunctions.jl#533 · 1 reazione ·
Tutte le issue di JuliaGaussianProcesses/KernelFunctions.jl
Issue simili
-
broken links in docsAperta
Difficoltà 1/5 Meno di un'ora Idoneità per principianti 78/100
I maintainer di solito rispondono entro 1 giorno
-
bug
Difficoltà 2/5 1-3 ore Idoneità per principianti 76/100
JuliaPhysics/BeamletOptics.jl#127 ·
I maintainer di solito rispondono entro 1 giorno
-
Chains resumed from `initial_state` take `num_warmup + 1` warm-up stepsForse già presa @thevolatilebit l’ha presa oggi. Aperta
Difficoltà 2/5 1-3 ore Idoneità per principianti 80/100
TuringLang/AbstractMCMC.jl#220 ·
-
found-by-agent
Difficoltà 2/5 1-3 ore Idoneità per principianti 68/100
exanauts/SparseDirectSolver.jl#92 ·
I maintainer di solito rispondono entro 1 giorno
-
`inv` of a dense matrix fails for arrays whose `parent` is not an array of the same kindForse già presa @devmotion l’ha presa oggi. Aperta
Difficoltà 2/5 1-3 ore Idoneità per principianti 76/100
JuliaLang/LinearAlgebra.jl#1740 ·
I maintainer di solito rispondono entro 2 giorni