josdejong/mathjs

Schur result not quasi upper triangular

Open

#3.600 geöffnet am 22. Nov. 2025

Auf GitHub ansehen
 (3 Kommentare) (0 Reaktionen) (0 zugewiesene Personen)JavaScript (13.832 Stars) (1.298 Forks)batch import
bugcategory:numericalhelp wanted

Beschreibung

Describe the bug A clear and concise description of what the bug is.

I think sometimes Schur returns a matrix that isn't quasi upper triangular.

To Reproduce Generic orthonormal matrices seem to create a problem. For example:

math.matrix([
      [-0.03591206220229135, -0.09100469507870354, 0.9952027277203429],
      [-0.3802068171617618, -0.9197139315803332, -0.09782157349362577],
      [0.9242040358990549, -0.38189583596928406, -0.0015717815434243287]
    ])

Possibly mathjs is running QR without a Wilkinson/Francis shift strategy?

Happy to contribute a PR on this if I'm not hallucinating the problem :).

Contributor Guide

Tech Stack
javascript
Domain
data
Issue Type
bug
SchwierigkeitGeschätzte Implementierungsschwierigkeit für neue Contributors, von 1 für sehr kleine Änderungen bis 5 für Expertenarbeit.
3
Geschätzte ZeitEin grober Zeitrahmen, um zu recherchieren, zu implementieren, zu testen und einen Pull Request vorzubereiten.
1-2 days
AktivitätsstatusWie verfügbar das Issue gerade wirkt: frisch, aktiv, stale, blockiert oder wartet auf Maintainer-Input.
active
KlarheitWie klar das Issue die erwartete Änderung, Akzeptanzkriterien und den nächsten Schritt erklärt.
clear
Voraussetzungen
Numerical linear algebramathjs code structure
EinsteigerfreundlichkeitEin Score von 1 bis 100, der schätzt, wie zugänglich dieses Issue für First-Time Contributors ist.
65
Research-Richtung
Investigate the Schur decomposition implementation in mathjs, likely in the linear algebra module (e.g., src/function/algebra/schur.js or similar). Check if the QR algorithm used includes a shift strategy (e.g., Wilkinson shift or Francis double shift). The issue author suggests that the absence of such shift may cause the result to not be quasi upper triangular. Consider implementing the shift or verifying the existing code. The linked issue comments may provide additional clues.