LinearAlgebra.jl
LinearAlgebra.jl copied to clipboard
JuliaLang
Julia comes with OpenBLAS and that's our default libblastrampoline backend. https://github.com/OpenMathLib/OpenBLAS/pull/4577 (which was a continuation of https://github.com/OpenMathLib/OpenBLAS/pull/2255 by @stevengj), which should come in the not-yet-released OpenBLAS v0.3.28, introduced callbacks to...
These are PRs that: 1. Have the linear algebra label 2. Only modify files in `stdlib/LinearAlgebra` ``` PR #56611: Update svd.jl documentation PR #56506: use LQ for "wide A" \...
``` 1412 -> 2 1629 -> 3 2193 -> 4 2398 -> 5 2929 -> 6 2934 -> 7 2936 -> 8 2941 -> 9 3003 -> 10 3014 ->...
In the ` aarch64-apple-darwin` job https://buildkite.com/julialang/julia-master/builds/41270#0192a51c-3afb-4a6f-a53f-3c3b9f537ed4, there is ```julia The global RNG seed was 0xf3fab14bb37f52ec47cdcf1f85e908bd. Error in testset LinearAlgebra/lapack: Test Failed at /Users/julia/.julia/scratchspaces/a66863c6-20e8-4ff4-8a62-49f30b1f605e/agent-cache/default-honeycrisp-HL2F7YQ3XH.0/build/default-honeycrisp-HL2F7YQ3XH-0/julialang/julia-master/julia-d3cd3f4a3a/share/julia/stdlib/v1.12/LinearAlgebra/test/lapack.jl:441 Expression: ≈(triu(inv(A)), triu(LAPACK.sytri!('U', B, ipiv)), rtol =...
I want to run an example in the [docs](https://docs.julialang.org/en/v1/stdlib/LinearAlgebra/#LinearAlgebra.SymTridiagonal): ```julia julia> using LinearAlgebra julia> A = SymTridiagonal(fill([1 2; 3 4], 3), fill([1 2; 3 4], 2)) 3×3 SymTridiagonal{Matrix{Int64}, Vector{Matrix{Int64}}}: Error...
in 1.10+, matrix decomposition components aren't always AbstractMatrix and don't behave like them
Tried upgrading Julia version from 1.9 in one of my envs, and noticed this: ``` julia> using LinearAlgebra julia> f(a::AbstractMatrix) = sum(a) # MWE julia> f(qr(rand(5,5)).Q) # 1.9: -1.4572588703587335 #...
These should probably just be removed https://github.com/JuliaLang/julia/blob/84a2458e0504d92b3db32bb367e449377c802593/stdlib/LinearAlgebra/src/factorization.jl#L102-L104 `convert` should be defined quite sparingly and these are different things.
This might be too ambitious, but broadcasting a function that may conditionally preserve zeros fails currently: ```julia julia> (x -> rand() < 0.5 ? 0 : 1).(Diagonal([1,2])) ERROR: ArgumentError: cannot...
It could help to determine whether a symmetric or Hermitian sparse matrix wrapped in `Symmetric` or `Hermitian` contains both triangles or not ('F' = full). With sparse matrices on GPUs,...
This is quite problematic: ``` julia> Any[824.9999999999999] ≈ [825] true julia> Number[824.9999999999999] ≈ [825] false julia> Union{Number,Missing}[824.9999999999999] ≈ [825] true julia> Real[824.9999999999999] ≈ [825] false julia> Union{Real,Missing}[824.9999999999999] ≈ [825] true...
