I am seeing crashing/segfaults from Pardiso.jl with (project) Pardiso.
I'm not sure if this is related, but the version of Pardiso that I got with my academic license was labelled as 7.0, though I get output that includes headers with PARDISO 6.0.0 in them.
The crashing is stochastic, the following code copied from one of the examples seems to triggers it:
using Pardiso
using SparseArrays
using Random
using Printf
using Test
verbose = true
n = 4 # The number of equations.
m = 3 # The number of right-hand sides.
A = sparse([ 1. 0 -2 3
0 5 1 2
-2 1 4 -7
3 2 -7 5 ])
B = rand(n,m)
ps = PardisoSolver()
if verbose
set_msglvl!(ps, Pardiso.MESSAGE_LEVEL_ON)
end
set_matrixtype!(ps, Pardiso.REAL_SYM_INDEF)
pardisoinit(ps)
fix_iparm!(ps, :N)
A_pardiso = get_matrix(ps, A, :N)
set_phase!(ps, Pardiso.ANALYSIS)
set_perm!(ps, randperm(n))
pardiso(ps, A_pardiso, B)
@printf("The factors have %d nonzero entries.\n", get_iparm(ps, 18))
set_phase!(ps, Pardiso.NUM_FACT)
pardiso(ps, A_pardiso, B)
set_phase!(ps, Pardiso.SOLVE_ITERATIVE_REFINE)
julia> B
malloc(): invalid size (unsorted)
signal (6): Aborted
in expression starting at REPL[24]:1I'm using Julia 1.6.0.
I am seeing crashing/segfaults from Pardiso.jl with (project) Pardiso.
I'm not sure if this is related, but the version of Pardiso that I got with my academic license was labelled as 7.0, though I get output that includes headers with
PARDISO 6.0.0in them.The crashing is stochastic, the following code copied from one of the examples seems to triggers it:
I'm using Julia 1.6.0.