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