topic Re:pardiso error with segmentation fault in Intel® oneAPI Math Kernel Library

pardiso error with segmentation fault

NEROKAI — Tue, 09 Apr 2024 07:00:48 GMT

I tried to encapsulate PARDISO with a simple wrapper, why doesn't this kind of encapsulation work? Not only in C++, but also in C language, errors occur when using `pt` across functions. This is the minimal unit for reproducing the error. Thank you.

Re:pardiso error with segmentation fault

Gennady_F_Intel — Tue, 09 Apr 2024 12:24:23 GMT

how do you allocate solution array ( which is double *result into your code)?

Re:pardiso error with segmentation fault

Gennady_F_Intel — Wed, 10 Apr 2024 06:17:52 GMT

I hardcoded your code by allocating the solution array and solution has been passed:

void NumericalAlgebra::solve(double* b, double* result)

{

MKL_INT phase=33;

double* x = (double*)mkl_malloc(N*sizeof(double), 64);

pardiso(m_pt,&m_maxfct,&m_mnum,&m_type,&phase,&N,values, rowindex,colindex,&idum,&nrhs,m_iparm,&m_msglvl, b, /*result*/ x, &m_error);

if( m_error != 0 ) {

printf("... Failure. Pardiso Fails at the Solution ( phase == 33 ) stage. m_error == %lld \n", m_error );

} else {

printf("... Passed - Pardiso at the Solution ( phase == 12) stage \n");

}

here is the log output:

....

number of non-zeros in U: 1

number of non-zeros in L+U: 6

gflop for the numerical factorization: 0.000000

gflop/s for the numerical factorization: 0.000027

... Passed - Pardiso at the Solution ( phase == 12) stage

Re: Re:pardiso error with segmentation fault

NEROKAI — Thu, 11 Apr 2024 07:13:04 GMT

Thank you very much. I made a very stupid mistake. I thought pardiso could allocate memory automatically. Thank you!