Validated transforms
Nessuno ha ancora preso questa issue.
Valutazione
- Difficoltà
- 5/5
- Tempo stimato
- Più di una settimana
- Idoneità per principianti
- 30/100
Direzione di ricerca
Inizia riproducendo l’esempio BigFloat RoundDown/RoundUp nell’issue, quindi esamina src/transforms_mpfr.c e la sua aritmetica ccall. Verifica le operazioni coinvolte, incluso il caso mpfr_neg segnalato, e assicurati che i transforms validati si comportino correttamente con entrambe le modalità di arrotondamento senza l’errore di token.
Scritto dal modello di indicizzazione a partire dal testo della issue.
Descrizione
This would be fun. As seen below, it almost "just works," but would require a careful investigation into https://github.com/MikaelSlevinsky/FastTransforms/blob/master/src/transforms_mpfr.c. The ccalls pass on Base.MPFR.RoundingMode[] to the C library, so it's just a matter of making sure every arithmetic operation is appropriate for RoundUp/Down. For example, I think the token error is due to using mpfr_neg.
julia> using FastTransforms
julia> setrounding(BigFloat, RoundDown)
MPFRRoundDown::MPFRRoundingMode = 3
julia> pd = plan_leg2cheb(BigFloat, 10)
FastTransforms Legendre--Chebyshev plan for 10-element array of BigFloat
julia> Pd = lmul!(pd, Matrix{BigFloat}(I, 10, 10));
julia> setrounding(BigFloat, RoundUp)
MPFRRoundUp::MPFRRoundingMode = 2
julia> pu = plan_leg2cheb(BigFloat, 10)
FastTransforms Legendre--Chebyshev plan for 10-element array of BigFloat
julia> Pu = lmul!(pu, Matrix{BigFloat}(I, 10, 10));
julia> Pu-Pd # rigorously off by only one entry
10×10 Array{BigFloat,2}:
0.0 0.0 -2.15904e-78 0.0 2.15904e-78 0.0 2.15904e-78 0.0 3.23856e-78 0.0
0.0 0.0 0.0 8.63617e-78 0.0 8.63617e-78 0.0 1.07952e-77 0.0 8.63617e-78
0.0 0.0 0.0 0.0 8.63617e-78 0.0 8.63617e-78 0.0 1.07952e-77 0.0
0.0 0.0 0.0 1.72723e-77 0.0 1.29543e-77 0.0 1.29543e-77 0.0 1.29543e-77
0.0 0.0 0.0 0.0 1.72723e-77 0.0 1.72723e-77 0.0 1.51133e-77 0.0
0.0 0.0 0.0 0.0 0.0 2.59085e-77 0.0 1.94314e-77 0.0 1.51133e-77
0.0 0.0 0.0 0.0 0.0 0.0 3.45447e-77 0.0 1.94314e-77 0.0
0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.88628e-77 0.0 2.15904e-77
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.88628e-77 0.0
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 4.31808e-77
julia>
- Lingua principale
- Julia
- Stelle
- 282
- Fork
- 27
- Merge medio
- 53m
- PR unite (30g)
- 1
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 JuliaApproximation/FastTransforms.jl
-
Loading FastTransforms.jl can make FFTs via FFTW.jl 100x slower due to threading conflictsForse già presa @dlfivefifty l’ha presa 90 giorni fa. Aperta
JuliaApproximation/FastTransforms.jl#267 · 2 reazioni · 2 assegnatari ·
-
Difficoltà 3/5 1-2 giorni Idoneità per principianti 38/100
JuliaApproximation/FastTransforms.jl#266 · 1 commento ·
-
Difficoltà 4/5 3-5 giorni Idoneità per principianti 42/100
JuliaApproximation/FastTransforms.jl#263 · 2 commenti ·
-
Difficoltà 4/5 3-5 giorni Idoneità per principianti 35/100
JuliaApproximation/FastTransforms.jl#257 · 1 commento ·
-
Allocating lmul!Aperta
Difficoltà 4/5 3-5 giorni Idoneità per principianti 25/100
JuliaApproximation/FastTransforms.jl#253 · 7 commenti ·
Tutte le issue di JuliaApproximation/FastTransforms.jl
Issue simili
-
Difficoltà 2/5 1-3 ore Idoneità per principianti 72/100
oxfordcontrol/COSMO.jl#211 ·
-
documentation
Difficoltà 2/5 Mezza giornata Idoneità per principianti 65/100
I maintainer di solito rispondono entro 6 giorni
-
Out-of-place JLArray/GPU problem with VectorContinuousCallback scalar-indexes (callback cache built with CPU zeros)Forse già presa @ChrisRackauckas-Claude l’ha presa oggi. Aperta
Difficoltà 2/5 1-3 ore Idoneità per principianti 74/100
SciML/OrdinaryDiffEq.jl#4813 ·
I maintainer di solito rispondono entro 1 giorno
-
ARKODE: callbacks that modify `u` throw MethodError on reinitForse già presa @devmotion l’ha presa 1 giorno fa. Aperta
Difficoltà 2/5 1-3 ore Idoneità per principianti 79/100
SciML/Sundials.jl#575 ·
-
broken links in docsAperta
Difficoltà 1/5 Meno di un'ora Idoneità per principianti 78/100
I maintainer di solito rispondono entro 1 giorno