Intel® oneAPI Math Kernel Library
Ask questions and share information with other developers who use Intel® Math Kernel Library.
Announcements
This community is designed for sharing of public information. Please do not share Intel or third-party confidential information here.
6590 Discussions

MKL Pardiso (version 11.2.3): wrong output of phase 331 with multiple rhs and Schur complement enabled

zampini__stefano
Beginner
124 Views

Hi,

I recently started using MKL_PARDISO. I noticed that phase 331 gives the wrong result if you want to solve for multiple right hand sides using the Schur complement feature.

Attached a code to reproduce the problem. I just copied the example you provide with mkl distribution for the Schur complement and added multiple rhs.

I'm using composer_xe_2015.3.187, with MKL 11.2.3

Stefano

0 Kudos
2 Replies
zampini__stefano
Beginner
124 Views

Just a clarification about the issue.

Since the 3 rhs I'm using are all equal to one, I'm expecting that the output of phase 331 will be the same for x[0...7],x[8...15] and x[16...23]. Instead, the numbers in x[8...15] and x[16...23] (that are the same) are not in the same order as in x[0..7].

asd__asdqwe
Beginner
124 Views

Hello Stefano,

Here are some problems in your piece of code :

  1. You need to set iparm(36) to 2 instead of 1
  2. If you want the same results from phase = 33, you need to call phase = 331, then 332 (or dsytrs as done in the MKL example), then 333.

Here is a program that returns the same results.

Reply