Add multivariate root finding (Newton and Cuyt's Halley)
I maintainer di solito rispondono entro 1 giorno
Valutazione
- Difficoltà
- 5/5
- Tempo stimato
- Più di una settimana
- Idoneità per principianti
- 12/100
Direzione di ricerca
Start from the scalar solvers in tools/roots.hpp (Newton, Halley, Schröder) to see the existing iteration and convergence conventions, then check autodiff.hpp for how forward-mode derivatives are exposed. Design the multivariate Newton step J(x)a = −F(x) and the Cuyt–Halley update first, since the globalization choice and the fallback near aᵢ + bᵢ/2 ≈ 0 are open questions. Done means a documented API with tests on known systems that converge to the expected roots and show cubic convergence for the Halley variant.
Scritto dal modello di indicizzazione a partire dal testo della issue.
Descrizione
tools/roots.hpp has Newton, Halley and Schröder for scalar equations, but there is no solver for a system F(x) = 0 with x in ℝⁿ.
Proposal:
- Multivariate Newton: solve J(x) a = −F(x), with a line search or trust region for globalization.
- Multivariate Halley via Cuyt's Padé approach (Cuyt & Van der Cruyssen 1983; Cuyt, ACM TOMS 1985): with a = −J⁻¹F and b = J⁻¹F″[a, a], update componentwise xᵢ ← xᵢ + aᵢ² / (aᵢ + bᵢ/2). It converges cubically, reuses one factorization of J for both solves, and needs only the directional second derivative F″[a, a], not the full tensor. Fall back to the Newton component where aᵢ + bᵢ/2 ≈ 0.
- Derivatives from a user-supplied Jacobian, or from autodiff (forward-mode
autodiff.hpphandles F″[a, a] directly).
Uses include equilibria, collocation, closest points on implicit curves and surfaces, and ray–surface intersection in rendering.
- Lingua principale
- C++
- Stelle
- 377
- Fork
- 269
- Merge medio
- 12h 52m
- PR unite (30g)
- 60
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 boostorg/math
-
mode_has_nan gives MSVC warning C4702 (unreachable code)Forse già presa @ibmibmibm l’ha presa oggi. Aperta
Difficoltà 2/5 1-3 ore Idoneità per principianti 78/100
I maintainer di solito rispondono entro 1 giorno
-
Add a trust-region Newton optimizerForse già presa @NAThompson l’ha presa oggi. Aperta
Difficoltà 5/5 Più di una settimana Idoneità per principianti 15/100
I maintainer di solito rispondono entro 1 giorno
-
hypergeometric_0F1 loses about 20 ulps for b = 1.5 and −10 < z < 0Forse già presa @NAThompson l’ha presa oggi. Aperta
Difficoltà 4/5 3-5 giorni Idoneità per principianti 22/100
I maintainer di solito rispondono entro 1 giorno
-
Feature WishlistApertaenhancement
Difficoltà 5/5 Più di una settimana Idoneità per principianti 20/100
boostorg/math#303 · 18 commenti ·
I maintainer di solito rispondono entro 1 giorno
Tutte le issue di boostorg/math
Issue simili
-
Difficoltà 2/5 1-3 ore Idoneità per principianti 80/100
-
Difficoltà 1/5 Meno di un'ora Idoneità per principianti 82/100
-
Difficoltà 2/5 1-3 ore Idoneità per principianti 72/100
eclipse-score/lifecycle#929 ·
I maintainer di solito rispondono entro 1 giorno
-
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 68/100
colbymchenry/codegraph#3223 ·
I maintainer di solito rispondono entro 1 giorno