SparseQR v.2019 inital release problem

[adlam]

I'am trying to execute the SparseQR methods, but get an unhandled exception when staring the either mkl_sparse_qr_reorder() or mkl_sparse_d_qr()

In both cases I construct the sparse_matrix_t object with one based indexing (SPARSE_INDEX_BASE_ONE). In order to verify the matrix I execute the matvec product (A*x) using sprase blas function mkl_sparse_d_mv() and in works properly (verify the result). Then when try to proceed with SparseQR routines I get error:

"Unhandled exception at 0x00007FFDA65867DF (mkl_intel_thread.dll) in ....,  Access violation reading location 0xFFFFFFFFFFFFFFFF."

The same error goes with zero-based matrices. 

In my program I use a lot of MKL functions, including sparse BLAS and Pardiso - everything else works fine. 

Any ideas what can be wrong? 

[MariaZh]

Hi,

Can you please provide a small reproducer for this case, so we will be able to check it on our side?

Best regards,
Maria

[adlam]

Here you have example dataset:

// Matrix data: assuming one-based indexing

nRows = 40;
nCols = 24;
nznum = 360; // number of non-zero elements;

pointerB = { 10, 19, 28, 37, 46, 55, 64, 73, 82, 91, 100, 109, 118, 127, 136, 145, 154, 163, 172, 181, 190, 199, 208, 217, 226, 235, 244, 253, 262, 271, 280, 289, 298, 307, 316, 325, 334, 343, 352 };

pointerE = { 19, 28, 37, 46, 55, 64, 73, 82, 91, 100, 109, 118, 127, 136, 145, 154, 163, 172, 181, 190, 199, 208, 217, 226, 235, 244, 253, 262, 271, 280, 289, 298, 307, 316, 325, 334, 343, 352, 361 };

columns = { 2, 3, 4, 5, 21, 22, 23, 24, 1, 2, 3, 4, 5, 21, 22, 23, 24, 1, 2, 3, 4, 5, 21, 22, 23, 24, 1, 2, 3, 4, 5, 21, 22, 23, 24, 1, 2, 3, 4, 5, 21, 22, 23, 24, 1, 2, 3, 4, 5, 21, 22, 23, 24, 1, 2, 3, 4, 5, 21, 22, 23, 24, 1, 2, 3, 4, 5, 21, 22, 23, 24, 1, 2, 3, 4, 5, 21, 22, 23, 24, 1, 2, 3, 4, 5, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24 };

values = { 0.007801517859, 0.000888034535, 0.035617131505, 1.000000000000, -0.876518082605, 1.282062390000, -0.019696949887, -0.019483465963, 0.008891192058, -0.109906056502, 0.001545121088, 0.045101165387, 1.000000000000, 0.124340980968, 0.010989533912, -0.033401421519, 0.000587524190, 0.002254107738, -0.051972579269, 0.003425019422, 0.064334367916, 1.000000000000, 0.042358852907, -0.039536523557, -0.030121979059, 0.052692523271, 0.001016243192, -0.032763315112, 0.013471489319, 0.121670743120, 1.000000000000, -0.016766350272, -0.030586529961, 0.055349274383, 0.107017826699, 0.000580904145, -0.023293925809, 1.020249499141, 0.005100992446, 1.000000000000, -0.014943391720, -0.011267140561, 0.481355569107, -0.433325394230, -0.007801517859, -1.560303571743, 0.023279990067, -0.000813437866, 0.000000000000, 1.290955777065, 0.889338706505, -0.012440921639, 0.029826764466, -0.124540671940, -0.008822560758, 0.033377431407, -0.001477758196, 0.000000000000, 0.011346925356, -0.109710569126, 0.000886209308, 0.045121543677, -0.065751959359, -0.002193263722, 0.053165934599, -0.003363886945, 0.000000000000, -0.050340130519, -0.033805660857, 0.043943301386, 0.037063791237, -0.045781897791, -0.000961932540, 0.111007386583, -0.013415761224, 0.000000000000, -0.042613359144, 0.011782241875, 0.097162455357, -0.059421326274, -0.035631325462, -0.000532128379, -0.005100992446, -1.020198489217, 0.000000000000, -0.017135127234, 0.009697009437, -0.438182282460, -0.485688823049, 1.560381586921, 0.007801517859, 0.000888034535, 0.035617131505, 1.000000000000, 0.000737568666, 0.000506329867, -0.000007081363, 0.000017040360, 0.008891192058, -0.109906056502, 0.001545121088, 0.045101165387, 1.000000000000, 0.000000065466, -0.000000633833, 0.000000005144, 0.000000260678, 0.002254107738, -0.051972579269, 0.003425019422, 0.064334367916, 1.000000000000, -0.000000768036, -0.000000515588, 0.000000670422, 0.000000565254, 0.001016243192, -0.032763315112, 0.013471489319, 0.121670743120, 1.000000000000, -0.000072264491, 0.000020055686, 0.000164733289, -0.000101045958, 0.000580904145, -0.023293925809, 1.020249499141, 0.005100992446, 1.000000000000, 0.020769765333, -0.011636781137, 0.526000389988, 0.588769249632, -0.007801517859, -1.560303571743, 0.023279990067, -0.000813437866, 0.000000000000, 0.000499004813, -0.000732505367, 0.000011252683, 0.000011090937, -0.124540671940, -0.008822560758, 0.033377431407, -0.001477758196, 0.000000000000, -0.000000718357, -0.000000063411, 0.000000192968, -0.000000003426, -0.065751959359, -0.002193263722, 0.053165934599, -0.003363886945, 0.000000000000, -0.000000646035, 0.000000603208, 0.000000459382, -0.000000803907, -0.045781897791, -0.000961932540, 0.111007386583, -0.013415761224, 0.000000000000, 0.000028537670, 0.000051869719, -0.000094116054, -0.000181445788, -0.035631325462, -0.000532128379, -0.005100992446, -1.020198489217, 0.000000000000, -0.017933993153, -0.013656317825, 0.583568122657, -0.520112697492, 1.560381586921, 0.007801517859, 0.000888034535, 0.035617131505, 1.000000000000, 0.000737568666, 0.000506329867, -0.000007081363, 0.000017040360, 0.008891192058, -0.109906056502, 0.001545121088, 0.045101165387, 1.000000000000, 0.000000065466, -0.000000633833, 0.000000005144, 0.000000260678, 0.002254107738, -0.051972579269, 0.003425019422, 0.064334367916, 1.000000000000, -0.000000768036, -0.000000515588, 0.000000670422, 0.000000565254, 0.001016243192, -0.032763315112, 0.013471489319, 0.121670743120, 1.000000000000, -0.000072264491, 0.000020055686, 0.000164733289, -0.000101045958, 0.000580904145, -0.023293925809, 1.020249499141, 0.005100992446, 1.000000000000, 0.020769765333, -0.011636781137, 0.526000389988, 0.588769249632, -0.007801517859, -1.560303571743, 0.023279990067, -0.000813437866, 0.000000000000, 0.000499004813, -0.000732505367, 0.000011252683, 0.000011090937, -0.124540671940, -0.008822560758, 0.033377431407, -0.001477758196, 0.000000000000, -0.000000718357, -0.000000063411, 0.000000192968, -0.000000003426, -0.065751959359, -0.002193263722, 0.053165934599, -0.003363886945, 0.000000000000, -0.000000646035, 0.000000603208, 0.000000459382, -0.000000803907, -0.045781897791, -0.000961932540, 0.111007386583, -0.013415761224, 0.000000000000, 0.000028537670, 0.000051869719, -0.000094116054, -0.000181445788, -0.035631325462, -0.000532128379, -0.005100992446, -1.020198489217, 0.000000000000, -0.017933993153, -0.013656317825, 0.583568122657, -0.520112697492, 1.560381586921, 0.007801517859, 0.000888034535, 0.035617131505, 1.000000000000, -0.872241323662, 1.284998321104, -0.019738010855, -0.019384658199, 0.008891192058, -0.109906056502, 0.001545121088, 0.045101165387, 1.000000000000, 0.124343772841, 0.010962503169, -0.033401202162, 0.000598641185, 0.002254107738, -0.051972579269, 0.003425019422, 0.064334367916, 1.000000000000, 0.042341057997, -0.039548469393, -0.030106445815, 0.052705619842, 0.001016243192, -0.032763315112, 0.013471489319, 0.121670743120, 1.000000000000, -0.016873473173, -0.030556799962, 0.055593470495, 0.106868038940, 0.000580904145, -0.023293925809, 1.020249499141, 0.005100992446, 1.000000000000, -0.014774064958, -0.011362010120, 0.485643818954, -0.428525417475, -0.007801517859, -1.560303571743, 0.023279990067, -0.000813437866, 0.000000000000, 1.293849234172, 0.885091306873, -0.012375673459, 0.029891074770, -0.124540671940, -0.008822560758, 0.033377431407, -0.001477758196, 0.000000000000, 0.011316289972, -0.109713273390, 0.000894438699, 0.045121397554, -0.065751959359, -0.002193263722, 0.053165934599, -0.003363886945, 0.000000000000, -0.050355098728, -0.033791684919, 0.043953944983, 0.037045165226, -0.045781897791, -0.000961932540, 0.111007386583, -0.013415761224, 0.000000000000, -0.042571055685, 0.011859132124, 0.097022940298, -0.059690296537, -0.035631325462, -0.000532128379, -0.005100992446, -1.020198489217, 0.000000000000, -0.017281335193, 0.009585675483, -0.433424708206, -0.489929073129 };

[MariaZh]

Hi,

I see in you reproducer, that the nRows is equal to 40, but the number of values in pointerB/pointerE arrays is 39.
If you're using 4arrays CSR representation, the size of pointerB/pointerE should be equal to nRows.

Please, refer to https://software.intel.com/en-us/mkl-developer-reference-c-sparse-blas-csr-matrix-storage-format:
"The length of the values and columns arrays is equal to the number of non-zero elements in A.The length of the pointerB and pointerE arrays is equal to the number of rows in A."

Best regards,
Maria

[adlam]

Hi,

It looks that I did not exported properly the first elements of my arrays shown in previous post. Here you have updated matrix information - the problem still remains:

nRows = 40;
nCols = 24;
nznum = 360; // non-zero element;

pointerB = { 1, 10, 19, 28, 37, 46, 55, 64, 73, 82, 91, 100, 109, 118, 127, 136, 145, 154, 163, 172, 181, 190, 199, 208, 217, 226, 235, 244, 253, 262, 271, 280, 289, 298, 307, 316, 325, 334, 343, 352 };

pointerE = { 10, 19, 28, 37, 46, 55, 64, 73, 82, 91, 100, 109, 118, 127, 136, 145, 154, 163, 172, 181, 190, 199, 208, 217, 226, 235, 244, 253, 262, 271, 280, 289, 298, 307, 316, 325, 334, 343, 352, 361 };

columns = { 1, 2, 3, 4, 5, 21, 22, 23, 24, 1, 2, 3, 4, 5, 21, 22, 23, 24, 1, 2, 3, 4, 5, 21, 22, 23, 24, 1, 2, 3, 4, 5, 21, 22, 23, 24, 1, 2, 3, 4, 5, 21, 22, 23, 24, 1, 2, 3, 4, 5, 21, 22, 23, 24, 1, 2, 3, 4, 5, 21, 22, 23, 24, 1, 2, 3, 4, 5, 21, 22, 23, 24, 1, 2, 3, 4, 5, 21, 22, 23, 24, 1, 2, 3, 4, 5, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 6, 7, 8, 9, 10, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 11, 12, 13, 14, 15, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24, 16, 17, 18, 19, 20, 21, 22, 23, 24 };

values = { 1.560381586921, 0.007801517859, 0.000888034535, 0.035617131505, 1.000000000000, -0.876518082605, 1.282062390000, -0.019696949887, -0.019483465963, 0.008891192058, -0.109906056502, 0.001545121088, 0.045101165387, 1.000000000000, 0.124340980968, 0.010989533912, -0.033401421519, 0.000587524190, 0.002254107738, -0.051972579269, 0.003425019422, 0.064334367916, 1.000000000000, 0.042358852907, -0.039536523557, -0.030121979059, 0.052692523271, 0.001016243192, -0.032763315112, 0.013471489319, 0.121670743120, 1.000000000000, -0.016766350272, -0.030586529961, 0.055349274383, 0.107017826699, 0.000580904145, -0.023293925809, 1.020249499141, 0.005100992446, 1.000000000000, -0.014943391720, -0.011267140561, 0.481355569107, -0.433325394230, -0.007801517859, -1.560303571743, 0.023279990067, -0.000813437866, 0.000000000000, 1.290955777065, 0.889338706505, -0.012440921639, 0.029826764466, -0.124540671940, -0.008822560758, 0.033377431407, -0.001477758196, 0.000000000000, 0.011346925356, -0.109710569126, 0.000886209308, 0.045121543677, -0.065751959359, -0.002193263722, 0.053165934599, -0.003363886945, 0.000000000000, -0.050340130519, -0.033805660857, 0.043943301386, 0.037063791237, -0.045781897791, -0.000961932540, 0.111007386583, -0.013415761224, 0.000000000000, -0.042613359144, 0.011782241875, 0.097162455357, -0.059421326274, -0.035631325462, -0.000532128379, -0.005100992446, -1.020198489217, 0.000000000000, -0.017135127234, 0.009697009437, -0.438182282460, -0.485688823049, 1.560381586921, 0.007801517859, 0.000888034535, 0.035617131505, 1.000000000000, 0.000737568666, 0.000506329867, -0.000007081363, 0.000017040360, 0.008891192058, -0.109906056502, 0.001545121088, 0.045101165387, 1.000000000000, 0.000000065466, -0.000000633833, 0.000000005144, 0.000000260678, 0.002254107738, -0.051972579269, 0.003425019422, 0.064334367916, 1.000000000000, -0.000000768036, -0.000000515588, 0.000000670422, 0.000000565254, 0.001016243192, -0.032763315112, 0.013471489319, 0.121670743120, 1.000000000000, -0.000072264491, 0.000020055686, 0.000164733289, -0.000101045958, 0.000580904145, -0.023293925809, 1.020249499141, 0.005100992446, 1.000000000000, 0.020769765333, -0.011636781137, 0.526000389988, 0.588769249632, -0.007801517859, -1.560303571743, 0.023279990067, -0.000813437866, 0.000000000000, 0.000499004813, -0.000732505367, 0.000011252683, 0.000011090937, -0.124540671940, -0.008822560758, 0.033377431407, -0.001477758196, 0.000000000000, -0.000000718357, -0.000000063411, 0.000000192968, -0.000000003426, -0.065751959359, -0.002193263722, 0.053165934599, -0.003363886945, 0.000000000000, -0.000000646035, 0.000000603208, 0.000000459382, -0.000000803907, -0.045781897791, -0.000961932540, 0.111007386583, -0.013415761224, 0.000000000000, 0.000028537670, 0.000051869719, -0.000094116054, -0.000181445788, -0.035631325462, -0.000532128379, -0.005100992446, -1.020198489217, 0.000000000000, -0.017933993153, -0.013656317825, 0.583568122657, -0.520112697492, 1.560381586921, 0.007801517859, 0.000888034535, 0.035617131505, 1.000000000000, 0.000737568666, 0.000506329867, -0.000007081363, 0.000017040360, 0.008891192058, -0.109906056502, 0.001545121088, 0.045101165387, 1.000000000000, 0.000000065466, -0.000000633833, 0.000000005144, 0.000000260678, 0.002254107738, -0.051972579269, 0.003425019422, 0.064334367916, 1.000000000000, -0.000000768036, -0.000000515588, 0.000000670422, 0.000000565254, 0.001016243192, -0.032763315112, 0.013471489319, 0.121670743120, 1.000000000000, -0.000072264491, 0.000020055686, 0.000164733289, -0.000101045958, 0.000580904145, -0.023293925809, 1.020249499141, 0.005100992446, 1.000000000000, 0.020769765333, -0.011636781137, 0.526000389988, 0.588769249632, -0.007801517859, -1.560303571743, 0.023279990067, -0.000813437866, 0.000000000000, 0.000499004813, -0.000732505367, 0.000011252683, 0.000011090937, -0.124540671940, -0.008822560758, 0.033377431407, -0.001477758196, 0.000000000000, -0.000000718357, -0.000000063411, 0.000000192968, -0.000000003426, -0.065751959359, -0.002193263722, 0.053165934599, -0.003363886945, 0.000000000000, -0.000000646035, 0.000000603208, 0.000000459382, -0.000000803907, -0.045781897791, -0.000961932540, 0.111007386583, -0.013415761224, 0.000000000000, 0.000028537670, 0.000051869719, -0.000094116054, -0.000181445788, -0.035631325462, -0.000532128379, -0.005100992446, -1.020198489217, 0.000000000000, -0.017933993153, -0.013656317825, 0.583568122657, -0.520112697492, 1.560381586921, 0.007801517859, 0.000888034535, 0.035617131505, 1.000000000000, -0.872241323662, 1.284998321104, -0.019738010855, -0.019384658199, 0.008891192058, -0.109906056502, 0.001545121088, 0.045101165387, 1.000000000000, 0.124343772841, 0.010962503169, -0.033401202162, 0.000598641185, 0.002254107738, -0.051972579269, 0.003425019422, 0.064334367916, 1.000000000000, 0.042341057997, -0.039548469393, -0.030106445815, 0.052705619842, 0.001016243192, -0.032763315112, 0.013471489319, 0.121670743120, 1.000000000000, -0.016873473173, -0.030556799962, 0.055593470495, 0.106868038940, 0.000580904145, -0.023293925809, 1.020249499141, 0.005100992446, 1.000000000000, -0.014774064958, -0.011362010120, 0.485643818954, -0.428525417475, -0.007801517859, -1.560303571743, 0.023279990067, -0.000813437866, 0.000000000000, 1.293849234172, 0.885091306873, -0.012375673459, 0.029891074770, -0.124540671940, -0.008822560758, 0.033377431407, -0.001477758196, 0.000000000000, 0.011316289972, -0.109713273390, 0.000894438699, 0.045121397554, -0.065751959359, -0.002193263722, 0.053165934599, -0.003363886945, 0.000000000000, -0.050355098728, -0.033791684919, 0.043953944983, 0.037045165226, -0.045781897791, -0.000961932540, 0.111007386583, -0.013415761224, 0.000000000000, -0.042571055685, 0.011859132124, 0.097022940298, -0.059690296537, -0.035631325462, -0.000532128379, -0.005100992446, -1.020198489217, 0.000000000000, -0.017281335193, 0.009585675483, -0.433424708206, -0.489929073129 };

[adlam]

Hi there, 

Just wondering if is there any follow up on this issue?

Regards,

[Gennady_F_Intel]

the fix of this issue will be added to the one of the next release. We will let you know when this release will be available.

[Gennady_F_Intel]

Hi, please check the latest MKL 2019 u3 and let us know if the problem is still there.