Hi! Could you set msglvl=1 before reordering stage and provide us with output info? This info will be useful for us in addition to input matrix and PARDISO parameters.
[cpp]for (int i = 0; i < 64; i++) pt = NULL; for (int i = 0; i < 64; i++) iparm = 0; pardisoinit (pt, &mtype, iparm); nrhs = 1; // Number of right hand sides. mnum = 1; // Which factorization to use. maxfct = 1; msglvl = 1; // Print statistical information in file iparm = 1; // 1: No solver default 0: default[/cpp]
Thank you very much forthe attached reproducer! We'll investigate the issue and come back to you as the first resultsare available.