Skip to content

Conversation

@ev-br
Copy link
Member

@ev-br ev-br commented Mar 6, 2024

Reference issue

closes #17125

What does this implement/fix?

The failure is relatively benign if annoying: the computation is most likely always correct, what fails is ordering of eigenvalues consistently on the complex plane, as Warren's analysis in #17125 (comment) shows.

So try a bit harder to trim numerical noise and account for complex conjugate pairs.

Have to admit I cannot repro locally. @h-vetinari you mentioned it's reproducible on conda-forge? #17125 (comment)

Additional information

@ev-br ev-br added scipy.linalg maintenance Items related to regular maintenance tasks labels Mar 6, 2024
@ev-br ev-br requested review from ilayn and larsoner as code owners March 6, 2024 12:33
@ev-br ev-br added defect A clear bug or issue that prevents SciPy from being installed or used as expected and removed maintenance Items related to regular maintenance tasks labels Mar 6, 2024
@h-vetinari
Copy link
Member

@h-vetinari you mentioned it's reproducible on conda-forge? #17125 (comment)

I realized that I never finished the cross-BLAS-flavour analysis for 1.12. I've restarted that now, and will try to write a summary later. If you want to check yourself, you should see1 it in the win_64_blas_implopenblas... builds, if the failure is still there.

Footnotes

  1. modulo getting the right agent in CI matching avx512 absence/presence; if the test suite didn't run, just chose any other job

@ev-br
Copy link
Member Author

ev-br commented Mar 16, 2024

Curious about the cross-flavor blas analysis. Is the intention to proceed in the 1.13rc phase, or is it realistically for 1.14 and beyond?

@ev-br ev-br force-pushed the test_bad_geneig branch from ac5641c to 85ba809 Compare March 18, 2024 07:36
@ev-br
Copy link
Member Author

ev-br commented Mar 18, 2024

@h-vetinari Any chance you'd be able to check that this PR fixes the failure seen on conda-forge? An alternative is of course to just land it and see what falls out in your next round of cross-blas analysis.

@h-vetinari
Copy link
Member

@h-vetinari Any chance you'd be able to check that this PR fixes the failure seen on conda-forge?

Gave it a shot in conda-forge/scipy-feedstock@c0a0f48. CI will run for 3-4 hours before we have results.

@h-vetinari
Copy link
Member

So it appear that this fixes the errors - I've started a second run to be sure. In any case thanks a lot for coming up with a fix! In this case I'd like to have this in 1.13, so I'm adding the backport-candidate label.

@h-vetinari h-vetinari added the backport-candidate This fix should be ported by a maintainer to previous SciPy versions. label Mar 18, 2024
@lucascolley lucascolley added this to the 1.13.0 milestone Mar 18, 2024
Copy link
Member

@h-vetinari h-vetinari left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I'm not familiar with the tests here, but I've tested this in conda-forge and would like to see this go in.

@ilayn
Copy link
Member

ilayn commented Mar 26, 2024

Always happy to remove the old crust anyways. Thanks @ev-br, @h-vetinari !

@ilayn ilayn merged commit 5a6af92 into scipy:main Mar 26, 2024
@tylerjereddy tylerjereddy removed the backport-candidate This fix should be ported by a maintainer to previous SciPy versions. label May 2, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

defect A clear bug or issue that prevents SciPy from being installed or used as expected scipy.linalg

Projects

None yet

Development

Successfully merging this pull request may close these issues.

BUG: osx-64 scipy 1.9.1 test_bad_geneig numerical error

5 participants