-
Notifications
You must be signed in to change notification settings - Fork 135
Allow more specialized linalg.solve
assume_a
cases
#1273
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
Conversation
55d0b79
to
41c2131
Compare
b70fca9
to
8002c87
Compare
Codecov ReportAttention: Patch coverage is
Additional details and impacted files@@ Coverage Diff @@
## main #1273 +/- ##
==========================================
- Coverage 82.01% 82.00% -0.01%
==========================================
Files 188 188
Lines 48442 48478 +36
Branches 8659 8665 +6
==========================================
+ Hits 39730 39755 +25
- Misses 6566 6575 +9
- Partials 2146 2148 +2
🚀 New features to boost your workflow:
|
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.
Minor suggestions
transposed=False, | ||
check_finite=True, | ||
lower: bool = False, | ||
overwrite_a: bool = False, |
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.
I thought we decided not to expose this option to users?
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.
We don't. The docstrings specify it's ignored and PyTensor will perform inplace if possible. This way the signature is the same for the user.
71e0518
to
0519629
Compare
linalg.solve
assume_a
cases
Accompanying notebook: https://gist.github.com/ricardoV94/005fb4b48d274bcf383a3eb5456462d3
Also fixes #1271
We want to implement solve_banded and solveh_banded, which allows working directly with the sparse representation of the matrix which is much more practical. Will do that in a follow up PR as it requires more work.
Would also be nice to implement the tridiagonal / banded solvers for numba.
📚 Documentation preview 📚: https://pytensor--1273.org.readthedocs.build/en/1273/