segmentation fault when mkl_sparse_d_mv for sparse matrix with many nonzeros

Hi,

I am computing y=Ax (mkl_sparse_d_mv) for a 81,000 x 81,000 sparse matrix A (csr format) with 2.3*10^9 nnzs, but MKL produce segmentation fault.

The sparse matrix is dumped from another program. I first read it from disk, and then call mkl_sparse_d_create_csr to create the matrix. mkl_sparse_d_create_csr returns SPARSE_STATUS_SUCCESS. But when I move to mkl_sparse_d_mv, it shows me segmentation fault.

My another observation is that when I decrease #nnz to be 1.5*10^9, MKL can finish the computation successfully.

My hypothesis is that MKL 2.3*10^9 > 2^31, and MKL may have constraints on #nnz in a sparse matrix.

Any help would be appreciated. Thanks!

Env: Linux, parallel_studio_xe_2019_update4.