mikmoore

`\` for square singular matrices

4

From this [Discourse discussion](https://discourse.julialang.org/t/possible-bug-with-operator/84753). ```julia using LinearAlgebra julia> ones(3,2) \ ones(3) 2-element Vector{Float64}: 0.5 0.4999999999999999 julia> ones(3,3) \ ones(3) ERROR: SingularException(2) julia> ones(3,4) \ ones(3) 4-element Vector{Float64}: 0.2500000000000001 0.25 0.25...

With `ENV["EDITOR"] = "code"` (i.e., VS Code), `edit` does not work from a *standalone* WSL terminal. It does work from a terminal *embedded in VS Code*. I've chased the issue...

`round(::Type{<:AbstractFloat}, x, ::RoundingMode)` violates docstring

5

*I don't have the development branch checked out right now, so forgive me for manually-interpreting this code and let me know if any of my concerns are due solely to...

StaticArrays at `v1.6.2` on Julia `v1.9.1`. StaticArrays attempts to index empty arrays, resulting in an error, when the inner dimension of a `SMatrix`-`SMatrix` multiply is empty. I.e., for `M x...

Docstring pollution

3

It appears that some `Base` functions extended to `::Quaternion` arguments contribute their own docstrings that are near-verbatim repeats of the existing versions from `Base`. They don't provide any new information...

```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...

Re-associating matrix-chain multiplication yields unpredictable results

6

Version 1.9.4. ``` *(A, B::AbstractMatrix, C) A * B * C * D Chained multiplication of 3 or 4 matrices is done in the most efficient sequence, based on the...

`\` for square singular matrices

6

From this [Discourse discussion](https://discourse.julialang.org/t/possible-bug-with-operator/84753). ```julia using LinearAlgebra julia> ones(3,2) \ ones(3) 2-element Vector{Float64}: 0.5 0.4999999999999999 julia> ones(3,3) \ ones(3) ERROR: SingularException(2) julia> ones(3,4) \ ones(3) 4-element Vector{Float64}: 0.2500000000000001 0.25 0.25...

2x2 hermitian `eigen` method prone to under/overflow

2

```julia-repl julia> eigen(Hermitian(SA[2 1; 1 1] * 1f-30)) Eigen{Float32, Float32, SMatrix{2, 2, Float32, 4}, SVector{2, Float32}} values: 2-element SVector{2, Float32} with indices SOneTo(2): 3.8196598f-31 2.618034f-30 vectors: 2×2 SMatrix{2, 2, Float32,...

Julia v1.11.1 with StaticArrays v1.9.8. Discovered while attempting ```julia julia> using StaticArrays julia> M = @SMatrix [0 0 1f1; 0 0 0; 0 0 0] 3×3 SMatrix{3, 3, Float32, 9}...