Quantcast
Channel: Intel® oneAPI Math Kernel Library & Intel® Math Kernel Library
Viewing all articles
Browse latest Browse all 2652

pardiso transpose solve vs non-transpose solve

$
0
0

Hello,

I am wondering what differences to expect between loading in a CSR3 matrix and doing a solve, and loading in a CSC3 matrix and doing a transpose solve.  The ia and ja matrices are swapped appropriately in calling the solver.  This is for an unsymmetric matrix, mtype=11.

In essence the 2 matrices have identical entries, except one is transposed.  I am seeing numerical differences in the convergence.  Using CSC3, all of my examples converge.  Using CSR3, one fails.

Is it reasonable to expect the same results.  This would then lead to issues with my CSR3 implementation.

 


Viewing all articles
Browse latest Browse all 2652

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>