Use pivoting as the default in LU regardless of the element type. #32989
Conversation
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
fredrikekre
approved these changes
Aug 20, 2019
For types that weren't subtypes of AbstractFloat, we used to try to LU factorize without pivoting and only use pivoting when it failed. This caused large numerical errors when computing the LU for element types which promoted to float like numbers such as most integers. The behavior was never documented and is error prone. Hence, this PR removes the behavior.
|
I think it would be good if we could try backporting this to 1.3 since many people have complained about JuliaLang/LinearAlgebra.jl#605. If PkgEval shows that it breaks too many packages (I don't expect that) then we can drop it. |
fredrikekre
added a commit
that referenced
this pull request
Aug 23, 2019
fredrikekre
added a commit
that referenced
this pull request
Aug 24, 2019
KristofferC
pushed a commit
that referenced
this pull request
Aug 25, 2019
…2989) For types that weren't subtypes of AbstractFloat, we used to try to LU factorize without pivoting and only use pivoting when it failed. This caused large numerical errors when computing the LU for element types which promoted to float like numbers such as most integers. The behavior was never documented and is error prone. Hence, this PR removes the behavior. (cherry picked from commit 5af3c2a)
KristofferC
pushed a commit
that referenced
this pull request
Aug 25, 2019
36 tasks
|
Note that this was merged even with a (correctly) red CI. |
This was referenced Sep 5, 2019
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
For types that weren't subtypes of AbstractFloat, we used to try
to LU factorize without pivoting and only use pivoting when it failed.
This caused large numerical errors when computing the LU for element
types which promoted to float like numbers such as most integers.
The behavior was never documented and is error prone. Hence, this
PR removes the behavior.
Fixes JuliaLang/LinearAlgebra.jl#605. Alternative to #30953