LinearAlgebra.jl
LinearAlgebra.jl copied to clipboard
JuliaLang
PkgEval log: https://s3.amazonaws.com/julialang-reports/nanosoldier/pkgeval/by_hash/7825364_vs_d63aded/PDMats.primary.log ``` * kron Kronecker product: Error During Test at /home/pkgeval/.julia/packages/PDMats/iX2rk/test/kron.jl:14 Got exception outside of a @test MethodError: Cannot `convert` an object of type Matrix{Float32} to an object...
Hi, The problem is as in the title. In short, `isposdef` method for `Eigen` factorization is incorrect. It does not check for hermiticity/symmetricity and only looks at the eigenvalues, which...
Is the following line in `lapack.jl` even reachable? https://github.com/JuliaLang/LinearAlgebra.jl/blob/2a1696a93683e947b46991ea454c4ffb1988e794/src/lapack.jl#L4141 The call to `chklapackerror(info[])` on line 4138 https://github.com/JuliaLang/LinearAlgebra.jl/blob/2a1696a93683e947b46991ea454c4ffb1988e794/src/lapack.jl#L4138 should throw an error if any of the eigenvectors failed to converge and...
I ran into an interesting numerical issue in Julia that seems to be related to a bug in LinearAlgebra. Specifically, I have a vector with norm $\approx 1$ that behaves...
Hello! I think that the symbol ⊗ is quite a common notation for denoting tensor products and the Kronecker product of matrices in linear algebra. I always use `⊗ =...
In Julia version 1.11.2, it appears as though I can speed up calculation of `opnorm` by taking advantage of the identity $$ \left\lVert A^\top A \right\rVert_{op} = \left\lVert A A^\top...
Consider a general and not necessarily normal matrix `A`. In MATLAB, one can do `[V,D,W] = eig(A)` where each column of `V` is a right eigenvector and each row of...
Coming from: https://github.com/JuliaLang/julia/pull/55830 Making this an issue rather than PR because there seems to be some question about what the solution should actually be, and it might be broader than...
The [JuliaLang/julia](https://github.com/JuliaLang/julia) repo has a `NEWS.md` which also documented changes in `LinearAlgebra`. We might need one for this repo
Running this piece of code, I get ```julia using LinearAlgebra, SparseArrays, BenchmarkTools function _spmatmul!(C, A, B, α, β) size(A, 2) == size(B, 1) || throw(DimensionMismatch()) size(A, 1) == size(C, 1)...
