Crash of eigenvalues calculations with arpack

New to Julia
, , ,
1

I have a problem of crash of eigenvalues calculations by Arpack.v0.5.4 only if the size of the array becomes very large

The bugging line of the program is
(line 99)
‘’‘Eigen = eigs(Sierp3, nev=1, which=:LM, sigma=1.7)[1]
‘’’

The error message is:
ERROR: LoadError: ArgumentError: new length must be ≥ 0
Stacktrace:
[1] throw_argerror(s::String)
@ Base ./array.jl:317
[2] resize!(a::Vector{Int64}, nl::Int64)
@ Base ./array.jl:1454
[3] _allocres
@ ~/.julia/juliaup/julia-1.11.6+0.x64.linux.gnu/share/julia/stdlib/v1.11/SparseArrays/src/higherorderfns.jl:254 [inlined]
[4] _noshapecheck_map(f::typeof(-), A::SparseMatrixCSC{ComplexF64, Int64}, Bs::SparseMatrixCSC{ComplexF64, Int64})
@ SparseArrays.HigherOrderFns ~/.julia/juliaup/julia-1.11.6+0.x64.linux.gnu/share/julia/stdlib/v1.11/SparseArrays/src/higherorderfns.jl:183
[5] map
@ ~/.julia/juliaup/julia-1.11.6+0.x64.linux.gnu/share/julia/stdlib/v1.11/SparseArrays/src/higherorderfns.jl:1187 [inlined]
[6] -
@ ~/.julia/juliaup/julia-1.11.6+0.x64.linux.gnu/share/julia/stdlib/v1.11/SparseArrays/src/sparsematrix.jl:2253 [inlined]
[7] -
@ ~/.julia/juliaup/julia-1.11.6+0.x64.linux.gnu/share/julia/stdlib/v1.11/SparseArrays/src/sparsematrix.jl:4341 [inlined]
[8] _eigs(A::SparseMatrixCSC{ComplexF64, Int64}, B::UniformScaling{Bool}; nev::Int64, ncv::Int64, which::Symbol, tol::Float64, maxiter::Int64, sigma::Float64,
v0::Vector{ComplexF64}, ritzvec::Bool, explicittransform::Symbol, check::Int64)
@ Arpack ~/.julia/packages/Arpack/FCvNd/src/Arpack.jl:208
[9] _eigs
@ ~/.julia/packages/Arpack/FCvNd/src/Arpack.jl:75 [inlined]
[10] eigs
@ ~/.julia/packages/Arpack/FCvNd/src/Arpack.jl:51 [inlined]
[11] eigs#3
@ ~/.julia/packages/Arpack/FCvNd/src/Arpack.jl:50 [inlined]
[12] top-level scope
@ ~/essai3.jl:99
in expression starting at /home/michel/essai3.jl:99

The program works for a dimension of the table Sierp3 up to 344373768, it fails for a dimension of
3099363912 There is no intermadiate values between these two because it corresponds to special structures of matrices generated elsewhere in the program

Is anybody can help me ?
Michel

2

Welcome to Julia Discourse. There are several experts here who may be able to help you. Responses may be a bit slow for the next few days since many of them are attending JuliaCon 2025 right now.

I note that you’ve provided the error message and the offending line from your code, but you haven’t supplied enough information for others to reproduce your error. If you can do this it will be very helpful. Please have a look at Please read: make it easier to help you, including the part about how to format your code for presentation on Discourse , and especially the discussion of how to create a MWE (minimum working example). Posting one will maximize your chances of receiving useful assistance. Oftentimes, just creating such an example will help one to diagnose the problem oneself.
Edit: I’m going to retract my comments about a MWE in this case. Upon reflection, it would be difficult to provide one for a matrix that large. The code (and error message) formatting discussed in the link I provided would make your post easier to read, though. You can edit the post if you like, to add the triple backtick formatting.