You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
I suspect mul! hits the generic fallback which would make it fascinatingly inefficient since getindex for Q has Θ(n²) complexity so the complete mul! computation is Θ(n⁵). Applying a compact Q matrix is really an in-place operation. We could hide that fact by adding a mul! method for Q matrices that just calls copy! and lmul!.
Following this Julia Discourse thread:
mul!
seems to be much slower thanlmul!
,rmul!
and*
in a denseQR
factorisation:Minimal example:
In the example above
mul!
runs in the order of seconds, while*
in the order of microseconds.Julia
version:1.0.1
macOS
10.14.2
The text was updated successfully, but these errors were encountered: