Hello,
I have an issue when using Parallel Direct Sparse Solver for Clusters. I modified a test case sample (cl_solver_sym_sp_0_based_c.c) in "mkl/example/examples_cluster_c.tgz" routines which includes a system of linear equation (AX=b) with a sparse symmetric matrix (A).
When I use the original file and run the code with different number of cores the code PASSES the test. But, when I modify just a member of A matrix (REPLACE THE DIAGONAL ELEMENT IN 2nd ROW from -4.0 to 0.0), the code fails when I use multiple cores. It works using 1 core. (I attached the modified file to this post)
The original shape of matrix is:
float a[18] = { 7.0, /*0*/ 1.0, /*0*/ /*0*/ 2.0, 7.0, /*0*/
-4.0, 8.0, /*0*/ 2.0, /*0*/ /*0*/ /*0*/
1.0, /*0*/ /*0*/ /*0*/ /*0*/ 5.0,
7.0, /*0*/ /*0*/ 9.0, /*0*/
5.0, 1.0, 5.0, /*0*/
-1.0, /*0*/ 5.0,
11.0, /*0*/
5.0
The modified
float a[18] = { 7.0, /*0*/ 1.0, /*0*/ /*0*/ 2.0, 7.0, /*0*/
-0.0, 8.0, /*0*/ 2.0, /*0*/ /*0*/ /*0*/
1.0, /*0*/ /*0*/ /*0*/ /*0*/ 5.0,
7.0, /*0*/ /*0*/ 9.0, /*0*/
5.0, 1.0, 5.0, /*0*/
-1.0, /*0*/ 5.0,
11.0, /*0*/
5.0
I am using these routines for couples of days and I may miss something.
Thank you for your help.
Nima Mansouri