yakovbraver

Occasional `LAPACKException(22)` in `eigsolve`

4

Consider diagonalisation of a real symmetric matrix in the following MWE: ```julia import BandedMatrices as bm import KrylovKit as kk import LinearAlgebra as la function test_krylov(a, b, N, n_j) h...

Performance regression in non-autonomous Linear ODE solvers

4

There are two performance issues when using state-independent solvers for non-autonomous linear ODEs: 1. A ~2x performance regression going from OrdinaryDiffEq v6.54.0 to v6.55.0 when the linear operator is either...

Consider summing a product of three `besselj(n, x)`s when the type of `n` is `Int64` (the default `Int` on a 64-bit machine) and type of `x` is `Float32`: ```julia import...

Consider calling `besselj(n, x)` where `x` is negative, while `n` is a large negative integer: ```julia julia> Bessels.besselj(-172, -2.0) NaN julia> Bessels.besselj(-37, -2.0f0) NaN32 ``` When at least one argument...

The [docs](https://julialinearalgebra.github.io/ArnoldiMethod.jl/stable/#Pre-allocating-and-custom-matrix-types) provide the following example: ```julia # allocates workspace for 20-dimensional Krylov subspace arnoldi = ArnoldiWorkspace(100, 20) ``` This gives a MethodError because the relevant constructor is `ArnoldiWorkspace(::Type{T}, ::Int64,...

`eigen` returns `ComplexF64` eigenvalues for a Hermitian matrix

6

Consider calling `eigen` for a `ComplexF64` matrix which happens to be real-diagonal: ```julia julia> D = ComplexF64[1 0; 0 2]; julia> eigen(D) Eigen{ComplexF64, ComplexF64, Matrix{ComplexF64}, Vector{ComplexF64}} values: 2-element Vector{ComplexF64}: 1.0...

- [ ] what version of Makie are you running? **v0.22.4** - [ ] can you reproduce the bug with a fresh environment ? **Yes** - [ ] What platform...