-
Notifications
You must be signed in to change notification settings - Fork 19
Open
Description
using SparseArrays, SuiteSparseGraphBLAS, KrylovKit
A = sprand(Float64, 100, 100, 0.05)
M = GBMatrix(A' + A)
eigsolve(M, size(M, 1), 3, :SR, ishermitian=true, maxiter=1000)
The error is
ERROR: UndefVarError: `f` not defined
Stacktrace:
[1] copyto!
@ ~/.julia/packages/SuiteSparseGraphBLAS/x2QAm/src/operations/broadcasts.jl:332 [inlined]
[2] materialize!
@ ./broadcast.jl:884 [inlined]
[3] materialize!
@ ./broadcast.jl:881 [inlined]
[4] mul!
@ /mmfs1/gscratch/cmt/wangc/.julia/juliaup/julia-1.9.3+0.x64.linux.gnu/share/julia/stdlib/v1.9/LinearAlgebra/src/generic.jl:142 [inlined]
[5] mul!
@ /mmfs1/gscratch/cmt/wangc/.julia/juliaup/julia-1.9.3+0.x64.linux.gnu/share/julia/stdlib/v1.9/LinearAlgebra/src/matmul.jl:276 [inlined]
[6] basistransform!(b::KrylovKit.OrthonormalBasis{GBVector{Float64, Float64}}, U::SubArray{Float64, 2, Matrix{Float64}, Tuple{UnitRange{Int64}, UnitRange{Int64}}, false})
@ KrylovKit ~/.julia/packages/KrylovKit/diNbc/src/orthonormal.jl:325
[7] eigsolve(A::GBMatrix{Float64, Float64}, x₀::Vector{Float64}, howmany::Int64, which::Symbol, alg::Lanczos{ModifiedGramSchmidt2, Float64})
@ KrylovKit ~/.julia/packages/KrylovKit/diNbc/src/eigsolve/lanczos.jl:107
[8] #eigsolve#38
@ ~/.julia/packages/KrylovKit/diNbc/src/eigsolve/eigsolve.jl:202 [inlined]
[9] eigsolve
@ ~/.julia/packages/KrylovKit/diNbc/src/eigsolve/eigsolve.jl:180 [inlined]
[10] #eigsolve#37
@ ~/.julia/packages/KrylovKit/diNbc/src/eigsolve/eigsolve.jl:178 [inlined]
[11] eigsolve
@ ~/.julia/packages/KrylovKit/diNbc/src/eigsolve/eigsolve.jl:178 [inlined]
Metadata
Metadata
Assignees
Labels
No labels