-
-
Notifications
You must be signed in to change notification settings - Fork 73
Fix handling for default algorithms for DualLinearProblems #775
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
base: main
Are you sure you want to change the base?
Conversation
function SciMLBase.init(prob::DualAbstractLinearProblem, alg::DefaultLinearSolver, args...; kwargs...) | ||
if alg.alg === DefaultAlgorithmChoice.GenericLUFactorization | ||
return __init(prob, alg, args...; kwargs...) | ||
else | ||
return __dual_init(prob, alg, args...; kwargs...) | ||
end |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This needs a JET test
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Tests added. Looks good on 1.11 but fails on LTS
tests fail |
CI is having some issues here, some are complaining about Pkg not being available, and one says that Julia isn't a recognized executable. |
NoPre is a real failure |
So the normal default algorithm But the default algorithm |
Checklist
contributor guidelines, in particular the SciML Style Guide and
COLPRAC.
Additional context
This makes it so that the default alg is chosen based on the primal A and b before going to make the Dual cache. This also means that for small problems GenericLU will be used, and the cache will be a normal LinearCache.