search

JuliaSparse / SuiteSparseGraphBLAS.jl

Sparse, General Linear Algebra for Graphs!
MIT License
102 stars 16 forks source link

UndefVarError in Base.copyto! #124

Open atbug opened 9 months ago

atbug commented 9 months ago 
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]