LinearAlgebra.jl
LinearAlgebra.jl copied to clipboard
JuliaLang
```julia-repl julia> VERSION v"1.10.3" julia> [1 1e-10; 0 1] ^ 1e20 ERROR: InexactError: Int64(1.0e20) ``` The issue is that there are multiple places where our exponentiation routines attempt to convert...
In `LinearAlgebra`, often there are shortcuts taken depending on whether a matrix has a particular structure (e.g. if it is `Diagonal`), which lets us evaluate certain functions by acting only...
`propertynames(F::SVD, true)` returns same set of elements as `propertynames(F::SVD)` https://github.com/JuliaLang/julia/blob/0e6d7972c62eb1dedaf00ea9b966a0060e2effb8/stdlib/LinearAlgebra/src/svd.jl#L209-L210
I haven't been able to replicate locally so can't bisect. But the regression is persistent on NanoSoldier so likely real. The bad commit should be one of these ones, https://github.com/JuliaLang/julia/compare/16d70b53f3e2493310b1bca61b3dcc861284a248...b8e714dfa3d81edd56ccda853322b6e42d94cf60....
As per the documentation for `mul!`, https://github.com/JuliaLang/julia/blob/263928f9ad450a28d601c3f185040987e0bc45d5/stdlib/LinearAlgebra/src/matmul.jl#L266-L268 The mwe ``` a = rand(3, 3); v = rand(3); mul!(a, a, a) ``` is 'protected' against wrong inputs for matrix-matrix product and...
MWE: ```julia using LinearAlgebra m = Symmetric(fill(ones(2,2), 2, 2)) eltype(m) n = Hermitian(fill(ones(2,2), 2, 2)) eltype(n) ``` On the other hand, if we do this with `UpperTriangular` or `LowerTriangular` the...
On the same line of [this issue](https://github.com/JuliaLang/LinearAlgebra.jl/issues/1052), it might be nice to have `nullspace` and `cond` functions defined for `QRPivoted` matrices.
```julia julia> D1 = Diagonal([1]); D2 = Diagonal([1,2]); julia> D1 .+ D2 2×2 Diagonal{Int64, Vector{Int64}}: 2 ⋅ ⋅ 3 julia> Matrix(D1) .+ D2 2×2 Matrix{Int64}: 2 1 1 3 ```...
Reference thread https://discourse.julialang.org/t/question-regarding-numerical-stability-among-julia-versions/112741 sqrt(A) is failing when A is non-symmetric with repeated eigenvalues.
I think it would be better to throw a more informative error message when `svd` is called with a matrix that contains `NaN` entries. ```julia julia> using LinearAlgebra julia> svd([0....
