bugcategory:numericalhelp wanted
仓库指标
- 星标
- (15,065 个星标)
- PR 合并指标
- (PR 指标待抓取)
描述
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 :).