Validated transforms

Offen
  • Stimmt etwas nicht?
#94 0 Kommentare 0 Reaktionen 0 zugewiesene Personen Auf GitHub ansehen

Dieses Issue hat noch niemand übernommen.

Bewertung

Schwierigkeit
5/5
Geschätzter Aufwand
Über eine Woche
Anfängerfreundlichkeit
30/100
Issue-Typ
Feature
Klarheit
Muss geklärt werden
Aktivitätsstatus
Veraltet
Tech-Stack
julia
Bereich
backend

Rechercherichtung

Beginne damit, das BigFloat RoundDown/RoundUp-Beispiel aus dem Issue zu reproduzieren, und untersuche anschließend src/transforms_mpfr.c und dessen ccall-Arithmetik. Prüfe die beteiligten Operationen, einschließlich des gemeldeten Falls mpfr_neg, und verifiziere, dass validierte Transforms unter beiden Rundungsmodi ohne den Token-Fehler korrekt funktionieren.

Vom Indexierungsmodell aus dem Issue-Text verfasst.

Beschreibung

beginner

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>
Vorherrschende Sprache
Julia
Sterne
284
Forks
28
Ø Merge
17 Std. 49 Min.
Gemergte PRs (30 T.)
1

Entwicklungsumgebung

Dieses Projekt bietet weder Dev-Container noch Dockerfile noch Beitragsleitfaden – die Einrichtung liegt bei Ihnen. Beginnen Sie mit der README; die allgemeinen Schritte stehen in unserem Leitfaden für den ersten Beitrag.

Erste Schritte

  1. Lesen Sie das ganze Issue und danach den Beitragsleitfaden des Projekts.
  2. Schreiben Sie ins Issue, dass Sie es übernehmen — das erspart doppelte Arbeit.
  3. Forken Sie das Repository und arbeiten Sie in einem Branch.
  4. Öffnen Sie einen Pull Request, der die Issue-Nummer nennt.

Mehr aus JuliaApproximation/FastTransforms.jl

Alle Issues in JuliaApproximation/FastTransforms.jl

Ähnliche Issues

Weitere Issues zu Julia

Neue Issues direkt in Ihr Postfach

Eine kurze Übersicht über anfängerfreundliche GitHub-Issues.