Intel® oneAPI Math Kernel Library
Ask questions and share information with other developers who use Intel® Math Kernel Library.
Announcements
FPGA community forums and blogs on community.intel.com are migrating to the new Altera Community and are read-only. For urgent support needs during this transition, please visit the FPGA Design Resources page or contact an Altera Authorized Distributor.
7234 Discussions

Changed behavior in Pardiso, MKL 11.3.3 vs MKL 11.1.3

Jens_E_
New Contributor I
424 Views

Hi,

I recently upgraded from Compser XE 2013SP1U3 to Composer XE 2016U3, and thus from MKL 11.1.3 to MKL 11.3.3. This has caused a slightly changed behavior in the output from Pardiso in the following situations:

  1. When trying to "solve" a system of equations of dimension 0 (I have a set of unit tests which for one test degenerates to this case).
    1. 11.1.3: the error code is 0 (no error), so that the solution may be interpreted as the empty vector.
    2. 11.3.3: the error code is -1, indicating "inconsistent input."
  2. When trying to solve a system Ax=b in which A is an 5x5 matrix of all zeros, and b is a 5x1 matrix of all zeros. 
    1. 11.1.3: the error code is -7
    2. 11.3.3: the error code is -4

Are these changes expected?

Of course, both cases are degenerate and not of any practical interest. I just need to confirm that this changed behavior does not suggest any errors in my (other) input to Pardiso.

Thanks!

Jens 

0 Kudos
1 Reply
Alexander_K_Intel2
424 Views

Hi Jens,

You are correct, we implement input parameters check and removed diagonal matrices branch so mentioned behaviour expected

Thanks,

Alex 

0 Kudos
Reply