-
-
Notifications
You must be signed in to change notification settings - Fork 35
Materialize adjoint in mul with HermOrSym
#1191
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
Codecov ReportAll modified and coverable lines are covered by tests ✅
Additional details and impacted files@@ Coverage Diff @@
## master #1191 +/- ##
=======================================
Coverage 91.89% 91.90%
=======================================
Files 34 34
Lines 15360 15366 +6
=======================================
+ Hits 14115 14122 +7
+ Misses 1245 1244 -1 ☔ View full report in Codecov by Sentry. |
|
x-ref JuliaLang/julia#52787 for the Bool case. |
|
Why not implement @mcabbott 's suggestion, though? It is faster for small matrices and avoids a copy. |
|
We seem to avoid wrapping results in A' * Q = copy((Q' * A)')We return wrapped results only when the first factor is an adjoint vector. |
|
Why, though? I don't see the problem with wrapping results in |
|
@andreasnoack Do you remember why we typically materialize products that we compute by the adjoint formula, except for |
|
I don't remember the specifics here. Only that the wrapped versions can sometimes lead to performance issues because they wrappers sometimes causes dispatch to slow fallback definitions for |
|
Another reason for going with this as is is that this currently returns a |
Fixes #850, at least the original issue. Promoting
Boolto someBlasFloatseems to be very tricky in terms method ambiguities.