2025.1-PARDISO: Sig11 in numerical factorization

[rbunger]

Hi,

unfortunately I do not have a MWE showing the problem: In some cases I see a segmentation fault in the 2025.1-PARDISO numerical factorization. It does not always crash. I do not see the problem with the 2024.2.1-PARDISO. The matrix type is complex unsymmetric. I know, the description is not specific enough...

Thanks,

Rainer

[Fengrui]

Hi Rainer, 

It will also be helpful if you could share the matrix that was used in the segmentation fault case.

Thanks,

Fengrui

[rbunger]

Hi Fengrui,

thanks for your interest. Up to now I can say that it happens with very large matrices so I think I can't upload them. Also the problem is not reproducible. It does not depend on the compiler I use for the built (GCC or icpx). I have the problem under observation, I will post an update when I can narrow it down a little bit further...

Thanks

Rainer

[c_sim]

Hi Rainer, 

Thank you for posting.

Are you working with LP64 or ILP64 interface? If you are using LP64, could you check the ILP64 interface? Because in 2025.1 version there was a slight change that might increase the non-zeros of LU factors slightly. Since your matrices are large, there might be an integer overflow with LP64.

Kind Regards,

Chris

[rbunger]

Hi Chris,

I'm using the ILP64 interface.

Rainer

[c_sim]

Thank you for the prompt response. In that case, please let us know when you can narrow down the problem a bit more. Ideally, it would be helpful to know the matrix and iparm settings.

[rbunger]

Hi,

I have updated my original posting:

community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/2025-1-PARDISO-Sig11-in-numerical-factorization/m-p/1682011#M37056 

The error has just occurred with a complex symmetric matrix in phase=11, iparm[1] = 3, iparm[34] = 1.

Thanks

Rainer

[rbunger]

Hi,

I don't know if I just mixed it up the first time when I reported the bug or if the Sig11 just occurs at different stages: I now see the Sig11 with a complex symmetric matrix in the "Analysis" step ("phase" = 11). I also have a complex unsymmetrical matrix to solve but the error just occurred while processing the symmetric one.  Also, I call PARDISO_64 instead of PARDISO.

Thanks

Rainer

"iparms" & "phase", apart from the default settings, are as follows (C-indexing):

iparm[1] = 3;

iparm[34] = 1;

MKL_INT64 phase = 11;