Tearing.Dispersion - Update the root polish's Jacobian by Broyden instead of two differences per step

Abierto
  • ¿Algo no está bien?
#604 0 comentarios 0 reacciones 0 asignados Ver en GitHub

Los mantenedores suelen responder en 3 días

Nadie ha tomado este issue todavía.

Evaluación

Dificultad
4/5
Tiempo estimado
3-5 días
Aptitud para principiantes
52/100
Tipo de issue
Refactorización
Claridad
Bien especificado
Estado de actividad
Activo
Stack tecnológico
julia
Área
performance

Línea de trabajo

Read _polish_root in Tearing/Dispersion/GrowthRateExtraction.jl and trace how it evaluates the residual, accepts trial steps, and handles up to three backtracks. Update the Jacobian after accepted steps using the step and residual change, and rebuild it from differences after a failed backtrack. Run the harness and compare polished roots with the existing tolerance.

Escrito por el modelo de indexación a partir del texto del issue.

Descripción

The root polish (_polish_root in Tearing/Dispersion/GrowthRateExtraction.jl) builds the real 2×2 Jacobian of (Re f, Im f) from two forward differences on every iteration. With the trial step and up to three backtracks that is 3–6 residual evaluations per iteration, and on the SLAYER path each evaluation is a stiff Riccati solve.

The cost is small today (diiid_slayer_n1 127 → 128 s), but it scales with the number of intersections polished and with anything that re-polishes, such as more surfaces or a finer scan.

Proposal: build the Jacobian once from two differences, then update it by Broyden's rank-one rule with the step and residual change already computed, so each iteration costs one evaluation. Rebuild it from differences only after a failed backtrack. Convergence and stopping criteria stay as they are, so roots should agree with the current polish to its tolerance; check with the harness.

Lenguaje dominante
Julia
Estrellas
11
Forks
7
Merge medio
3 d 14 h
PR fusionados (30 d)
69

Preparar el entorno

  • Sin Dockerfile ni archivo de Docker Compose
  • Tiene una plantilla de pull request
  • Sin guía de contribución

Primeros pasos

  1. Lee el issue completo y luego la guía de contribución del proyecto.
  2. Comenta en el issue que vas a ocuparte — evita que dos personas hagan lo mismo.
  3. Haz un fork del repositorio y trabaja en una rama.
  4. Abre un pull request que haga referencia al número del issue.

Más de OpenFUSIONToolkit/GPEC

Todos los issues de OpenFUSIONToolkit/GPEC

Issues similares

Más issues de Julia

Recibe los nuevos issues en tu correo

Un resumen breve de issues de GitHub para principiantes.