KrylovKit.jl
KrylovKit.jl copied to clipboard
Jutho
Krylov methods for linear problems, eigenvalues, singular values and matrix functions
KrylovKit does not precompile with Julia 1.2/1.3 and it is only tested with Julia >= 1.4.
Hi, I'm a big fan of the package. Being able to solve eigenvalue problems without vectors and matrices is awesome. Thanks you very much! I was wondering about the possibility...
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...
Almost working version of the eigsolve rrules. Seems to be failing only for real matrices, where comparison with rrule for eigen does not give equal results. Specifically, it seems that...
Using the `exponentiate` interface to `expintegrator`, I have been finding that the `tol` keyword does not behave as expected. Please see the minimal sample code to reproduce below. I was...
With v0.5, there is a method error if the initial vector is Abstract, e.g. ```julia using KrylovKit n = 10 KrylovKit.eigsolve(rand(n, n), 1:n, 1, :LM) # ERROR: MethodError: no method...
When using Arnoldi on some (badly scaled) problems, `trexc!` throws `LAPACKException(1)` with the following stacktrace ``` [1] trexc!(::Int64, ::Int64, ::SubArray{Float64,2,Array{Float64,2},Tuple{UnitRange{Int64},UnitRange{Int64}},false}, ::SubArray{Float64,2,Array{Float64,2},Tuple{UnitRange{Int64},UnitRange{Int64}},false}) at /Users/nrontsis/KrylovKit.jl/src/dense/linalg.jl:415 [2] permuteschur!(::SubArray{Float64,2,Array{Float64,2},Tuple{UnitRange{Int64},UnitRange{Int64}},false}, ::SubArray{Float64,2,Array{Float64,2},Tuple{UnitRange{Int64},UnitRange{Int64}},false}, ::Array{Int64,1}) at /Users/nrontsis/KrylovKit.jl/src/dense/linalg.jl:335 [3]...
Hi, I am still on KK0.4.2. Somehow, this used to work but not anymore ``` using KrylovKit using CuArrays # version 0.2.2 CuArrays.allowscalar(false) import LinearAlgebra: mul!, axpby! mul!(x::CuArray, y::CuArray, α::T)...
I suggest to add KrylovKit.jl to https://github.com/JuliaLinearAlgebra/
I am trying to compute eigenvalues of matrices with extended precision but the function `eigsolve` return an error. For example with the following matrix ```julia using LinearAlgebra, LinearMaps, Quadmath, KrylovKit...
