josdejong / josdejong/mathjs

Schur result not quasi upper triangular

Open
#3,600 3 comments 0 reactions 0 assignees View on GitHub
bug category:numerical help wanted
Dominant language
JavaScript
Stars
15.1k
Forks
1.3k
PR merge metrics
No merged PRs in 30d

Description

**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

Open the contributing guide

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.