Daniel Karrasch
Daniel Karrasch
Would be interesting to know what other languages do about it (Python, Matlab)...
This "has always been like this": ```julia _ _ _ _(_)_ | Documentation: https://docs.julialang.org (_) | (_) (_) | _ _ _| |_ __ _ | Type "?" for help,...
... and in `gemm_wrapper!`. Maybe those aliasing tests got lost a bit over the recent rearrangements?
I'm not sure how these functions would be used. Symmetry (or Hermitian-ness) is checked exactly, based on the components, and we don't provide an "almost-symmetric check", for good reasons, because...
Somebody said (perhaps in Slack or somewhere in the issues here), that MATLAB dropped its dependency on Arpack and has its own `eigs`.
This was also changed in JuliaLang/julia#46196 on purpose. `Q`s are not readily accessible in memory, their components need to be computed one by one, so printing them like a matrix...
Can this be closed?
I've looked into it. It's a pretty big deal, after all. We have some matrix printing code that does all computations regarding alignment of elements, fill-up with ellipses etc. That's...
> Is that more expensive than before? Of course. You'd allocate an entirely new matrix in memory. > Was the matrix previously allocated when printing Q? No, for big Q's,...
What if we change the output to ```julia julia> qr(randn(2,2)) LinearAlgebra.QRCompactWY{Float64, Matrix{Float64}, Matrix{Float64}} Q factor: 2×2 LinearAlgebra.QRCompactWYQ{Float64, Matrix{Float64}, Matrix{Float64}} R factor: 2×2 Matrix{Float64}: -2.3358 -0.951623 0.0 -0.608775 ``` ? Like...